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Построение и исследование корректности математической модели транспорта и 


осаждения взвесей с учетом изменения рельефа дна" 


А. И. Сухинов', В. В. Сидорякина””" 


' Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация 
* Таганрогский институт имени А. П. Чехова (филиал) РГЭУ (РИНХ), г. Таганрог, Российская Федерация 


Пеу@ортепЕ апд соггес@те$$ апа1у$1$ оЁГ Фе та Фетайса|! плоде] оЁ гапзрогЕ апд зизрепуоп зедппетайоп 


дерепате оп Бобот геНеЁ уатаНоп "^^" 


А. Г. ЗаКБштот', У. У. ЗаогуаКта”” 


'роп б{ае Тесбшса! Оптуегзиу, Козюу-оп-Роп, Киззлап Недегайоп 


* Тазапгое СвеКвоу №ш$Н иже, Возюу Зе ОшуегзИу оЁ Есопопус$ (В ПМН) БгапсВ, Тазапгоэ, Виз$1ап Недеганоп 


Введение. Настоящая работа посвящена исследованию 
пространственно-трехмерной модели транспорта и оса- 
ждения взвеси в прибрежной зоне с учетом изменения 
рельефа дна. Модель учитывает следующие процессы: 
адвективный перенос, обусловленный движением водной 
среды, микротурбулентную диффузию и гравитационное 
осаждение частиц взвеси, а также изменение геометрии 
дна, вызванное осаждением частиц взвеси или подъемом 
частиц донных отложений. 

Целью работы являлось проведение аналитического ис- 
следования корректности начально-краевой задачи, соот- 
ветствующей построенной модели. 

Материалы и методы. Изменение рельефа дна приводит к 
необходимости решать начально-краевую задачу для 
уравнения параболического типа с младшими производ- 
ными в области, геометрия которой зависит от искомой 
функции решения, что приводит, в общем случае, к нели- 
нейной постановке задачи. Выполнена линеаризация мо- 
дели на временной сетке за счет «замораживания» релье- 
фа дна в пределах одного шага по времени и последующе- 
го пересчета функции поверхности дна на основе изме- 
нившейся функции концентрации взвешенного вещества, 
а также возможного изменения вектора скорости движе- 
ния водной среды. 

Результаты исследования. Для линеаризованной задачи 
построен квадратичный функционал и энергетическим 
методом доказана единственность решения соответству- 
ющей начально-краевой задачи в пределах произвольного 
шага по времени. На основе преобразования квадратично- 
го функционала получена априорная оценка нормы реше- 
ния в функциональном пространстве Г. в зависимости от 
интегральных оценок по времени правой части, гранич- 
ных условий и начального условия, и, таким образом, 
доказана устойчивость решения исходной задачи при из- 


тго4исноп. ТВе рарег 15 4еуое4 ю Фе зшау оп Фе @тгее- 
Читепзлопа! плоде] о? шапзрог апа зизреп$1юп зедппешайоп 
ш Фе соа$а| агеа 4ие ю срБапоез ш Фе Боцот гейеГ. ТБе 
то4е] соп$14дегз Фе ГоПо\лше ргосез$ез: адуесиуе тапзЁег 
саизе4 Бу Ше адцайс тедтлт тойоп, писго-вагщет а - 
$1оп, ап огаупу зедптещаНоп оР зизрепае4 рагас]ез, аз мей 
аз фе Боцот соеотешу уапайоп саизеа Бу Фе рагасе зе{- 
(пох ог Бобошт зедитепЕ г1$11>. ТБе у’огК об]есйуе \аз ю 
соп4исЕ ап апа!уйса| заду оГ Фе согтестез$ оГ Фе шиаа|- 
Боип4агу уаше ргоет согезроп4аше № Ше сопзигацщеа 
штоб]. 

Маепа[5 апа Мейодфб. 'ТБе сБапхе ш фе Боцот гейеТ а145 
Ш зошйоп {о Фе пийа1-Боипдагу уаше ргоет Гог а рагаБо]- 
1с едиайоп уф Фе 1о\езё депуайуез ш а 4оташ \Возе 
сеотехшу 4ереп4$ оп Фе 4езиед Гапсйоп оЁГ Фе зомйоп, 
\ ср т сепега| ]еа4$ © а попппеаг Гогпу]айоп оГ Фе ргоБ- 
ет. ТЬе то4е] 15 Ппеатте4 оп Фе ите о114 дие ю Фе “Ёее7- 
112” ог Фе Бобот гейе? уп а зае $ер ш Ите ап Фе 
заб зеаиепЕ гесаси]аНоп о? фе Боцот зигЁасе Гапсйоп оп Фе 
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менении начального и граничных условий, функции пра- 
вой части. 

Обсуждение и заключения. Модель может представлять 
ценность при прогнозе распространения загрязнений и 
изменения рельефа дна, как при антропогенном воздей- 
ствии, так и в силу естественно протекающих природных 
процессов в прибрежной зоне. 


Ключевые слова: прибрежные системы, математическая 
модель, задачи диффузии-конвекции осаждения взвешен- 
ного вещества, изменение рельефа дна, единственность 
решения и устойчивость начально-краевой задачи. 
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Введение. Среди факторов, определяющих комплексное развитие исследований прибрежных террито- 
рий, большое значение принадлежит охране водной среды [1-2]. Устранение последствий естественных про- 
цессов, таких как загрязнение, заиление и истощение акваторий приводит к необходимости изучения всех ас- 
пектов, влияющих на изменение прибрежных акваторий. Поддержание водных объектов в надлежащем состоя- 
нии и своевременное вмешательство в режим его функционирования напрямую связано с увеличением порто- 
вых мощностей и обеспечением эффективного развития прибрежной инфраструктуры (обеспечение подхода к 
причалам судов с низкой посадкой, очистка от ила и водной растительности прибрежной полосы и пр.) [3—3]. 
Как правило, проведение исследований в данной области требует построения математических моделей, макси- 
мально приближенных к реальным процессам [6-11]. 

В работе рассматривается непрерывная математическая модель, описывающая пространственно- 
трехмерные процессы, связанные с транспортом и гравитационным осаждением взвесей в водной среде с изме- 
няющимся рельефом дна. Эта модель учитывает микротурбулентую диффузию и адвективный перенос взвесей, 
действие на частицы взвеси силы тяжести, наличие дна и свободной поверхности, изменение рельефа дна. Мо- 
дель транспорта взвеси позволяет проводить исследования гидрофизических процессов водных систем, осу- 
ществлять прогноз динамики изменения донной поверхности на основе описания процессов подъема, переноса, 
осаждения, изменения концентрации взвеси [12—13]. Доказана единственность решения соответствующей 
начально-краевой задачи и получена априорная оценка нормы решения в зависимости от интегральных оценок 
правой части, граничных условий и начального условия. 

Материалы и методы. Непрерывная ЗО модель диффузии-конвекции взвеси и соответствующая 
ей начально-краевая задача. Рассмотрим непрерывную математическую модель распространения взвеси в 
водной среде, учитывающую диффузию и конвекцию взвеси, действие на взвеси силы тяжести, наличие дна и 
свободной поверхности. Будем использовать прямоугольную декартовую систему координат Охуз , где ось Ох 


проходит по невозмущенной водной поверхности и направлена в сторону моря, ось О2 направлена вертикаль- 
но вниз. Пусть А=Н+п — общая глубина акватории, м; Н— глубина при невозмущенной поверхности водоема, 
м; | — возвышение свободной поверхности относительно геоида (уровня моря), м (рис. 1). 





Свободная поверхность 


Рис. 1. Введение системы координат Охуз 
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Пусть в замыкании области С’= 10 <х<[,,0<у<Г,,0<&<Н(х, У) находятся частицы взвеси, ко- 


торые в точке (х, у, =) и в момент времени { имеют концентрацию с = с(х, У, 5,1), мг/л; { — временная пере- 


менная, сек. Также далее будем использовать обозначение Г, = шах Н (х, У) | 
0<х<Г,, 0<у<Г, 


Поведение частиц взвеси будет описано следующей системой уравнений: 


О[и-+и ]с 2 2 
бе, О(ие) _9(ис) , Э((т+н, с) _ ы а 9 м 
О Ох ду [бе 0х ду 0% 05 (1) 
ОН 5 
Е--ИС, 
0: р ° 


— 


где и, у, и — компоненты вектора И скорости движения жидкости, м/сек; у, — гидравлическая крупность или 
скорость осаждения частиц, м/сек; ци», и, — коэффициенты горизонтальной и вертикальной турбулентной диф- 
фузии частиц соответственно, м/сек; ЕЁ — мощность источников частиц; & — пористость донных материалов. 

Слагаемые в левой части (кроме производной по времени) первого уравнения системы (1) описывают 
адвективный перенос частиц, обусловленный инерционным движением водной среды, а также осаждение под 
действием силы тяжести. Слагаемые в правой части описывают диффузию взвеси. Коэффициент вертикальной 
диффузии выбран отличным от коэффициента горизонтальной диффузии в связи с тем, что эффект различия 
этих коэффициентов часто наблюдается в различных средах и может быть вызван различными факторами. 

В качестве области С рассматриваем «скошенный» к берегу «параллелепипед» АВСРАОС О, ‚, верх- 


нее основание АОСО, которого лежит на свободной поверхности (= = 0) ‚ анижним основанием АВСО явля- 


ется часть поверхности дна (= = Н(х, У)). Пусть 5 — поверхность С, И — внешняя нормаль к поверхности 


к * 
«скошенного Параллелепипеда». Будем считать заданнои и — скорость движения жидкости на боковых по- 


верхностях С’. В совокупности с граничными условиями первого рода для функции концентрации частиц это 
позволяет определять потоки взвешенного вещества как по направлению к берегу, так и вдоль берега (рис. 2). 


Ст Берег 






Поток 
частиц етекает 


А 


Рис. 2. Область решения задачи транспорта взвеси 


Добавим к системе (1) начальные и граничные условия (предполагая, что осаждение частиц на дно не- 


обратимо). 
В качестве начальных условий при времени 1=0 принимаем 
с(х, у, 2,0) = с, (х,у,2); (2) 
Н(х,у,0)=Н, (х, у). (3) 


Граничные условия задаем на гранях АВСРАОС О, (задаем потоки взвесей как по направлению к бе- 


регу, так и вдоль берега): 
— на гранях 5 =ААОВ (х=0, О<у<1,, 0<:<1[,} , 5, ЕААБЬ (у=1,, О<х<(,, 0<5<[.) и 


5, =ВОСС (у=0, О<х<Д, 0<5<[) 
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С=С „где с. со [= [0,1]; (4) 


с=0; (5) 
— на поверхности 5, = АВСО (==Н (х,у,!), О<х<Ё, О<у<Г ) 


— на гранях 5, =РБСС (х=(,, 0<у<1,, 0<:<1,) и 5, =АОСЬ, (==0, О<х<1,, 0<у<[,) 


х 


и? и? 
т ИЛИ в (6) 
и Ц, 0% Ц, 








Следующее условие невырожденности области решения задаем для всех (х, У,Е) при которых постав- 


лена начально-краевая задача: 
Н(х, у,1) > ®, = соп5Р> 0  0<#ЕЗТ. (7) 


При исследовании объединенных моделей транспорта наносов и взвесей возможно увеличение концен- 
трации взвешенных частиц в придонном слое за счет подъема частиц донных отложений при превышении 
сдвигового напряжения некоторого критического значения [13-16]. Тогда вместо граничного условия (6) будем 
рассматривать граничное условие вида 
З- =@с, @= соп51 > 0. (8) 

0% 
Линеаризация начально-краевой задачи транспорта и осаждения взвесей. Для того чтобы создать 
линеаризованную модель на временном отрезке 0 <Е<Т построим равномерную сетку ©. с шагом т, т. е. мно- 


жество точек ©. = {1 = их, п = 0,1....,№М, № =Т}. 





о х 


Рис. 3. Построение временной сетки 


Функции с“ (ху) иН (") (х, у, „} определяем на каждом шаге временной сетки ©... Если п=1, 
(1) ни 
то в качестве с |х,у, <, ), х, У, достаточно ВЗЯТЬ функции начального условия, т. е. 


с) (х, у, 5,0) =. (х, у, 5) , Н® [о У, ) =Н) (х, У) соответственно. Если же п=2.....№ ‚ то функции 


с) ( у, т) = с" [5 у аа предполагаются известными, поскольку предполагается решенной задача 
(1)-(6) для предыдущего временного промежутка @_, <<. 


Систему (1) на промежутке Е << запишем в виде: 














дс) 9 (ис") | д (ег д ((* +, ) с") | 22) дс) 5 ас") 
аи, + |+— и, +Р, 
(611 Ох ду [бе Ох ду 05 {ба (9) 
ан“) [3 (п) 
=--тс 
0: р ° 
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и дополним ее начальными условиями: 
с о 2), с” = с" и 
Н® а У, ел [х у), Н“” (х, У) = Ни) [ь. У, 1), Яо. №. (11) 
Граничные условия (4)—(6) предполагаются выполненными для всех промежутков времени #_ <1<Е. 
Определив функцию с") [5 и =с" 9 (х, Ут) на временном промежутке (_ <#<@, можно 
найти функцию Н (и) (х, А) С этой целью проинтегрируем обе части второго уравнения системы (9) по пе- 


ременной {_, <1<[. Получим 





ГА 1 
‚ ОН“ о 
| @=-2 и, са. (12) 
Г 0 Р г 
Из равенства (12) нетрудно получить 
М т 
Н®) = Н® Пи, [о (13) 
п= а. - 
Введем на каждом временном шаге ГРЫ Е, область 


п 
С ={0<х<1,,0< у<[,,0<&<Н (х, у,1,,)} 
Имеем цепочку линейных начально-краевых задач для каждого временного слоя, где для промежутка 


1 <#<[ рассматривается система вида 








дс) (ис) 9(»с°) 0((»+ и, }е® 920) ‚92 п) р. ас(") 
АА — Е, 
а аа м пы а |* 4) 
(х, у, 2 ) ЕС, ии О<у<Ё,, 0<&<Н"" (х, у, р 
Н® =Н(9 3 [еар, п=1ь2,.М. 
ы (15) 
с начальными условиями: 
с" | У, вы = Г (х, у, 1 ). (16) 
н“) 1 У) — Ни (х, Ут) | (17) 


Отметим, что на каждом временном шаге граничные поверхности будут меняться (исключение грань 


5. ). Рассматривая временной промежуток #1_ <<, ‚ осуществляем задание граничных условий на гранях об- 


ласти С, _ 


— на гранях 9. „.(х=0, 0<у<1,, 0<=<Н“" 9 (0,у,1,,)} $. (У=Ё, 05х51, 052$ НН" (ху, )} 


и 5. (У=0, О<х<Г,, 0<5<Н"' (х,0,1, 


и 
с") =с где с. =с (х,у, 1), ей |; (18) 
— на гранях 
5.1 (х=1, О<У<Е,, 0<=<Н® (1,1, )) 
и 5, (5=0, О<х<Г,, О<у<Е,)=АОСВ, 
с 0 (19) 
— на поверхности 5 (- =Н"” (ху), 0<х<Ё,, 0<У<Е, 











Ее не ыы (20) 
дп и, 0% и, 
Граничное условие (8) заменится на следующее 
(") 
— = ас”, а = соп54 > 0. (21) 
0% 
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Таким образом, предполагается, что рельеф дна в пределах данного временного шага при расчете рас- 
пределения концентраций взвешенного вещества не меняется и берется из предыдущего временного слоя. В 


первую очередь на данном временном шаге 1_ <#1<Е решается начально-краевая задача для уравнения кон- 
> > —1 
векции-диффузии (14) с зафиксированной функцией рельефа дна Н (9 


) 


единственности иИ непрерывной зависимости решения от входных данных задачи проводится на фиксированном 


‚ а уже затем выполняется обновление 


(пересчет) функции рельефа Н () в соответствии с равенством (15). Определение условий существования, 
временном слое в этих предположениях и при условии выполнения условия (7). 

Авторы не планируют исследовать в данной работе существование решений начально-краевых задач 
(14)—(20) и (14)—(19), (21). Вопросы существования решений начально-краевых задач для уравнений параболи- 
ческого типа с младшими производными (уравнений диффузии-конвекции) рассмотрены, например, в моно- 
графиях [17—18]. 

Результаты исследования. Исследование единственности решения начально-краевой задачи 
транспорта взвеси. 

Рассмотрим начально-краевую задачу (14)-(20), сформулированную для произвольного временного 
слоя (1 <{<6. 


Умножим левую и правую часть уравнения (14) на функцию с) и получим: 


2" + ый и + (=, =") — се Е р | с") 2 =” сЕ .(22) 
х . 


О Ох ду {ба 05 
Левая часть равенства (22) может быть преобразована следующим образом: 
2 

д с | 








. (") (") (") 
п ее но 
01 Ох ду 0% 2 (23) 
9[с(®) 
= г +4 (<®) я 
где (И = |“. у, жи |. 
С учетом (23) уравнение (22) запишется в виде 
2 
ое” а _ [52° 29\ са[ &®\ © 
(а | ов им '8 м р. (24) 


Затем обе части уравнения (24) проинтегрируем по времени на интервале #_ <1<Е и, после этого — 
по пространственным переменным в области С, . В первом слагаемом порядок интегрирования изменён в си- 


лу теоремы Фубини [19]. Получим 


О (”) . р. — 
Г. ] р Ве. | (“бис ‚_ 


п-1 


р . 92с®  62е“) т, „д АЕ 
-1| еее в, -] ее. ‚+ 5 


Г [сэ ас, , 
С 


1—1 


Первое слагаемое из левой части равенства (25), очевидно, равно 


и 


(„) \* 
> ВыИЕ С. = ПТ”) се.) (6%) без, ав. (26) 


—1 


Далее обратимся к преобразованию второго слагаемого левой части равенства (25). С учетом формулы 
Остроградского-Гаусса и граничных условий (18)—(20) его можно записать в виде [20]: 
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| Гао свв. РЕЛЕ } (© пе Поль а 


—1 


(С) (и по Де де --] (< ви (7) 
Ил 


(с . п [ | Г (с ран [ | |] ей )) к | 
где И’ — известная скорость водной среды на гранях, где заданы граничные условия первого рода; фактически 


И 5б, п-1 


это все боковые грани, кроме 5,,„, и верхней крышки 5. ‚, на которых концентрация взвеси равна нулю, по- 


тому и потоки через них равны нулю. 
Обратимся к преобразованию правой части равенства (25). Имеет место равенство 


п) (п) (п) 
= в, |6, - 
ах. ду| ду 0 0% 
дс" 9[ дс” |] 9 дс” 
-] в в + | 9,7 | аб, - (28) 
Ох Ох ду ду 05 05 
дс") | дс") | дс") 
-[| И | ТИ || ТЫ, ас. 
Г. Ох ду 0% 


ИИ „д (п) „д (п) „ а (п) 
Пусть О = 10,.0,,0. | Не Ирис к ‚с ут —— . Тогда в силу теоремы Остроградского- 
х и 























Гаусса имеем: 











1] и, 2 8 нь 2 и. |+ > ты я | 5. = 14» дас = т 


= | Оха + 1 О, дас + Л О, ха + 1 О, дас + И Оку + ‚Гомь = (29) 
и -1 О. ака + И О. ака + Л оду Л Оф. = 


2,п-1 3,п-1 


Преобразуя каждое слагаемое из правой части (29) с учетом условий на границе (18)—(20), получаем 


Я (”) (п) 
1] = п с, 9 |6, = 
Ох Ох ду ду 0% 05 (30) 


-П м. И сы + р пы }» м ›) АхаУ. 











С учетом (06), (28), (29) и (30) равенство (25) принимает вид 


ео бозныаа, ле (изень > со 
Ил [Е Ус", ети = пе ш 
дем 


(31) 











‚| 2" | дс") | 
5 +ы, : А(,_, | 4 = 
у с 


=" У (аи, р 


Тождество (31) будет основным при исследовании единственности и получения априорной оценки 
нормы решения начально-краевой задачи (14)—(20). В случае замены граничного условия (20) на граничное 
условие (21), квадратичный функционал (31) изменяется следующим образом: 
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Песен | [ Д(Неуынсьда- 


п-1 —1 51 1 


| [5 УСН, Ни Ту рее 
"Л (ви, 9 раша | 0“ "|. „(“”] "|=. 46. а - 
Е“) боны раб + | ОИС, 


Предположим, что уравнению (14) с одними и теми же условиями (16)—(20) удовлетворяют два различ- 


(32) 











ных решения задачи с! = с, (х, у, #1 с, = С, (х, у, 2,1) . Для их разности с=с —с, справедлива следующая 


начально-краевая задача: 


БР 9(ис) д(ус) 0((»+ и, )с) 92е 6е\ д дс 
И, (33) 
О Ох ду 05 Ох ду 05 05 
с(х, у, 5,0) = 0, (х, у, 5) ЕС. р (34) 
— на гранях 5 2 5 и 7 9 7 и ? 95 
Сс =6 =0: (35) 
— на поверхности 55 „| 
и? и? 
[а =-—(с-с,)=--*&. (36) 
в, м, 


Для функции с равенство (33) примет вид с учетом равенств (34)—(36) 


52 ху, 51,) 46, За [» даа 


, | [0 (®. “(®) (м. #0. (7) 


Поскольку и’, > 0 и другие известные величины под знаком интегралов положительны [„>?0, и,>0 то 


равенство (36) выполняется лишь при условии 
<(х, у,2,1) =0, (х,у,2)еС., В. <ЕЗЬ, (38) 


что и завершает доказательство единственности решения начально-краевой задачи (14)—(20). 
В случае замены граничного условия (20) на соотношение (21) получаем вместо выражения (37) равен- 


ство вида 


52 х, у, <, ) АВ... [1 [ль Ро 


м. — — . (39) 
[> +Ы, — +1 “ АС, | 4 =0. 
Ох ду 05 


Потребуем выполнения неравенства 


1 
г. —“щ, >20, (х, у, 5) = 6—1 к Ро 


п 


ИЛИ 
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и 
‚ (хуебнаь В. << (40) 
ы, 


тогда все слагаемые в равенстве (39) неотрицательны и равенство нулю возможно тогда и только тогда, когда 


п? 


С (х, у, 2,1) = (), Ге у, 5) ЕС, Е <Е<Ь, что и означает единственность решения и в этом случае. 
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Рассуждения аналогичным образом повторяются для всех слоев временной сетки ‹.. Изменение гра- 
ничных условий, связанное с непрерывным изменением рельефа дна в зависимости от временной переменной 
требует дополнительного исследования и выходит за рамки данной статьи. 

Теорема. Пусть дана система уравнений 

дс) 9(ис®) 9(»с") 0((»+ и, ти 525) 229) а ас") 
к | СО Е ЗОИЕНИИ 

0: дх ду 0% дх” ду” 05| ^ д 
(х,у,=)} ЕС, С ={0< хер, 0<у<Ь, 05 «НР (хуи, 


п 








+Р, 


м & 
Н® = НН а >. | саг, п=1,2..... № 
р 


1 


1 
ГА 

8 
пП=| 3 


в односвязной области ОФ _=С х(Г<1<Ё)}, а =0<х<Ё,0<у<Ё .0<5<Н"% х, у )], Сс ДО- 
п—1 п-| п-| п 1 х у п-| 


п 
“ > —1 
статочно гладкой границей, определяемой гладкостью функции &=Н (9 ра У), О<х<Ё,, О<у<Ё, с 


начальными и граничными условиями (16)—(20). Пусть функции решения с (х, у, а вектора скорости 


Т (п-1) $ 
иИ,У, и- Я. | ‚ Начального УСЛОВИЯ С [% у, ь | ) 5 правои части К (х, у, о Г) ь граничного усло- 





водной среды 


вия С” (х, у, 5,!), коэффициента вертикального турбулентного обмена и, =н, (=), я у, =) еС._ удовлетворяют 








следующим условиям гладкости:  с”(х,у,21,.)еС*(9,,)ОС(О,.). та с" еС(О,_). 
[музу | еС'(@,)оС(О,_), с" (х,у, и.) = С(б,.) , Е(ху,=,/)еС(0,) , 
и, (х, у,=)=С' (6) с С(б,.) с" (х, у, 2,1) Е С(5, [<< |, бы 
ес <, ОЗУ, =" (ху) {ь а 515], алаюже условиям согласованности граничных 
и начальных условий с’(х,у, 5,0) =с, (х,у,=), (х,у, <) Е 5 \(0 вой беузь, вн” )). 
ы =- й с", (0<х<1,0<у<ЕЁ, ==Н"\(жу)), тогда решение данной задачи существует и един- 
ственно. | 


Замечание. В случае замены граничного условия (20) на граничное условие (21), в качестве достаточно- 
го условия выполнения предыдущей теоремы следует добавить неравенство (40). 


Исследование непрерывной зависимости решений начально-краевой задачи транспорта взвеси 
от начальных, граничных условий и функции правой части. Следующий этап связан с исследованием не- 
прерывной зависимости решения от функций правой части, граничных и начальных условий для системы (14) 
(15). 

Будем предполагать, что 

с > ео = сои > 0, 
О 


п-1 — 


(41) 


Введем для удобства обозначения: объединение всех частей боковой цилиндрической поверхности 


(границы области С„./) обозначим как 5 а нижнее основание области С -— 5,,,. В силу условий глад- 


с.п-1? 
КОСТИ, перечисленных В условиях вышеупомянутой теоремы, достигаются экстремумы функций На ограничен- 


ных замкнутых множествах: 

















(") (") 
М = пах |с° } М, = тах — — 
- т || @х || ду (42) 
Му = пах, }, М — ‚тах, м ? у }, М — ини, } 


Будем ориентироваться на уравнение (31), если используется граничное условие (20), и равенство (32) 
— в случае граничного условия (21). Привлекая неравенство Фридрихса, имеем цепочку неравенств: 
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ас") ас) у ас) у 
С С С 
—— ас. > 
И а _. ду — д ой 
с) (ас [Га 
> пиши ,, ас. > 43 
орт (р |] мы № (43) 


1 1 1 з 
Витино] 06 6. 
Обратимся к уравнению (26), из которого в силу (42) и (43) получим неравенство: 
у 11 1 | | 
ПГ ') а +2М 5 п т? т ге С | у й Ум ›) АС, + 
з| [| В }* Ахау ых Ми | Ле Е (44) 


Ш 5ь, п-[ ба 


+2М 





›, М 3,п-1 





гам, | Ге. р 


Из неравенства (44) следуют два неравенства 


ем ас < регао нм. [ ев 


(45) 
+2 М „Мэ ГДЕ "45, нам. ПЛАС, й 
И 
ем) ав, «м. ДТезав,. + 
‚, (46) 
М РДС ы нам... | Пе аз и +2м „| | [[Е|4в,, & |. 

бел Аб 

= 
где р п” р 


Из полученных неравенств следует непрерывная зависимость (устойчивость) решения задачи (14)—(20) 


от функций: начального условия, граничных условий и правой части, в норме Г, для любого момента времени 
0 < 7 <-к®, а также в интегральной по времени норме [.. 


Очевидно, что при выполнении неравенства (45) и условий теоремы начально-краевая задача (14)—(19), 
(20) также будет иметь решение, непрерывным образом зависящее от функций: начального условия, граничных 
условий и правой части в соответствующих нормах. 

Обсуждение и заключения. Новизна предложенной нестационарной пространственно-трехмерной ма- 
тематической модели транспорта взвесей состоит в том, что наряду с учетом процессов адвективного переноса, 
микротурбулентной диффузии и гравитационного осаждения частиц взвеси модель описывает изменение гео- 
метрии дна, вызванное осаждением частиц взвеси или подъемом частиц донных отложений. Выполнена линеа- 
ризация соответствующей начально-краевой задачи на временной сетке и получены для произвольного времен- 


ного шага Г_ << условия единственности решения начально-краевой задачи и непрерывной зависимости 


от входных данных — от функций начального условия, граничных условий и правой части в норме гильбертова 


пространства Г. в интегральной по времени норме Г, для двух вариантов граничных условий. 
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Введение. Верифицирована методика численного модели- 
рования поперечного обтекания пролетных строений мо- 
стов на основе нестационарного решения Рейнольдса для 
уравнений Навье — Стокса (ОКАМ$, Опуеа4у Кеупо!&$- 
ауегаге4 Ма\мег — зюКез). Данный двумерный подход 
используется в современных методах и пакетах приклад- 
ных программ вычислительной гидроаэродинамики. Цели 
работы — отладка и экспериментальное обоснование ука- 
занной методики. Для реализации поставленной цели ис- 
пользована ранее разработанная авторами база данных по 
аэродинамическим характеристикам поперечных сечений 
пролетных строений балочных мостов типовых форм. 
Материалы и методы. Проведено численное моделирова- 
ние поперечного обтекания мостовых строений низкотур- 
булентными (гладкими) и турбулентными воздушными 
потоками в диапазоне практически интересных углов ата- 
ки. В итоге использовалась модель турбулентности Э5Т 
К — «. Методика предварительно отработана на тестовой 
задаче обтекания балок прямоугольного поперечного сече- 
ния. Расчеты проводились с помощью лицензионного про- 
граммного комплекса АМЗУЪ. 

Результаты исследования. Показано, каким образом угол 
атаки определяет силы (подъемную и лобового сопротив- 
ления) и момент поперечных сечений балочных мостов 
типовых форм. Полученные расчетные зависимости отно- 
сятся к пролетным строениям на стадиях монтажа (без 
плиты перекрытия и ограждений, без ограждений) и экс- 
плуатации в условиях модельных набегающих потоков — 
гладкого и турбулентного. Это позволяет очертить грани- 
цы для более взвешенных оценок аэродинамических ха- 
рактеристик балочных мостов в реальном ветровом потоке. 
Лобовое сопротивление сечению демонстрирует наилуч- 
шее согласование с опытными данными. Величина подъ- 
емной силы более чувствительна к наличию и протяжен- 
ности отрывных зон, поэтому ее расчетное определение 
менее точно. Наиболее проблемным для большинства 
конфигураций является воспроизведение влияния угла 
атаки на аэродинамический момент сечения. 
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Обсуждение и заключения. Сравнение расчетных и опыт- 
ных данных свидетельствует о применимости ОКАХ$- 
подхода к оперативному предсказанию аэродинамических 
характеристик однобалочных пролетных строений. Если 
же речь идет о многобалочных пролетных строениях, су- 
щественную роль играет аэродинамическая интерференция 
между отдельными балками. В этом случае вместо 
ОВАЮ-подхода следует применять более точные вихре- 
разрешающие методы. Полученные результаты могут быть 
использованы в исследованиях аэродинамики сооружений 
и в практике проектных организаций в сфере транспортно- 
го строительства. 


Ключевые слова: механика жидкости, газа и плазмы; 
математическое моделирование; вычислительная гидро- 
аэродинамика; ОКАМ№-подход; пролетные строения мо- 
стов; аэродинамические характеристики. 
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Введение. Известно, что обтекание ветром инженерных сооружений носит, как правило, нестационар- 
ный турбулентный характер, в потоке наблюдаются разномасштабные вихревые структуры [1, 2]. Вблизи пло- 
хообтекаемых тел (к которым, в частности, относятся пролетные строения мостов) возникают нестационарные 
отрывные зоны [3]. В связи с этим адекватное моделирование эффектов турбулентности становится важным 
требованием к методике вычислительного эксперимента. 

Для оценки аэродинамики сооружений используется моделирование крупных вихрей (Гагое Еа4у 
эпишаНоп, ГЕ) и гибридное описании потока (Реасвеа Едду Зипайоп, РЕЗ). Однако применение указанных 
методов осложнено их повышенной ресурсоемкостью, причины которой: 

— трехмерность задачи; 

— высокие требования к плотности расчетной сетки в пристеночной области и в области «фокуса» [4]; 

— ограничения на шаг интегрирования по времени; 

— относительно большая длительность временного «окна» для сбора нестационарной статистики в установив- 
шемся режиме. 

В то же время известно [5], что для цилиндрических призм, близких по форме к балочным мостовым 
пролетным строениям, двумерный подход воспроизводит основные свойства течения (первичная неустойчивая 
мода в следе за телом существенно двумерна). Так, в работе [6] выполнен РОЛ-анализ (Ргорег Огфогопа| Ое- 
сотроз1Яоп — истинное ортогональное разложение) течения вблизи призмы с относительной шириной сечения 
В/Н =5 (Н — высота сечения). В результате установлено, что 1-я и 2-я моды возмущений двумерны (постоян- 
ны по пролету) и соответствуют переносу вихрей вдоль поверхности. Трехмерные моды изменяются по пролету 


на характерной длине, которая не меньше ширины сечения Б. 


Для оперативного предсказания аэродинамических характеристик (АДХ) мостовых сооружений и со- 
провождения продувок в аэродинамической трубе авторы настоящей работы использовали нестационарное 
двумерное моделирование на основе осредненных по Рейнольдсу уравнений Навье — Стокса (ОКАХ-подход, 
Опзеаду Кеупо!4$-ауегагеа Мау1ег — 5юКез$). Его применимость к определению АДХ плохообтекаемых тел 


(неподвижных и колеблющихся) изучалась в ряде работ зарубежных авторов (см., например, [7]). 


Механика 


Материалы и методы. При постановке компьютерных экспериментов учитывались приведенные в [8— 


10] рекомендации. Расчеты проводились в программе Ншеп( комплекса АМЗУЪ. 


363 


Бер://уезбК.Аопза.га 


364 


Вестник Донского государственного технического университета. 2016. Т. 18, №4. С. 362-378. 155М№ 1992-5980 е155М 1992-6006 
Уезник о} оп Экие Тесйтса! Отуегзйу. 2018. Уч. 16, по. 4, рр. 362-378. 1585М 1992-5980 е155М 1992-6006 





Предварительно методика была отработана на тестовой задаче обтекания балок прямоугольного попе- 
речного сечения. В итоге для описания турбулентности потока была выбрана модель переноса сдвиговых 
напряжений (ЭВеаг 5гез5 Тгапзрой, э5Т) К — ®, подобраны параметры сетки и численного алгоритма. 

Протяженность прямоугольной расчетной области в продольном направлении — (30...40)Н, в попе- 
речном — (14... 20) Н. Передняя грань обтекаемого тела отстояла от входной границы на (8 ...12)Н. 

Использовались низкорейнольдсовые сетки (безразмерное расстояние первого узла до стенки у* < 
4...5), позволяющие рассчитывать отрыв и присоединение пограничного слоя. Ввиду сложности контуров об- 
текаемых тел строились многоблочные сетки. Внутренний, охватывающий тело, блок состоял из четырехуголь- 
ных элементов, плотность которых ближе к поверхности тела увеличивалась. Непосредственно у стенки созда- 
вался слой со структурированной ортогональной четырехугольной сеткой. Область следа покрывалась сеткой 
из квадратных ячеек размером не более Н/15...Н/10. К внешним границам размер ячеек возрастал до Н/ 
4...Н/З. На периметр сечения в зависимости от его формы приходилось порядка 10°...10? ячеек. Общее коли- 
чество ячеек варьировалось от 40—50 тысяч (для сечений простых форм) до 250-300 тысяч (для сложных). При- 


мер расчетной сетки возле балки трапециевидного сечения с перекрытием и ограждениями приведен на рис 1. 


Рис. 1. Пример расчетной сетки (фрагмент) 


При решении уравнений Навье — Стокса связь скорость — давление реализовывалась с помощью ал- 
горитма 5 1МРГЕ. Конвективные и вязкостные члены уравнений движения жидкости и переноса турбулентных 
параметров аппроксимировались схемами второго порядка точности. 

Численное интегрирование осуществлялось неявной по времени схемой 2-го порядка точности. Шаг 
интегрирования ДЕ составлял (0,02... 0,04) Н/У» (У, — скорость набегающего потока), т. е. при сходе вихрей с 
безразмерной частотой ДН /У., = 0,1 был приблизительно в 250—300 раз меньше периода 1/{, и это обеспечи- 
вало приемлемое разрешение нестационарных параметров потока. Установившаяся вихревая дорожка обычно 
формировалась к моменту (60...120)Н/ У. Таким образом, общее число шагов интегрирования составляло в 
среднем 6000-10000. Для сбора нестационарной статистики использовался отрезок времени, не меньший 5 пе- 
риодов. 

Пример качественного сопоставления расчетной и опытной картин течения вблизи пролетного строе- 


ния приведен на рис. 2. 
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Рис. 2. Пример картины обтекания пролетного строения: эксперимент (а), расчет (6) 


Результаты исследования. Подробные сведения по АДХ типовых поперечных сечений можно найти 
в [11]. На рис. 3-14 представлено сравнение расчетных и опытных данных по коэффициентам осредненных 
аэродинамических сил (лобового сопротивления, подъемной силы) и момента для некоторых характерных се- 
чений. 





Рис. 3. Коэффициенты осредненных аэродинамических сил. 
Здесь: Сха — лобовое сопротивление; Са — подъемная сила; Си — момент; 


Ви Н — продольный и поперечный размеры сечения (без учета ограждений); а — угол атаки 
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Набегающий гладкий поток характеризуется интенсивностью 0,5 %, турбулентный — 8 %. Расчетные 
данные представлены сплошными линиями. 
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Рис. 4. Узкая одиночная прямоугольная балка (соотношение В/Н 0,75): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 5. Узкая одиночная прямоугольная балка (соотношение В/Н 1,24): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 6. Узкая одиночная балка с перекрытием (соотношение В/Н 2,3): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 7. Узкая одиночная балка с перекрытием и ограждениями (соотношение В/Н 2,3): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 369 
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Рис. 8. Широкая одиночная прямоугольная балка (соотношение В/Н 1,85): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 9. Широкая одиночная балка с перекрытием (соотношение В/Н 3,9): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 37] 
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Рис. 10. Трапециевидная балка (соотношение В/Н 3,09): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 11. Трапециевидная балка с перекрытием (соотношение В/Н 5,1): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 


Механика 


ЭЙ 


Вестник Донского государственного технического университета. 2016. Т. 18, № 4. С. 362-378. 155М 1992-5980 е155М 1992-6006 
Уезник о} оп Экие Тесйтса! Отуегзйу. 2018. Уч. 16, по. 4, рр. 362-378. 1585М 1992-5980 е155М 1992-6006 











В/Н = 2/15 
а) 








= 
5 
77 08 
= У № 
= 0.6 
У 0.4 
о = 
|= я 0.2 
К ш р 
2 0 
> 0 - 475 |0) 
8 -0.4 ыы 
= ). 
в) 


Рис. 12. Двухбалочное строение (соотношение В/Н 2,75): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 13. Двухбалочное строениес перекрытием (соотношение В/Н 3,85): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 14. Многобалочное строение с перекрытием (соотношение В/Н 9): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Обсуждение и заключения. Анализируя результаты, можно отметить, что за некоторыми исключени- 
ями с увеличением относительной ширины В/Н сечения однобалочного строения точность расчетного предска- 
зания его АДХ несколько возрастает. Как правило, наилучшее согласование отмечается по лобовому сопротив- 
лению сечения. Для большинства конфигураций расчетные данные несколько выше коэффициента лобового 
сопротивления, полученного экспериментально. Поясним, что для рассматриваемых плохообтекаемых тел пре- 
обладающий вклад в лобовое сопротивление сечения вносит сопротивление формы (давления), в основном 
определяемое разностью сил давления на наветренной и заветренной сторонах сечения. Принятый расчетный 
подход огрубляет динамику вихревых структур в области за телом, что приводит к заниженному восстановле- 
нию давления в указанной области. 

Величина подъемной силы более чувствительна к наличию, протяженности и типу (открыта/замкнутая) 
отрывных зон. Особенно это касается оснащенного плитой пролетного строения — в этом случае возможно 
повторное присоединение потока к верхней стороне плиты с образованием замкнутой зоны отрыва (примерно 
при В/Н >5). Поэтому по сравнению с лобовым сопротивлением расчетное определение подъемной силы ме- 
нее точно, особенно для пролетных строений с плитой перекрытия. 

Наиболее проблемным для большинства конфигураций является воспроизведение влияния угла атаки 
на аэродинамический момент сечения. 

Если при поперечном обтекании многобалочных пролетных строений между балками возникает аэро- 
динамическая интерференция [12], точность предсказания АДХ падает с увеличением числа балок (относитель- 
ной габаритной ширины сечения). В этом случае вместо ОКА\$-подхода целесообразно использовать более 
точные вихреразрешающие методы РЕЗ и [Е$. 

Библиографический список 

1. Симиу, Э. Воздействие ветра на здания и сооружения / Э. Симиу, Р. Сканлан. — Москва: Стройиздат, 
1984. — 360 с. 

2. Ра4ои$$1$, М.-Р. Н/и19-згасваге ищегасноп$: сго$з-Йо\-шисеа шуаБШаез / М.-Р. Раё14ои515, 5.-. Рисе, Е. 4е 
Гапоте. — Ме\м УоК : Сатбпае Ошуегзиу Ргез$, 2011. — 402 р. 

3. Каугакоу, [. А зупеге1ас зва4у оРа СЕШО апа зеп-апа!уйса| плоде]$ ог аего@аз@с апа]уз1$ оЁ Без ш цифи- 
еп улиа сопа1воп$ / 1. КаугаКоу, (@. Могоепа! // оигпа] оЁ Е14$ апа эгасаге$. — 2018. — Уо|. 82. — Р. 59-85. 

4. Зраац, Р.-ВК. Уоцие регзоп’$ олае ю деасве4-е44у зпищайоп о4$ [Электронный ресурс] / Р.-К. Зраак ; 
МАЗА СВ-2001-211032 ; Воешо Соттегсла] Апр/апез. — Режим доступа: 
№ рз://пг5.паза.соу/агсыуе/паза/саз1.п5.паза.оу/20010080473.р4Е (дата обращения: 05.09.18). 

5. Мопкемиы, Р. А. АБзое шза Шу ш Фе пеаг-\уаКе оЁ (\о-4итепзюпа|! ЫюЁ Боез / Р. А. МопКезли, 
Г. М. Мепуеп // оигпа| оЁ 4$ апа 5гасвгтез. — 1987. — Уо1. 1, 155. 2. —Р. 165—184. 

6. Вгопо, Г.. Апа1уз1$ оЁ Фе зерагае По\ агоип4 а 5:1 тебапещаг суппаег @гоцеВ сотршайопа! зпишабоп 
[Электронный ресурс] / Г.. Вгапо, М. Со%е, О. Егапзо$ // ЕВ Еогореап & АЁлсап Сошегепсе оп УУта Епетеепипо : ргос. 
— Ногепсе, 2009. — Режим доступа: ВИр://\у\/\и ла\е.оге/Ргосее411$/5ЕАСУУЕ/1 15.р4Е (дата обращения: 05.09.18). 

7. Оп Фе еуаабоп оЁЬп4ее деск Ницег депуайуе$ изшо КАМ пигбщепсе тоде[5 / Е. Вгазлапи [е( а1.] // Хоигпа] оЁ 
У/та Епетееппс апд п4изла! Аегодупапис$. — 2013. — \Уо]. 119. — Р. 39-47. 

8. Кесоттепаайоп$ оп Фе зе о СЕО ш за епотееппе [Электронный ресурс] / 7. Егапке [её а1.|// Ве- 
зеагсВОае. — Режим доступа: 6р5://\у\у уу гезеагсВоае.пе/рибЙсаНоп/25'7762096_КВесоттепдайоп$ 
_оп_Фе_изе_оЁ_ СЕШ)_ш_рге@сипе_редезвлап_улид_епупоптейе (дата обращения: 05.09.18). 

9. Белостоцкий, А. М. Численное моделирование задач строительной аэродинамики. Разработка методик и 
исследования реальных объектов / А. М. Белостоцкий, С. И. Дубинский, И. Н. Афанасьева // Пиегпайопа! Лоигпа! оЁ 
Сотриабнопа| СУ апа эгасага| Епошеепие. — 2010. — Т. 6, № 1/2. — С. 67-69. 

10. Вычислительная аэродинамика строительных сооружений. Задачи и методы/ С. В. Гувернюк [и др.] // 
Вестник МГСУ. — 2011. — Т.2, №2. — С. 113-19. 

11. ОДМ 218.2.040—2014. Отраслевой дорожный методический документ. Методические рекомендации по 
оценке аэродинамических характеристик сечений пролетных строений мостов / Новосибир. гос. техн. ун-т ; Федераль- 
ное дорожное агентство (Росавтодор). — Москва : Изд-во ФГУП «Информавтодор», 2014. — 87 с. 

12. Саленко, С. Д. Аэродинамические исследования типовых многобалочных конструкций / С. Д. Саленко, 
Ю. А. Гостеев, А. Д. Обуховский // Теплофизика и аэромеханика. — 2013. — Т. 20, № 4. — С. 451-460. 


Поступила в редакцию 30.08.2018 Кесеуед 30.08 .2018 
Сдана в редакцию 03.09.2018 Бибииеа 03.09.2018 
Запланирована в номер 01.10.2018 эсредще4 11 Ше 1$54е 01.10.2018 


Механика 


577 


Бир://уезбиК.Аопза.га 


378 


Вестник Донского государственного технического университета. 2016. Т. 18, №4. С. 362-378. 155М№ 1992-5980 е155М 1992-6006 
Уезник о} оп Экие Тесйтса! Отуъегзйу. 2018. Уч. 16, по. 4, рр. 362-378. 1585М 1992-5980 е155М 1992-6006 





Об авторах: 


Гостеев Юрий Анатольевич, 

доцент кафедры «Аэрогидродинамика» 
Новосибирского государственного технического 
университета (РФ, 630073, г. Новосибирск, 

пр. К. Маркса, 20), 

кандидат физико-математических наук, доцент, 


ОКСО: Бр://огс14.0г>/0000-0002-4979-3174 
созбееу @ согр.пз@и.га 


Обуховский Александр Дмитриевич, 

доцент кафедры «Аэрогидродинамика» 
Новосибирского государственного технического 
университета (РФ, 630073, г. Новосибирск, 

пр. К. Маркса, 20), 

кандидат технических наук, доцент, 


ОКСТ:ВЕф://огс14.0г®/0000-0002-8842-0178, 
обихоуз КИ @ согр.пз@а.га 


Саленко Сергей Дмитриевич, 

заведующий кафедрой «Аэрогидродинамика» 
Новосибирского государственного технического 
университета (РФ, 630073, г. Новосибирск, 

пр. К. Маркса, 20), 

доктор технических наук, профессор, 


ОКСО: БЕф://огс14.0ге/0000-0003-8885-1885 
заепКо$4 @ уапдех.га 


Аи!шог5: 


Соз{ееу, Уигу А.., 

аззослае рго{еззог оЁ Фе Аегорудгод4упаис$ 
Рерагитет, МоуозизК 5 ае Тесбса! Отуегзиу (20, 
К. Магк$ Рг., Моуозфизк, 630073, ВЕ), Сапд.$са. 
(Рвуз.-Ма®.), аззослае рго{еззог, 


ОКСО: Бр://огс14.0г®/0000-0002-4979-3174; 
созбееу @ согр.пз@а.га 


Ороквоу$Юу, Мехапдег О., 

аззослае рго{еззог оЁ Фе Аегорудгод4упаис$ 
Рерагитет, МоуозизК 5 ае Тесбса! Отуегзиу (20, 
К. Магк$ Рг., Моуозризк, 630073, ВН), Сапа.Зс1. (Епе.), 
аззослае рго{еззог, 


ОКСО: Б@ф://огс14.0г$/0000-0002-8842-0178 
обихоуз КИ @ согр.пз@а.га 


заепКо, эегоеу О., 

Беа4 оГ Фе Аеговудго4упаптс$ Оерагаптеп, МоуозизК 
ое ТесЬтса| Фшуегзиу (20, К. МагК$ Рг., Моуоз1и К, 
630073, ВВ), Ог.5с1. (Еп.), рго{еззог, 


ОВКСТШО :Вф://огс14.0г</0000-0003-8885-1885 
заепКоз$Ч @ уапаех.га 


Вестник Донского государственного технического университета. 2016. Т. 18, № 4. С. 379-384. 155М№ 1992-5980 е155Мм 1992-6006 
Уезник о} Ооп Ма Тесйтса! Отуегзйу. 2018. У. 18, по. 4, рр. 379-384. 185М№ 1992-5980 е155М 1992-6006 





МАШИНОСТРОЕНИЕ И МАШИНОВЕДЕНИЕ 
МАСШМАЕ ВОПОИУС АХО МАСШУАЕ 5СТЕМСЕ 


УДК 62.593 





рз://401.0г2/10.23947/1992-5980-2018-18-4-379-384 


Влияние параметров рекуперативного объема на динамические характеристики 


пневмопривода при торможении' 


А. Н. Сиротенко", С. А. Партко”, Ваел Саллум”” 


1? Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация. 


Алеппо университет, г. Алеппо, Сирия 


ЕНесЕ оЁ гесирегайуе уоише рагатеег оп дупапис свагасег15с$ оЁ рпеитайс атуе ипаег БгаКшо 


А. М. Знобепко', 5. А. РагКо”, УУае! байопи?” 


12 роп 5ае Тесвлиса! Оптуегзиу, Козюу-оп-Роп, Ки$$1ап НедегаНоп. 


ЗОлмуегзИу оЁ АДерро. Аерро, Зупа 


Введение. Рассмотрены способы энергосбережения в 
пневматическом приводе. Интерес представляет способ 
торможения созданием противодавления в выхлопной 
полости пневмопривода. При этом энергия сжатого 
воздуха накапливается в рекуперативном объеме. Задавая 
начальные параметры рекуперативного объема, возможно 
управлять динамикой торможения. Цель работы — 
создание — математической модели, — описывающей 
динамические процессы, происходящие в пневматическом 
приводе при торможении — противодавлением, с 
ПОСТОЯННОЙ массой, заключенного В полостях 
пневмодвигателя воздуха, и с учетом изменения 
начальных параметров тормозного объема. 

Материалы и методы. Предложена математическая 
модель, описывающая изменение скорости движения 
выходного звена, давлений и температур в полостях 
пневмопривода в зависимости от начальных параметров 


рекуперативного объема. Решение математической 
модели осуществлялось методом численного 
интегрирования. 

Результаты — исследования. Получены зависимости 


скорости выходного звена, давлений и температур в 
полостях пневмопривода от начальных параметров 
рекуперативного объема. Адекватность полученной 
математической модели подтверждена по критерию 
Фишера. 

Обсуждения и заключения. Полученные результаты могут 
быть полезны для решения задач энергосбережения в 
пневмоприводах при организации торможения 
противодавлением. Применение рекуперативного объема 
повышает технологическую гибкость привода при его 


переналадке и расширяет возможности энергосбережения. 


Ключевые слова: пневматический привод, рекуперация, 
энергосбережение, противодавление, торможение, 
рекуперативный объем. 


* <“ 
Работа выполнена в рамках инициативной НИР. 
” Е-шай: апазио{епко @ уапдех.ги, раЙапа @ зраК-тай.га 
—” ТЬе гезеагсВ 15 допе уф шт 1е йаше оЁ фе шдерепдеги В &). 


тподисноп. Мефод$ о? епегоу зауте ш рпеитайс апуе аге 
сопз14егеа. ТЬе тефо4 о? БгаКте Бу сгеайпе БасК ргеззиге 
1 Фе ехБаи$( саущу оГ Фе рпеитайс асаагог 1$ оЁ ищегез. 
Опа4ег БгаКт>, Фе сотргеззе4 аш епегоу 15$ юге ш Фе 
тесирегануе уоште. [ 1$ роззе 10 сопго|] Фе БгаКшо 
Чупапс$ ШФгоиеб зеште Фе шша| рагатаег$ оЁ Фе 
тесирегайуе уоате. ТБе \отК оБесйуе 1$ № стеме а 
тафетайса| то4е| аезстте фе дупапис ргосез$ез (акт? 
расе ш Фе рпеитайс апуе ип4дег БгаКте Бу БасКргеззиге, 
УИ а соп%апе та$$ епс]озе4 ш Фе сам@ез оЁ Фе аш тогюг, 
ап соп$14етте уапайоп ог Фе ша! рагатеег$ оГ Фе 
БтаЮтс уоате. 

Маетпа!5 апа Мейо45. А тафетайса| то4е! 1$ ргорозеа 
ШФаЕ 4езстЬез Фе зрее4 свапое оГ Фе ошриЕ ПпК, ргеззигез 
ап {етрегаагез ш Ше сау@ез оГ Ше рпеитайс апуе 
дереп4т>® оп Фе шша| рагатеег$ оЁ Фе гесирегайуе 
уо[ате. ТБе зоааоп ю Фе тафетайса! тоде] 1$ саглей ош 
Бу Ше питепса! ицезогайоп тефо4. 

Кезеатсй Кезий5. 'ТБе дереп4депсез$ оЁ Фе ошриЕ Пик уе]оспу, 
ргеззигез ап4 {етрегаагез ш фе рпеитайс апуе саутае$ оп 
Фе ша] рагатеег$ оР Фе гесирегайуе уоате аге оатед. 
Адедиасу оГ Фе Би тафетайса| то4е!| 1$ сопйгте4 Бу 
ЕзБег$ сгцепоп. 

Р5сизяопту апа Сопсияоп5. ТБе геза $ оБаше сап Бе 
изе4 ю зоуе Фе ргоетз оР епегсу зауше ш рпеитайс 
Чтуез ипаег Фе огоап17айоп о? БасКргеззиге БгаК1п>. ТБе изе 
о гесирегайуе уоите шсгеазез Фе (есбпо]о21са| НехЬ|иу 
ог Фе апуе Чигтпе 15$ теадуазитепе ап ежеп4$ Фе 
роз Шве$ оЁ епегоу заупо. 


Кеужога: рпеитайс Апуе, тесиреганоп, епегеу зауше, 
БасКргеззиге, БгаК1п®, гесирегайуе уоате 


Машиностроение и машиноведение 


502 


Бир://уезбК.Аопза.га 


380 


Вестник Донского государственного технического университета. 2018. Т. 18, № 4. С. 379-384. 155М№ 1992-5980 е15 М 1992-6006 
Уезник о} оп кие Тесйтса! Отуегзиу. 2018. У. 18, по. 4, рр. 379-384. 155М№ 1992-5980 е155М 1992-6006 





Образеи для иитирования. Сиротенко, А. Н. Влияние Еог сйапоп: А.М. опоепко, 5.А. Райко, У/’ае|] заПоит. 
параметров рекуперативного объема на динамические ЕНесЕ оГ тесирегайуе уоте рагатаег$ оп Дупапис 
характеристики пневмопривода при торможении / свагасегтзИс$ оЁ рпеитайс апуе ипдег БгаК1п®. УезаиК оЁ 
А. Н. Сиротенко, С. А. Парко, Ваел Саллум // Вестник ОТО, 2018, уо. 18, по. 4, рр. 379-384. 
Дон. гос. техн. ун-та. — 2018. — Т. 18, № 4. — С. 379-384. Б@рз://401.ого/10.23947/1992-5980-2018-18-4-379-384. 


Ьрз://401.0г2/10.23947/1992-5980-2018-18-4-379-384. 


Введение. Пневматические приводы широко применяются для автоматизации и механизации 
вспомогательных технологических процессов. Энергопотребление пневмофицированного оборудования может 
составлять более 20% общего потребления предприятия. Поэтому вопросы энергосбережения в пневмоприводах 
являются весьма актуальными [1-3]. Такие особенности, как сжимаемость рабочей среды и инерционность 
выходных звеньев, затрудняют безударное торможение, усложняют управление и конструкцию привода [1, 3, 4]. 
Управлять законом торможения пневмопривода можно как воздействием на регулирующее устройство 
[1, 2, 5, 6], так и выбором рационального способа торможения [7-9]. 

Перечисленные особенности пневмопривода позволяют использовать сжатый воздух, как тормозной 
демпфер и накапливать энергию торможения, что эффективно реализуется при торможении пневмопривода 
способом противодавления [7]. Изменяя координату переключения на торможение, давление в тормозной и 
нагнетательной полости, подключая к выхлопной полости дополнительный объем, возможно не только 
обеспечить безударное торможение, но и рекуперировать энергию сжатого при торможении воздуха [2, 7, 10]. 
Задавая начальные параметры рекуперативного объема, можно влиять на тормозные и энергоскоростные 
параметры пневмопривода. Однако данный вопрос исследован в настоящее время недостаточно. 

Цель исследования — создание математической модели, описывающей динамические процессы, 
происходящие в пневматическом приводе при торможении противодавлением, с постоянной массой, 
заключенного в полостях пневмодвигателя воздуха и с учетом изменения начальных параметров тормозного 
объема. 

Постановка задачи. Требуется математически описать зависимость динамических характеристик 
пневматического привода при торможении противодавлением от начальных параметров рекуперативного 
объема. 

Математические зависимости. Рассмотрим торможение способом противодавления путем полного 
перекрытия магистралей нагнетания и выхлопа. Недостатком способа является то, что при изменении рабочего 
перемещения и внешней нагрузки на выходное звено необходимо заново настраивать координату включения 
торможения. Альтернативой является подключение в момент торможения к тормозной полости 
пневмодвигателя рекуперативного объема. Задавая его начальные параметры до торможения, возможно влиять 
на конечное давление в тормозном объеме и тормозной путь выходного звена пневматического привода [8]. 

Для математического описания газодинамических процессов, протекающих в полостях пневмопривода, 
были сделаны следующие допущения [10, 12]: давление в магистралях принимаем постоянным; 
термодинамические процессы, протекающие в полостях пневмопривода, рассматриваем как адиабатные; 
рабочую среду в полостях пневмопривода, сжатую при торможении, рассматриваем как идеальный газ. 

В основе торможения противодавлением лежит принцип создания силы сопротивления движению 
выходного звена пневмодвигателя. Это достигается частичным или полным перекрытием каналов нагнетания и 
выхлопа. Полное перекрытие более эффективно, так как отсутствует сброс сжатого при торможении воздуха из 
тормозной полости. Динамика пневмопривода до торможения определяется известной системой уравнений, 
учитывающей подкритический и надкритический режимы истечения [10]. 

Параметры пневмопривода показаны на рис. 1. Параметры рабочей среды обозначены как: У р ТГ, — 


объем, давление и температура сжатого в магистрали воздуха соответственно; У, р, Т» — объем, давление, 
температура воздуха, сжатого в поршневой полости пневмодвигателя; У„ р, Ти — объем, давление, 
температура воздуха сжатого штоковой или тормозной полости пневмодвигателя; У», Р»» Трек — Объем, 


давление, температура сжатого воздуха в рекуперативном объеме; Уж, рат, Т.т — объем, давление, температура 
атмосферного воздуха. Геометрические параметры пневмодвигателя имеют следующее обозначение: РЕ), Е, — 
поршневая и штоковая рабочие площади пневмодвигателя; }, /, — площади проходных сечений магистралей 
нагнетания и выхлопа соответственно; х, хо, Хо — координаты текущего перемещения, постоянного поршневого и 
штокового «пассивных» объемов соответственно; $ — максимальный рабочий ход; Р — полезная нагрузка. 


Сиротенко А. Н. и др. Влияние параметров рекуперативного объема на динамические характеристики пневмопривода при торможении 
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Рис. 1. Параметры пневмопривода при торможении противодавлением в рекуперативный объем 


При торможении в полостях пневмодвигателя находится постоянное количество сжатого воздуха. Считаем, что 
теплообмен с окружающей средой незначителен, поэтому термодинамический процесс считаем адиабатным, с 
показателем адиабаты К= [,4 [10]. Для торможения золотник распределителя переключится в нейтральное положение и 
перекроет проходные каналы магистралей нагнетания и выхлопа (},=0, р=0). Время переключения золотника 
распределителя не учитываем, каналы нагнетания и выхлопа перекрываются одновременно. Соединение рекуперативного 
объема с штоковой полостью пневмодвигателя считаем мгновенным. 


С учетом данных допущений полости, давление рабочей среды в поршневой полости пневмодвигателя 
представим в следующем виде: 


К 
Ра=(у У) ° Рипь (1.1) 
где Ри, Рит текущее и начальное Давление в момент переключения На торможение В поршневой ПОЛОСТИ 
пневмодвигателя; у, у» — текущий и начальный удельные объемы поршневой полости пневмодвигателя; К — 


показатель адиабаты. 
Сократим уравнение (1.1), разделив его на площадь поршня пневмодвигателя: 

Р"=((хол Нхи)/ (о хр», (1 .2) 
где хо», х„, Хх — координаты: начальная, включения торможения и текущего положения поршня пневмодвигателя 
соответственно. 

В момент переключения распределителя на торможении рекуперативный и тормозной объемы 
объединяются, в этом случае параметры сжатого воздуха определятся следующие системой (1.3): 


рой ри (1.3.1) 
13) шт шт 
| в (1.3.2) 
Р рек * рек > 


ГДе Ушт» Урек — Удельные объемы сжатого в начале торможения в штоковой полости и рекуперативном объеме 
воздуха; у/, у› — «условные» удельные объемы воздуха в штоковой полости и рекуперативном объеме; рить Рьек 
р — давления сжатого вначале торможения воздуха в штоковой полости, рекуперативном объеме и в 
«объединенном» объеме соответственно. 

В системе уравнений (1.3) выразим удельные объемы ъ,, у: 


К К 
у у 
Е ет А Е 
шт шт 1 шт | т — 
шт шт 
1 | 
к Е 
в м й (2) (1.4.1) 
т \ шт 1 
(1.4) р т 
у 
К К рек '. 
й —=рР'У_ —>Р -р —> 
рек рек й. Рек | т я 
рек рек 
1 1 
ь Е 
рек \° рек 2 (1.4.2) 
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где У„„» У»» — Штоковый и рекуперативный объемы соответственно на момент начала торможения,; У‚, У› — 


штоковый и рекуперативный объемы, соединенные в «объединенный» объем; т„и, Ш» — масса воздуха, 
сжатого в штоковом и рекуперативном объёме на момент начала торможения; р — давление в «объединенном» 
объеме. 
Для определения давления в «объединенном» объеме просуммируем уравнения (1.4.1) и (1.4.2): 
К 
1 1 
т К 
У | р +У РР 
шт \ о шт рек \` рек 
Г - У +у | 
го 


Используя уравнение Клапейрона, найдем температуру рабочей среды в «объединенном» объеме: 


"У и ек 
рУ=КТеТ=^- - м — х (1.6) 
_шт , РЕК 
у у 
шт рек 


где Т — температура сжатого при торможении в «объединенном» объеме воздуха; К — универсальная газовая 
постоянная... 

Окончательную зависимость температуры воздуха в тормозной полости пневмопривода от начальных 
параметров рекуперативного объема получим, раскрыв значения удельных объемов в уравнении (1.6) и 


сократив его на площадь штоковой ПОЛОСТИ пневмодвигателя: 


Г 
шт ек 
ТЕР: (5+ хо, +® РИН... 31 _->ЗНИИИИИНИНИЕ (1.7) 
рек (+0 -х).Т в +й Е 
Р шт рек рек 
где Ик — приведенный к штоковой площади пневмодвигателя, рекуперативный объем; Т, Тт Г Ре 


абсолютные температуры в «объединенном», штоковом и рекуперативном объемах соответственно. 











я 
а7х 
тории Ри" (1.8.1) 
Ти 
+х . 
Хх 
_| би т |. (1.8.2) 
| 
Ой ь 
фа шие Аб (рак 
5+ 0р -Х|(Ршт К +йрек | (Ррек ь 
Г (1.8.3) (1.8) 
: р- 
Е (5+%07 +Йрек —х 
— КЬ—1 
= р К 
= 
ре р в (1.8.4) 
— р 
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Движение поршня пневмодвигателя описывается уравнением (1.8.1). Общая система уравнений процесса 
торможения пневматического устройства будет иметь вид (1.8). Здесь: т — приведенная к штоку пневмодвигателя 
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масса рабочего органа; Т, — абсолютная температура воздуха в магистрали нагнетания; Т„ — начальное значение 
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температуры воздуха в поршневой полости пневмодвигателя в момент переключения на торможения. 

Уравнения (1.8.2) и (1.8.3) описывают изменение давления, а уравнения (1.8.4) и (1.8.5) — изменение 
температуры в поршневой и тормозной полостях пневмодвигателя соответственно. 

Совместное решение системы уравнений, описывающей динамику торможения пневмопривода (1.8) и 
его разгон методом численного интегрирования Рунге-Кутта [14], позволяет изучить зависимость динамических 
характеристик пневмопривода от начальных параметров рекуперативного объема. Стендовая оценка 
характеристического совпадения теоретических и практических зависимостей показала удовлетворительную 
сходимость [13], что позволило провести проверку адекватности полученной математической модели по 
критерию Фишера. Для этого был проведен двухфакторный эксперимент. Независимыми факторами были 
выбраны начальное давление и величина рекуперативного объема. Факторы изменялись согласно трем уровням 
варьирования, что позволило использовать девять комбинаций сочетания факторов при повторении каждого 
опыта трижды. 

Математическая модель адекватна практическим результатам, полученным при экспериментальной 
проверке, так как расчетное значение критерия Фишера составило 2,67, что меньше табличного (2,7) [15]. 

Полученные результаты позволили перейти к вычислительному эксперименту, который позволит 
определить рациональные комбинаций начальных параметров рекуперативного объема для заданного закона 
торможения пневмопривода и обеспечения при этом высоких энергоскоростных показателей. 


Выводы: 

1. Разработана математическая модель, адекватно описывающая динамические процессы, 
протекающие в полостях пневматического привода при торможении противодавлением с накоплением энергии 
сжатого воздуха в рекуперативный объем. 

2. Получены математические зависимости параметров пневматического привода от начальных 
параметров рекуперативного объема при торможении противодавлением. 
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Линейно-логический алгоритм принятия решения при обработке сигналов. 


жж 
В. С. Плаксиенко" 


"Институт радиотехнических систем и управления Южного федерального университета, г. Таганрог, Российская Федерация 
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у. $. Ракчепко"” 


Газе цие Рог Вадюесыиса! бузетл$ апа Сопи`о1, Зошфеги Еедега] Оиуег$йу, Тазапгое, Визз1ап Еедеганоп 


Введение. Для повышения эффективности приема и обра- 
ботки дискретных сигналов в условиях дефицита априор- 
ных сведений применяют эвристический синтез. Пред- 
ставлен анализ алгоритма принятия решения при линейно- 
логической обработке дискретных сигналов в случае, если 
априорные данные об их параметрах неполны. Цель рабо- 
ты — построение и анализ эффективности линейно- 
логических алгоритмов. 

Материалы и методы. Предложены эффективные в усло- 
виях априорной неопределенности математические алго- 
ритмы приема и обработки сигналов. Они основаны на 
учете структуры выбросов и превышений процессов в ка- 
налах обработки. 

Результаты исследования. Созданы линейно-логические 
алгоритмы обработки дискретных сигналов. Они основаны 
на учете одной, двух и более детальных характеристик 
выбросов или превышений случайных процессов. 
Обсуждение и заключения. Полученные результаты могут 
быть полезны в процессе синтеза алгоритмов и устройств 
приема и обработки сигналов. Алгоритмы и устройства 
реализуемы как в аналоговом виде, так и в виде алгорит- 
мов для ЭВМ. Созданы программы моделирования при 
обработке сигналов в условиях значительной априорной 
неопределенности информации о сигналах и о каналах их 
распространения. 


Ключевые слова: решающие правила, эвристический син- 
тез, прием и обработка дискретных сигналов, плотности 
распределения вероятностей, выбросы и превышения слу- 
чайных процессов. 
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Введение. Для оптимизации обработки сигналов в условиях значительной априорной неопределенности ис- 
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сигналов [1—7]. Теория выбросов широко используется в инженерной практике. Однако анализ закономерностей и де- 
тальных характеристик выбросов — сложная аналитическая задача даже при постоянном или медленно меняющемся 
пороге. При аналитическом подходе получаются соотношения, приводящие к несходящимся рядам, и этим объясняется 
отсутствие физически значимых результатов. 

В настоящее время и в ближайшем будущем перспективно применение дискретных многопозиционных 
сигналов. Это и дискретные адресные системы, и системы многократной телеграфии (частотная телеграфия, 
многопозиционная частотная телеграфия), и системы с Р и Е кодами, в которых последовательно-параллельно 
во времени передаются элементарные посылки длительностью Тх на одной или разных частотах [5, 8—10]. 

При приеме бинарных сигналов на входе решающего устройства присутствуют два случайных процес- 
са, и процедура принятия решения сводится к задаче статистической проверки гипотез. Необходимо опреде- 
лить, какой случайный процесс (из сравниваемых) обладает большей энергией на интервале наблюдения. Про- 
цедура принятия такого решения может быть сведена к сравнению величины разностного сигнала на выходе 
приемника с постоянным нулевым порогом. Таким образом, необходимо сравнивать минимум два случайных 
процесса (в лучшем случае анализируется случайный процесс с выхода приемника дискретных сигналов при 
медленно меняющемся пороге). 

Статистическая проверка гипотез сводится к анализу взаимных превышений двух и более процессов. 
Аналитическое представление данной задачи бывает громоздко и не дает инженерных решений [1, 11-15]. Эв- 
ристический синтез и статистическое моделирование на ЭВМ позволяют получить значимые инженерные при- 
ложения. 

При приеме многопозиционных сигналов можно упрощенно представить: 

— их спектры 5 (Г); 
— амплитудно-частотные характеристики (АЧХ) фильтров приемников К (}) (рис. 1). 
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Рис. 1. АЧХ фильтров приемников 


Для линейного приемника многопозиционных сигналов обязательно наличие: 

— общего усилителя промежуточной частоты (УПЧ), 
— общей схемы принятия решения (СПР). 

Необходимы также наборы разделительных фильтров (см. рис. 1) и амплитудных детекторов (АД). 

Анализ таких приемников показывает [1]: чем больше число позиций сигнала т, тем ниже их помехо- 
устойчивость. Требования к АЧХ фильтров (см. рис. 1) достаточно жесткие. АЧХ не должны перекрываться для 
обеспечения частотной ортогональности. При этом шумы на выходах фильтров будут независимы. Не должно 
быть областей перекрытия между АЧХ фильтров, настроенных на частоты }. (см. рис. 1). Форма АЧХ должна 
иметь плоский участок в окрестности частоты настройки, чтобы спектры сигналов не искажались. 

Следует отметить имеющиеся противоречия. Так, узкополосные фильтры ограничивают скорость рабо- 
ты. При замираниях сигналов и эффекте Доплера возможно ухудшение качества и даже срыв связи. Широкопо- 
лосные фильтры приводят к взаимопроникновению сигналов соседних по частоте каналов, т. е. нарушается ор- 
тогональность и, соответственно, падает качество приема. 

Для процессов на выходах фильтров соседних каналов в системах приема и обработки сигналов харак- 
терны статистические взаимосвязи, которые увеличиваются при увеличении взаимных перекрытий спектров 
сигналов либо АЧХ фильтров. Детальное рассмотрение процедуры разделения сигналов по частотному призна- 
ку показывает, что при обработке многопозиционных сигналов в случае перекрывающихся АЧХ фильтров 
можно говорить о попарно зависимых процессах в соседних по частоте каналах. Это имеет особое значение при 
использовании для оптимизации приема многопозиционных сигналов нелинейных либо линейно-логических 
процедур [1, 11-14]. 

Постановка задачи. В условиях значительной априорной неопределенности частота сигналов известна 
с ограниченной точностью, а диапазон варьирования длительностей элементарного символа может изменяться 


Плаксиенко В. С. Линейно-логический алгоритм принятия решения при обработке сигналов 


РакяепКо У. 5. [лпеаг-оэ1са] аестя1оп-тают2 @еотгийт ог я2па[ ргосе;ятя 





до тысячи раз. Дополнительное требование: работа в реальном масштабе времени. Ограничимся для простоты 
случаем обработки бинарных сигналов. 

Синтез алгоритма принятия решения. В работах [1, 11-13] рассматриваются особенности обработки 
дискретных сигналов в условиях значительной априорной неопределенности информации о длительности сим- 
вола. Задачу можно решить линейными процедурами путем параллелизма алгоритма принятия решения. Про- 
блема усложняется при повышении требуемой точности и расширении диапазона изменения длительностей 
элементарного символа. 

В таких условиях применение адаптивных процедур требует значительного времени на адаптацию или 
практически нереализуемо из-за высокой априорной неопределенности. Проанализируем возможности эври- 
стического синтеза алгоритма, основанного на учете статистических свойств взаимных превышений двух и бо- 
лее случайных процессов [1, 11-12]. 

Рассмотрим широкополосный прием с интегрированием или фильтрацией после детектора, когда 
АЕТ>> 1 (АР, — полоса пропускания приемника, Т — длительность элементарного символа). В этом случае от- 
счеты процесса на выходе решающего устройства можно считать распределенными по закону, близкому к нор- 
мальному. Математическое ожидание и дисперсия этого процесса определяются соотношениями [4]: 

АМ[х] = М[х] — М[х>], Ох] = Бр + В]. 

Вероятность ошибочного приема рассчитывают по формуле [1-3, 12]: 


Р= 5П- Фо) 


— отношение постоянной со- 


2 о. АМ й 
Здесь Ф(а.) = ,|-— | ехр(——)Ах — функция Крампа; @ = — — = 
п. 2 УВ 2+2АУТ/ 
2 


а 

<“ я 2 

ставляющей к эффективному значению переменной; И“ = р — отношение энергии элемента сигнала к спек- 
У 


тральной плотности помехи, где а — нормированная амплитуда сигнала, У” — спектральная плотность помехи. 

При обработке дискретных сигналов в указанных условиях возможно применение только непараметри- 
ческих алгоритмов принятия решения. Сами по себе интересны статистические характеристики выбросов слу- 
чайных процессов: длительность выбросов (превышений) #1, длительность интервалов между выбросами Ти, 
значения максимумов выбросов &„ ит. д. [3]. Алгоритмы оптимизации принятия решения могут быть основаны 
на учете одной из детальных характеристик, например: число пересечений порога за время наблюдения, дли- 
тельность интервалов превышения порога и т. д. Необходимо определить информативные признаки таких де- 
тальных характеристик взаимных превышений случайных процессов. Технически реализуемый инженерный 
учет этих данных позволит оптимизировать алгоритм приема и обработки сигналов в указанных условиях. 

С энергетической точки зрения наибольшую точность в условиях значительной априорной неопреде- 
ленности обеспечит учет площадей взаимных превышений $п (в пределе энергий). Однако это сложно реализо- 
вать в реальном масштабе времени алгоритмически и тем более технически. 

Понятие выбросов является частным случаем понятия взаимных превышений для двух и более процес- 
сов, в том числе случайных. При обработке бинарных сигналов ставится задача анализа выбросов разностного 
процесса относительно нулевого порога. Она адекватна задаче анализа взаимных превышений огибающих двух 
процессов в каналах обработки сигналов [3, 11]. 

Для оптимизации задачи принятия решения в реальном масштабе времени целесообразно учитывать 
одновременно две и более детальные характеристики выбросов — например, длительность и величину превы- 
шения. Произведение длительности на уровень (амплитуда) превышения является оценкой площади превыше- 
ния (с учетом коэффициента формы превышения) [11]. 

Воспользуемся двумя детальными характеристиками. 

1. Относительный уровень превышения 4: 

а = Х;(® / Х,(®, где Х, (1 и Х, (0) — анализируемые процессы. 

2. Относительная длительность превышений ®: 

О = ш / тк, где лк — интервал корреляции процессов на выходе линейной части приемника, определяемый его 
полосой пропускания. 

Совместный учет указанных характеристик технически сравнительно не труден. При этом его досто- 
верность близка к показателям, получаемым при учете площадей превышений 5п. Учитывая значение термина 
«превышение», 4 должно быть больше 1, поэтому в бинарной ситуации при формировании двухпараметриче- 
ских распределений 4 определяется соотношением [11]: 
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а х, (В) > х, (1); 
х. (Г) (1) 


Ч = 
Р- 
х, (1) 

Такой подход позволяет анализировать и многопараметрические распределения взаимных превышений 
двух и более случайных процессов. 

Для повышения качества принятия решения необходимо сократить общее количество взаимных пре- 
вышений процессов. Это следует из результатов анализа двумерных законов распределения взаимных превы- 
шений огибающей смеси сигнала и шума и огибающей шума [11]. Необходимо преобразовать процессы таким 
образом, чтобы, не нарушая отношения правдоподобия, получить формы двумерных распределений, легко раз- 


личаемые решающим устройством [12]. 
Обе задачи решаются применением модифицированного алгоритма комбинированного сложения [1, 11, 
14]: 
Хш (0 = [Х® - К Хх @а)] ИХ! - К Х› 1], 
(2) 
Хоп (® = [Х>® - К Хи] ИХ) - КХ, (1)], 


где К — коэффициент, принимающий значения от 0 до 1; 1[7(] — единичная ступенчатая функция, причем 
1172@)] = 1 при Г@ > 0и 117@)]| = 0 при (© < 0. 

Спектры процессов Хип (1) и Хоп (®, сформированных после обработки по алгоритму (2), расширяются. 
Их необходимо ограничивать до ширины спектров исходных информационных процессов Х1(1) и Х>(®. 

Значение К определяет реализуемые модификации алгоритма и устройств комбинированного сложе- 


ния. Речь идет об алгоритме взаимного преобразования, методе комбинированного сложения, системе пере- 
крестного запирания [1, 11]. 

При К = 1 из детальных характеристик превышений учитывается только факт превышений процессов. 
Реализуется алгоритм автовыбора при разнесенном приеме сигналов [4] и взаимное преобразование при приня- 
тии решения [11]. 

При К = 0,414 учитываются и факт превышения одного процесса другим, и уровень превышения а [4]. 
Реализуется комбинированное сложение в разнесенном приеме. При 1 > К > 0 учитывается не только факт пре- 
вышения, но и уровень превышения а. Реализуется модифицированный метод комбинированного сложения при 
принятии решения. 

Линейно-логические процедуры, реализуемые в результате эвристического синтеза, позволяют стати- 
стически преобразовать исходные случайные процессы. При К = 1 в (2) реализуется автовыбор процессов во 
времени. Непрерывные во времени процессы (рис. 2, а) превращаются в последовательности импульсов со слу- 


чайными амплитудами Х — их — (Г) (рис. 2, би2, в). 
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Рис. 2. Комбинированное сложение при К = 1 


На рис. 2, г представлены разностные процессы 7 = Х, - Х› без обработки и Й® =хи(? -— х›(? после об- 
работки по правилу (2). 


Плотности распределения вероятностей процессов ж > их — (1): 


У? < = Аб + ОЕ, (<) (3) 
И? © = 4,560 +, СОЕ (©). (4) 
Коэффициенты А и А, определяются из условия нормировки: 
А =1- | СОБ Сдах, (5) 
А, =1- | У, СЕ Сдах, (6) 


где Е(х) и ЕЁ (х) — интегральные функции распределения. 


Из рассмотрения рис. 2, а — 2, в следует, что часть времени процессы я” фи Хх _ (Г) равны нулю. 


Поэтому их плотности распределения вероятностей А (х) и АМ (х) будут содержать дельта-функции — (616) 


(рис. 3). 


му 


Ах] 





а) 6) 


Рис. 3. Плотности распределения вероятностей процессов 


На рис. 4 представлены плотности распределения вероятностей разностных процессов без обработки по 
правилу (2) \(<) (рис. 4, а) и после обработки \‘(7) (рис. 4, 6). Здесь А — зона стирания решающего 


устройства. 
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№ ма) 


Рис. 4. Плотности распределения вероятностей процессов (В, И“? (1) 


В соответствии с рис. 2, би 2, в процессы после обработки по правилу (2) — импульсные. Их спектры 
будут шире спектров исходных процессов. Таким образом, обработка по правилу (2) не приведет к повышению 
качества приема, а только улучшит пороговые свойства [11—12, 16]. 

Полная процедура обработки включает две операции: 

— линейно-логическую по правилу (2) (она нелинейна по сути, так как в спектре появляются составляющие, 
которых не было на входе); 
— фильтрацию. 

Фильтрация учитывает второй параметр превышений @® — относительную длительность. Расширивши- 
еся спектры процессов остаются в границах, задаваемых при формировании исходных процессов Х!() и 
Х>(1) (см. рис. 2, а). 

Обработка по правилу (2) в зависимости от диапазона изменения длительностей ожидаемых сигналов и 
условий приема [12] может повторяться. При этом значения параметра К различны, всегда меньше | и увеличи- 
ваются в последующих циклах обработки. 

Задачи автоподстройки частоты решаются более эффективно при использовании устройств, реализую- 
щих линейно-логическую обработку [14, 17—18]. Сравнительно несложными техническими средствами расши- 
ряются полосы захвата и удержания. Обеспечивается определенная устойчивость параметров автоподстройки в 
условиях воздействия аддитивных помех. 

Результаты исследования. Результаты исследования алгоритма (2) в случае, когда сигналы зависимы, 
выходят за рамки настоящей статьи. Однако уже полученные данные показали эффективность алгоритма 
вплоть до значений коэффициентов взаимной корреляции процессов р = 0,5-0,6. Это хорошо согласуется с ре- 
зультатами решения задач разнесенного приема [2—4]. Алгоритмы модифицированного комбинированного 
сложения сигналов эффективны и при разнесенном приеме, и в задачах принятия решения. Однако рассмотре- 
ние особенностей их реализации выходит за рамки статьи. 

Итак, оптимизация процедуры принятия решения при обработке дискретных сигналов в условиях зна- 
чительной априорной неопределенности может быть выполнена на основе непараметрических алгоритмов, ос- 
нованных на эвристическом учете детальных характеристик взаимных превышений случайных процессов. 
Большее количество учитываемых детальных характеристик усложняет алгоритм. Однако именно такой алго- 
ритм обеспечивает большую инвариантность по помехоустойчивости в условиях изменения длительности ожи- 
даемых сигналов. Синтезированные алгоритмы могут быть реализованы в виде процедур дополнительной об- 
работки на ЭВМ и технически — в виде аналоговых устройств [14, 19]. 
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Введение. В статье рассматриваются вопросы нормирования 
показателей надежности на ранних стадиях проектирования 
машин. Такой подход обеспечивает соблюдение заданного 
уровня надежности при разработке деталей, критичных с 
точки зрения стабильной эксплуатации. Цель работы — ис- 
следовать проблемы проектирования в указанной сфере. 
Анализ методов нормирования надежности позволяет утвер- 
ждать, что их недостаточно для проектирования. Кроме того, 
отмечены противоречия, связанные с использованием в про- 
ектировании экспоненциального закона надежности. 
Материалы и методы. На стадии технического задания 
определены численные значения наработки и безотказности 
машины. Затем значения вероятностей безотказной работы 
или отказов системы и ее элементов представлены степен- 
ными выражениями. В них степени являются параметрами 
ранжирования показателей безотказности элементов систе- 
МЫ. 

Результаты исследования. Синтезирован альтернативный 
подход к нормированию, который позволяет выполнить пол- 
ный структурный анализ проекта. Таким образом, могут быть 
оценены показатели надежности всей системы или деталей, 
безотказность которых определяет надежность машин. Дета- 
ли и другие элементы, отказы которых не внезапны, рассмат- 
риваются без применения экспоненциального закона. При 
этом сохраняется свойственная ему простота математических 
операций. 

Обсуждение и заключения. Для проектирования с заданным 
уровнем надежности недостаточно числового значения веро- 
ятности безотказной работы (ВБР) машины в целом, прини- 
маемого на стадии технического задания. Необходимы тре- 
буемые ВБР и наработки элементов системы, которые явля- 
ются источниками отказов, развивающихся по разным зако- 
нам. Полученные результаты могут быть использованы как 
при проектировании новых механических систем с заданным 
уровнем надежности, так и при модернизации машин. 


Ключевые слова: надежность, вероятность, безотказность, 
наработка, машина, техническая система, элемент. 
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Введение. Условия постановки на производство продукции машиностроения и единая система кон- 
структорской документации (ЕСКД) [1-2] требуют, чтобы при разработке проекта машины учитывались сле- 
дующие параметры надежности технических систем (ТС): 

— долговечность — календарный срок службы Ту; 
— безотказность — вероятность безотказной работы Р(|...); 
— наработка { = Т часов (на этот период должна быть гарантирована ВБР). 

Данные показатели согласуются на стадии технического задания (ТЗ) и характеризуют машину в це- 
лом. Разработчик должен повлиять на безотказность, которая является главным свойством надежности невос- 
станавливаемых ТС и в межремонтные периоды работы формирует долговечность восстанавливаемых ТС. 

Наработка { является аргументом функции надежности Р(®. По любому закону надежности Р(Ф моно- 
тонно убывает с наращиванием 1. В силу этого для достижения требуемой ВБР наработка должна быть рацио- 
нально минимизирована и выбрана в соответствии с логикой работы ТС [3]. 

Детали и невосстанавливаемые изделия, вырабатывая свой ресурс, определяют наработку системы на 
отказ. В идеале наработки элементов на отказ должны совпадать или хотя бы быть кратными выбранной нара- 
ботке ТС. При этом календарный срок службы Ту и требуемая наработка Т связываются соотношением: 

Т=Т, -365-К „-24-К- ПВ, ч, (1) 


где Т, — календарный срок службы в годах; К». — коэффициент возможного использования ТС в год; Кеш — 
коэффициент использования ТС в сутки; /1В — относительная продолжительность включения, которая пред- 
ставляет собой осредненное отношение времени включения (время разгона и время установившегося движе- 
ния) ко времени рабочего цикла ТС. 

Высокая вероятность безотказной работы (ВБР) возможна при снижении требуемой наработки. Если ее 
необходимо существенно увеличить, безотказность придется поддерживать следующим образом: 
— чаще выполнять техническое обслуживание (ТО); 
— при проектировании ключевых элементов закладывать более высокие затраты на обеспечение качества изде- 
лий [4,5]. 

ВБР проекта может быть выбрана из существующих отраслевых нормативов, из конкурентных условий 
и по другим основаниям, в том числе из условных классов надежности, принятых для изделий машинострое- 
ния (табл. 1). 

Таблица 1 
Классы надежности изделий машиностроения 


ОТ ОО ОО ПОЕМ ПО ОЕ 





Допустимое значение ВБР Р(1) < > 0,9 > 0,999 | >0,9999 


При рассмотрении классов надежности учитываются перечисленные ниже группы характери- 
стик ТС [3, 6, 7]. 

1. По структурному типу изделий. Рассматриваются технологические комплексы, машины, агрегаты, 
механизмы и иерархия (уровни) их сборочных единиц и деталей. При последовательном взаимодействии эле- 
ментов в большинстве машин ВБР растет с переходом к нижним уровням структурной схемы надежности. То- 
гда ВБР деталей при структурном анализе может приближаться к показателям 3—5 классов надежности. 

2. По видам ТС. Планируемый уровень ВБР напрямую зависит от уровня ответственности производи- 
теля перед потребителем ТС. Безусловно, высочайший уровень ВБР планируется, например, для летательных 
аппаратов, машин химических производств, медоборудования, военной техники ит. п. 

3. По видам последствий отказов. Планируемый уровень ВБР напрямую зависит от потенциального 
уровня ущерба при отказе машины (экономического, экологического, репутационного и пр.). 

Общепринято [5, 6, 8—13] представление отказов ТС как внезапных. В этом случае ВБР описывается 
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Здесь интенсивность отказов Л. в период нормального функционирования после приработки связывается со 
средней наработкой на отказ Т системы 


Л. =1/Т (3) 
с параметрами безотказности 
Л. =(-ш РСГУ/Т. (4) 
Отсюда может быть рассчитана средняя наработка на отказ ТС 
Т=-Т/( т РСГ)). (5) 


Этот алгоритм формирования параметров надежности ТС имеет как сторонников, так и крити- 
ков [7, 10, 14]. Все ли отказы происходят внезапно? Насколько верны параметры надежности ТС (3)—(5), выте- 
кающие из экспоненциального закона? Эти вопросы возникают при определении интенсивности отказов (4) и 
соответствующей ВБР ТС (3) вне зависимости от методики оценки: 

— по заданным наработке Ти ВБР Р(Т); 
— по средней наработке на отказ 7; 
— по результатам испытаний 7. 

В литературе структурный анализ надежности ТС, преобразующих входное воздействие (ВВ) в выход- 
ные параметры (ВП) принято представлять как восходящий схемный анализ от ВБР элементов к ВБР систем. 
При этом рассматриваются структурные схемы: последовательные (рис. 1, а), параллельные (рис. 1, 6) и комби- 


нированные (рис. 1, в). 








вв>|пс1 |>пс2 |>...>|ПС:> > псп |>ВП 


в) 


Рис. 1. Структурные схемы надежности технических систем при различных видах взаимодействия элементов 


При последовательном взаимодействии элементов ТС (см. рис. 1, а), отказы которых независимы, экс- 
поненциальный закон (2) представляет удобную математическую модель. Если ВБР системы Р(® и элементов 


— подсистем Р (1) находятся в соотношении 


п 
РО =| РО, (6) 
1=1 
то интенсивности отказов системы А, И подсистем А; находятся в соотношении 
п 
л=>^, (7) 
1=1 


и ранжируются в зависимости от принятого принципа распределения интенсивностей по элементам. 

Иерархическая структурная схема, сформированная из вложенных уровней подсистем и элементов [3, 6, 8-12], 
позволяет выборочно или по всей структуре выполнять нисходящий анализ [3, 9] (от ВБР ТС к ВБР элементов). Требу- 
емые значения ВБР элементов вместе с требуемой наработкой являются исходными данными для проектирования де- 
талей с заданным уровнем надежности [15, 16]. 

Алгоритм нисходящего анализа на базе (2) позволяет одинаково успешно распределять интенсивности и ВБР 
элементов по всем вложенным уровням структурных схем. Это справедливо для радиоэлектронных систем [10, 11]. Но 
на уровне деталей механической системы [3] возникает противоречие: нельзя применять экспоненциальный закон, если 
развитие постепенного отказа протекает по другому закону. Отказ системы произойдет в результате последовательного 
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взаимодействия, не связанного с экспоненциальным законом. В то же время разработка структурной схемы 
надежности ТС имеет смысл, если она служит основанием для проектирования деталей и подбора стандартизо- 
ванных элементов с требуемой ВБР на выбранную наработку. 

При анализе ВБР ТС с параллельным взаимодействием элементов (см. рис. 1, 6) учитываются следую- 
щие условия: 
— элементы включены постоянно, 
— их отказы независимы, 
— каждый из т элементов имеет ВБР Р,@), 


— каждый из т элементов способен принять входное воздействие и преобразовать его в выходной параметр 
ТС (см. рис. 1,6). 

В этом случае отказ ТС произойдет после отказа последнего работоспособного элемента. Вероятность 
отказа системы через вероятности отказов элементов: 


ЕЮ =| Е. (8) 
7=1 
Из свойства 
РО-ЕСа)=1 (9) 
вероятность безотказной работы системы через ВБР элементов: 
Р(® =1-] [П-Р,/(@]. (10) 
1=1 
Для однородных элементов 
Ра) =1-[1-Р(0 |. (11) 


Если цель — обеспечение заданной ВБР ТС и учитываются однородные резервированные элементы, то 
использование выражения экспоненциального закона (2) дает простой расчет интенсивности отказов элементов 
или числа элементов т только в соотношении (11). Для более сложных соотношений (10) и комбинированных 
структур (см. рис. 1, в) примеры использования экспоненциального закона не приводятся. 

Простое разделение на элементы невозможно, если речь идет о сложных структурах с разнонаправлен- 


ным взаимодействием элементов, имеющих зависимые отказы. Параметры функционирования и надежности 
должны определяться для системы в целом, а для этого, возможно, потребуется выполнить большой объем ана- 
литических и экспериментальных работ. Графическое представление таких схем [7, 9] сопровождается упоми- 
нанием о сложности модели функционирования и громоздкости расчетов надежности (например, сложное за- 
мкнутое взаимозависимое функционирование блоков двигателя внутреннего сгорания или реактивного двига- 
теля). 

Основная часть. Для конструктивной реализации ТС с любым видом взаимодействия необходимо из- 
готовить детали, из которых составляют сборочные единицы, а затем функциональные модули (механизмы, агрега- 
ты, блоки питания, блоки управления и т. п.). Современный модульный принцип построения машин предусматрива- 
ет модульные операции: разработку, сборку и отладку, модификацию и модернизацию, ремонт, замену и т. п. Детали 
и прочие элементы находятся на нижнем уровне иерархической конструктивной системы, состоящей из вложенных 
блоков. Используя ее в качестве структурной схемы надежности, необходимо найти критерий распределения ВБР по 
элементам. Следует рассмотреть не функцию Р(®, а ее численное значение Р(Т). В этом случае можно подставить 
заданную наработку Ти при распределении численных значений ВБР по элементам использовать в качестве крите- 
рия стоимости восстановления структурных единиц в случае их отказов. В стоимость восстановления могут быть 
включены: 

— стоимости используемых для ремонта материалов, изделий, диагностики и ремонтных работ; 
— оцененные в денежных единицах последствия отказа (восстановление других отказавших элементов при зависи- 
мых отказах, потери от простоя ТС, страховые выплаты на возмещение ущерба и др.). 

Задачей настоящего исследования является разработка метода нисходящего анализа структурных схем 
надежности, устраняющего отмеченные выше противоречия. Структурная схема надежности ТС должна быть ин- 
струментом распределения ВБР элементов по заданной ВБР ТС на начальных стадиях проекта. Тогда выбранная 
наработка и сбалансированные ВБР элементов будут исходными данными для проектирования с заданным уровнем 
надежности. Для этого структурная схема должна удовлетворять ряду условий. 

1. Структурная схема должна основываться на структуре конструкции. Это упростит сложные функциональ- 
ные связи между элементами 

2. Структурная схема должна содержать математическую модель надежности, доступную для расчетов на 
ранних стадиях проекта. 
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3. Модель надежности по структурной схеме должна опираться не на экспоненциальный закон, а только 
на фундаментальные свойства надежности технических объектов и систем (6), (8)—(11). 

4. Критерием ранжирования ВБР элементов должна быть стоимость восстановления элемента в случае 
его отказа. Ее можно рассматривать как денежный эквивалент возмещения ущерба от отказа. 

5. Нисходящий структурный анализ должен применяться одинаково эффективно как для последователь- 
но взаимодействующих, так и для резервированных и комбинированных схем надежности. 

Альтернативное представление численного значения ВБР объекта. Пока время { имеет значение 
аргумента функции безотказности Р(1), она представляет собой функцию, убывающую по любому из известных 
законов или полученную статистическим путем. После выбора требуемой наработки 1 = Т ВБР получает значе- 
ние Р(г) = Р(Т) в границах 0 < Р(Т) < 1. Такое число может быть представлено множеством способов, из кото- 
рых выбирается показательное выражение 


Р(Г) = В^. (12) 
Х определяется из соотношения 
Х = 18 (Р(Г)) / [в В. (13) 
Приняв значение основания степени В = 10, получим выражение численного значения ВБР: 
РГ) =10^, (14) 
откуда 
Х = №® РТ). (15) 


Показатель степени Х назовем параметром ранжирования (ПР) ВБР ТС. В приложении к элементам си- 
стемы показатели степени Х, также именуются параметрами ранжирования ВБР элементов. Далее в этом 


представлении рассматривается «самый простой и самый важный случай» [13] надежности систем. 

ВБР элементов ТС при их последовательном взаимодействии. Рассмотрим структурную схему 
надежности ТС с последовательным взаимодействием элементов (см. рис. 1, а), отказы которых независимы. В 
этом случае отказ системы при отказе любого из и элементов в момент времени { = Т выражается через ВБР 
элементов согласно (6): 


Рт)= Ра). 


В (6) произведение может быть получено множеством комбинаций Р(Г) и видов их представления. 
Применив (14) для численных значений ВБР системы Р(Г) и ее элементов, произведение 
РГ) =В (Г)хР(Г)х...х Р(Т) можно представить: 


10° =10^' х10^? х...х10^", 
откуда связь ПР системы и элементов: 
ово (16) 


Из множества возможных комбинаций Х, в (16) выбираются вытекающие из стоимостей восстановле- 


ния отказавших элементов # : При этом элементы с большей стоимостью восстановления ДОЛЖНЫ обладать 


большими ВБР Р(1). То есть при стоимостях восстановления элементов Ст, С, ..., Спи набору ПР элементов Х, 
,Х,,.... Х, должен быть сопоставлен набор обратных величин стоимостей 1/С!, 1/С,, ..., 1/С.. Указанное 
соотношение может быть записано суммами 
1 1 1 
1= я + в О 
о 
С С С. 


в которых принятые условия для ПР Х, получаются из почленного равенства слагаемых: 


1 


в 


Хх 
С, 
Правую часть этого равенства Назовем «весовой коэффициент стоимости восстановления последова- 
тельно взаимодействующих элементов»: 








1 


а. = ы (177) 


Ой 1 
2. 


1 
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Значения ПР Х, ‚ отвечающие условию ранжирования: 
Х, = Ха,. (18) 
Единица измерения стоимости восстановления не имеет значения, так как в (17) используются отноше- 
ния стоимостей. Нисходящий анализ ВБР элементов ТС с последовательным взаимодействием рассмотрен в 
примере 1. Представленная ниже система содержит для простоты три элемента. Однако для действия алгорит- 
ма, основанного на (14—18), возможно любое число элементов. 
Пример 1. Расчет вероятностей безотказной работы элементов ТС с последовательным взаимодействи- 
ем (см. рис. 1, а). Исходные данные: 
— ВБР системы Р(Т) = 0,9; 
— число элементов и = 3; 
— стоимости восстановления элементов (в у. е.): С! = 5000, С, = 3000, С, = 2000. 
ПР для ВБРТС (15): Х = е РСТ) = 1 0,9 =-0,04576. 
Результаты пооперационного расчета ВБР элементов представлены в табл. 2. 
Таблица 2 


Расчет ВБР элементов ТС с последовательным взаимодействием 


| 06 | був | Ма | & 0 | Х9| ва — 
С О О ОИ ОС О 
лет [500 [00% [0103588 200086 [о9жи  — 
Элмет?_ [300 [00003 [одоэы 00476 0564 
лем 


Элемент 3 2000 0,0005 0,483871 | —0,02214 | 0,950297 


Королыыесуиыя | [О [Ха [Ух | "РП 


А 1000033 [1 | 10.0456 109 


ВБР элементов ТС при их параллельном взаимодействии. Рассмотрим следующую ситуацию: отказ 





системы при отказе всех элементов с вероятностью А(Т) в момент времени {1 = Т, соответствующий выбранной 
наработке. В этом случае основное свойство параллельно взаимодействующих элементов (см. рис. 1, 6) соглас- 
но (8) выражается через вероятность отказов элементов Р,(Т): 


Ег)= ПЕ). 


= 
Принимая по аналогии с (14) 


ЕСГ) =10”, (19) 
где У определяется из логарифмического уравнения [е (Е(Т)) = У 1210: 
У=№ю ЕТ). (20) 


Показатель степени У назовем параметром ранжирования вероятности отказа ТС. Если речь идет об 
элементах системы, показатели степени У, называются параметрами ранжирования вероятности отказа элемен- 


тов. Условия (8, 19) могут быть представлены произведением 10^ =10'' х10’” х...х10" , откуда 
У ЕР (21) 
Из множества возможных комбинаций У, в (21) выбираются вытекающие из стоимостей восстановле- 
ния отказавших элементов С,. При этом элементы с большей стоимостью восстановления должны обладать 
меньшими вероятностями отказа (что также означает большие ВБР). То есть при стоимостях восстановления 


элементов С, С,, ..., С, ряду показателей степени У,У,,....У, должен быть сопоставлен ряд стоимостей 


1? 2 5***9 т 


С,С.,....С„ . Указанное соотношение может быть записано суммами 











У у, 7 С С, С 
1+ —+—+...+—“ и 1= + + м. 
о 


Здесь выбранные условия для показателей степени при У, получаются из почленного равенства слагае- 


) С, 


МЫХ: у —= уха 


раллельном взаимодействии элементов»: 





. Правую часть равенства назовем «весовой коэффициент стоимости восстановления при па- 
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С. 
р (22) 


, 
2.С, 
Значения ПР У, отвечающие сформулированному выше соотношению: 
7, =УЬ,. (23) 


Нисходящий анализ вероятностей отказов и ВБР элементов ТС с параллельным взаимодействием рас- 





смотрены в примере 2. Система состоит из трех элементов. 
Пример 2. Расчет вероятностей отказа и ВБР элементов технической системы с параллельным взаимо- 
действием (см. рис. 1, 6). Исходные данные: 
— ВБР системы Р(Т) = 0,9; 
— число элементов 1 = 3; 
— стоимости восстановления элементов (в у. е.): С, = 5000, С, = 3000 и С, = 2000. 


Рассчитаем показатели надежности системы. Вероятность отказа системы из (9): 
ЕСГ) =1-Р(Г) =1- 0,9 =0,1. 
Значения ПР для вероятности отказа системы из (20): 
У = е ЕТ) = 180,1 =-1. 
Результаты пооперационного расчета приводятся в табл. 3. 
Таблица 3 


Расчет ВБР элементов ТС с параллельным взаимодействием 


о 


5000 031628 05583775 
3000 05501187 0498813 


2000 0.530957 0369043 


Контрольные суммы ЕСГ) = И ЕСГ) | РОГ) = 1- [а —Р.Г)) 
и произведения 
С ОО ООО ООО ООО 
Пример 3. Расчет вероятностей безотказной работы элементов комбинированной ТС (см. рис. 1, в). Ис- 
ходные данные: 
— число последовательно взаимодействующих подсистем и = 3; 
— число параллельно взаимодействующих элементов подсистемы ПСЗ т = 4; 
— ВБР системы Р(Т) = 0,95; 
— стоимости восстановления подсистем (в у. е.): С; = 5000, С, = 8000, С, = 12000; 





— стоимости восстановления параллельно взаимодействующих элементов (ву. е.): С,, = 3000, С,, = 4000, 


С,; =5000, С, . = 6000. 


Рассчитаем параметры надежности системы. ПР для ВБР ТС из (15): 
Х = ю РСТ) = 41° 0,95 =-0,022276. 
Результаты пооперационного расчета приводятся в табл. 4 и 5. 
Таблица 4 


Расчет ВБР подсистем ТС с последовательным взаимодействием 


обет | бе | Иа [ аб | ха | ва | | 
О ПО ПО О МС ОСИ О 
[о | я [обожю [одвузе [| вожюи [| обв | 
п | вю [обоз [ов | бювю | он | 
| 

Г 
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Из табл. 4 ВБР подсистемы ПСЗ Р.(Т) = 0,989587. Вероятность отказа ПСЗ (8) 
ЕСГ) =1-В(Г)=1-0,989587 = 0,010413 . ПР для вероятности отказа ПСЗ (22): 
У = (Т) = 180,010413 = —1,982407 . 


Результаты пооперационного расчета приводятся в табл. 5. 
Таблица 5 


Расчет ВБР элементов ТСЗ с параллельным взаимодействием 


| | 


3 


5000 


—_ [жми | об | 0959 
05532697 
4000 0.637369 
0718593 


Контрольные суммы и уз С, >. У, Е(Г)= П ЕАГ) | Рб)= 1- [а —Р‚(Г)) 


произведения 
18000 —1,982407 0,010413 0,989587 


Следует принять требуемую наработку ТС и всех ее элементов (1), чтобы ВБР элементов, рассчитан- 





ных в табл. 2—5, были исходными данными для расчета и подбора элементов с заданным уровнем надежности. 

Необходимо определить полученные ресурсы по соответствующим моделям надежности с расчетными 
ВБР (табл. 2, 3, 5). Для этого учитываются размеры и конструктивные параметры элементов (деталей, стан- 
дартных и других невосстанавливаемых изделий), полученные на стадии эскизного проекта. Затем могут быть 
приняты меры для сближения полученных и заданных ресурсов [15, 16]. 

Заключение. Разработан и проверен на численных примерах метод нисходящего анализа структурных 
схем надежности машин. Его алгоритм совпадает с последовательностью стадий проектирования: от ВБР ма- 
шины — к ВБР ее деталей. Выбор требуемой наработки и распределение ВБР по деталям на начальных стадиях 
проекта обеспечивают возможность проектирования с заданным уровнем безотказности. Структурная схема 
может базироваться на структуре конструкции. Таким образом, обеспечивается универсализация подхода к 
формированию структурных схем надежности при различных видах взаимодействия элементов. Критерий ран- 
жирования ВБР элементов — стоимость восстановления элемента в случае его отказа. В нее на любом уровне 
структурной схемы включается денежный эквивалент стоимостей материалов, изделий, ремонтных работ, 
ущерба от отказа. Они с достаточной точностью известны на начальной стадии разработки машин. 

Устранена необходимость использования экспоненциального закона при нормировании надежности и 
разработке ее структурных схем. Анализ ВБР опирается только на фундаментальные свойства надежности, об- 
щие для всех технических объектов и систем. Алгоритм анализа сохраняется на всех уровнях структурных 
схем, включая уровни деталей. 

Математические модели структурных схем просты, доступны для расчетов на ранних стадиях проекта, 
одинаково эффективны для последовательно взаимодействующих, резервированных и комбинированных схем 
надежности. Они дают удобные возможности алгоритмизации, программирования и коррекции данных. 
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Введение. Исследовано действие варьируемых параметров 
процесса электроакустического напыления (ЭЛАН) на 
характеристики кристаллической структуры упрочняю- 
щих покрытий. Выявлены значения указанных парамет- 
ров, обеспечивающие получение наноструктурных защит- 
ных покрытий деталей машин и режущего инструмента. 
Упрочнение при помощи подобных покрытий позволяет 
достичь значительного (в 5—10 раз) повышения ресурса 
работы упрочняемых деталей и инструментов, предназна- 
ченных для механообработки. Для получения покрытий с 
заданными свойствами поверхностного слоя следует вы- 
бирать нанокристаллические материалы. При этом допу- 
стимо определенное содержание аморфной фазы. 
Материалы и методы. Использован рентгеноструктурный 
анализ, который проводился на рентгеновском дифракци- 
онном аппарате отечественного производства «ДРОН- 
ЗМ». Метод Шеррера — Вилсона применен с целью опре- 
деления зернистости блоков частиц по значению физиче- 
ского уширения пиков дифрактограммы. Выводы по ре- 
зультатам работы основаны на методике разделения вкла- 
дов значащих факторов в уширение пиков дифракцион- 
ных отражений (метод Уоррена — Авербаха). 
Результаты исследования. В зависимости от технологи- 
ческих режимов и метода получения наноструктурных 
материалов происходит неоднозначно интерпретируемое 
изменение показателей уширения дифракционных пиков, 
что, в общем, характерно для нанокристаллических ме- 
таллов. Это явление можно объяснить, в частности, нали- 
чием наноразмерного эффекта в упрочненном слое. При- 
сутствие нанокристаллической структуры в напыленном 
слое напрямую подтверждается вычисленными значения- 
ми размеров областей когерентного рассеивания (ОКР). 
Наличие значимых величин угла разориентирования 
структуры кристаллов подтверждается величиной ОКР 
для исследованных рефлексов 110 и 220, что подкрепляет- 
ся высоким процентным содержанием аморфной фазы. 
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Обсуждение и заключения. Метод электроакустического 
напыления перспективен для получения нанокристалличе- 
ских структур в поверхностном и приповерхностном сло- 
ях напыленных образцов. Изменение технологических 
параметров ЭЛАН приводит к разбросу значений пара- 
метров кристаллической решетки и областей когерентного 
рассеивания. В этом случае не выявляется определенная 
тенденция. В дальнейшем предполагается решение данной 
проблемы. Во-первых, будут проведены эксперименты с 
целю определения оптимальных режимов напыления, 
способствующих образованию — нанокристаллических 
структур. Во-вторых, планируется визуальное наблюдение 
и оценка структуры напыляемого слоя при помощи элек- 


тронной микроскопии. 


Ключевые слова: нанокристаллические 


структуры, 
упрочняющие покрытия, кристаллическая решетка, ди- 
фрактометрические исследования, электроакустическое 


напыление. 
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Введение. Получение наноструктурных упрочняющих покрытий является актуальной и перспективной 
технической задачей [1, 2]. Многочисленные исследования в области материаловедения показали, что ощути- 
мое (в несколько раз) изменение прочности, твердости, износостойкости материалов возможно при уменьше- 
нии кристаллических зерен до определенного значения. 

В рамках данного исследования предполагается изучение кристаллической микроструктуры упрочня- 
ющих покрытий с использованием способа электроакустического напыления (ЭЛАН) [3, 4, 5]. Такое упрочне- 
ние предполагает использование электрической искры (высококонцентрированный поток энергии) и ультразву- 
ковых колебаний продольно-крутильной направленности. ЭЛАН-упрочнение с использованием специализиро- 
ванной установки позволяет в 5—10 раз повысить ресурс работы деталей и инструментов. 

В работе технологические параметры ЭЛАН варьируются с целью выявления их значений, способ- 
ствующих получению наноструктурных защитных покрытий деталей и режущего инструмента. Для создания 
защитных покрытий с требуемыми свойствами нередко используют сплавы в аморфном состоянии. Однако по 
распространенности нанокристаллические материалы в некоторых областях уже опередили аморфные. Нано- 
кристаллы в гораздо меньшей степени релаксируют в процессе использования покрытия и по свойствам не 
уступают аморфным материалам. Очевидно, что для получения покрытий с заданными и стабильными свой- 
ствами поверхностного слоя следует выбирать нанокристаллические материалы, пусть и с определенным со- 
держанием аморфной фазы. 

Материалы и методы. Для выявления в поверхностном слое нанокристаллических материалов актив- 
но используются резонансные методы, методы рентгеноструктурного анализа, электронная микроскопия и 
др. [2]. В представленной работе для оценки параметров микроструктуры применяется наиболее доступный 
метод, основанный на выявлении уширений дифракционных отражений при рентгеноструктурном анализе об- 
разцов. Для рентгеноструктурного анализа использован рентгеновский дифракционный аппарат российского 
производства «ДРОН-3М». Анализировалось характеристическое излучение железного анода с выделением 
ЛК. — спектральной линии образца, а также детектора. Характеристики рентгеновской трубки: 25 кВ, 5 мА. 
Использовался метод фокусировки излучения по Брэггу — Брентано [6]. Выбранные параметры детектора: ско- 
рость движения — 1 град/мин, время интегрирования — 5 с. При определении расстояний между плоскостями 
измерения ошибка составила + 0,001 А, что аналогично определению положения дифракционных максимумов с 
точностью + 0,02 градуса. Исследовались образцы цилиндрической формы из стали 45. Диаметр и высота об- 
разца — 6 мм. 
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Для идентификации образцов на лицевой поверхности каждого из них наносились специальные метки 


при помощи ЭЛАН. В качестве эталона был выбран не близкий к чистому @-Ее, а образец без напыления из ста- 


ли 45. 


Результаты исследования. Полученные в ходе работы профили представлены на рис. 1. 
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Рис. 1. Дифракционные отражения с профилями: эталон с крупным размером зерен (а); 
образец для напыления (И = 16,9 В, А = 4,9 мкм) (6); И = 13,1 В, А = 15,2 мкм (в); И = 12,9 В, А =9,9 мкм (г); 
(И = 12,9 В, А=5,1 мкм (д); 0 = 9,1 В, А = 15,1 мкм (е) 


Форма дифракционных профилей для эталонного образца приведена на рис. 1 (а), для образцов с нане- 


сенным покрытием на приведенных режимах — на рис. 1 (6) — (е). Варьируемые параметры ЭЛАН — напря- 


жение электрода-анода (7) и величина амплитуды (А) колебаний ультразвуковой частоты. 


Дифрактограммы эталонного и исследуемого образца демонстрируют отношение значений интенсив- 


ностей. По этим данным рассчитывалось соотношение кристаллизовавшегося и аморфного материала. Вклад 


текстурной составляющей не учитывался. 


Анализ экспериментальных данных позволяет сделать вывод, что ОСНОВНОЙ фазой в большинстве ис- 


следуемых образцов является а — Ее. Метод Шеррера — Вилсона (наиболее подходящий в данном случае) был 


использован с целью определения зернистости блоков частиц по значению физического уширения пиков дДи- 


фрактограммы [7]. В табл. 1 представлены сводные данные по результатам исследования: 


— физическое уширение В, 


— _ межплоскостное расстояние Ч, 


— размерный параметр кристаллической решетки а, 


— размер областей когерентного рассеивания. 
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Таблица 1 


Расчетные параметры кристаллической структуры и данные дифрактограмм 


Параметр / об- ны 550 


— яз т ПОЕМ ЗИ МАЕ ЗИ 


рока [тя [алло | [а [я лы 


20, град 26,916 | 56,86 | 56, 926 56,876. [56.876 876 | 56,976 Е: 144,66 44 44 [144.78 [144.78 144,91 1144,53 | 144,4 


2,026 | 2,028 | 2,025 | 2,027 2,019 1,0129 | 1,0126 | 1,0122 | 1,0132 | 1,0129 


| «А [28553 2,8581 | 2,8539 | 2,8567 | 2,8525 Ета 2,8553 | 2,8793 | 2,8571 1 2,8571 
ан. 


Содержание 
аморфной 59 55 27, 53,5 51 
компоненты, % 


В качестве основных выделены рефлексы 110 и 220. 

Данные табл. 1 ясно говорят о значительном уширении дифракционных пиков. Логично предположить, 
что такая ширина в нашем случае обусловлена в большей мере мелкодисперсной структурой кристаллитов, чем 
дислокационными микронапряжениями поверхностного слоя. Для подтверждения данного вывода вклады от 
указанных факторов оценены методом Уоррена — Авербаха [8]. В пользу сформулированного выше предполо- 
жения косвенно свидетельствует смещение на дифрактограммах центра тяжести пиков с одновременным уве- 
личением параметра решетки а до 1,15-—1,17 %. Само смещение невелико, но ощутимо (выражается в сотых 





долях градуса). 
Рассмотрим явления, происходящие при изменении значений технологических режимов (электрическо- 
го напряжения и амплитуды колебаний), и их влияние на величину уширения пиков дифрактограммы в зависи- 


мости от размерного параметра решетки (рис. 2). 
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Рис. 2. Влияние варьирования режимов ЭЛАН на зависимость параметра решетки а от величины 
уширения пиков дифрактограмм В 
Как и предполагалось, изменение режимов напыления влияет на параметр решетки а без определенной 
прослеживаемой зависимости. При этом изменения значений а прямо пропорциональны величине уширения 
линий дифрактограммы. В работе [9] сравнимые пропорциональные зависимости можно наблюдать для раз- 
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личных оксидов. 
Такие неоднозначно интерполируемые зависимости свойственны наноструктурным материалам (в 


частности, металлам) и определяются физическим способом их синтеза [1, 10] (рис. 3). 
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Рис. 3. Влияние варьируемых параметров ЭЛАН на параметры микроструктуры 
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Возможно, это явление объясняется наличием наноразмерного эффекта в упрочненном слое. Присут- 
ствие нанокристаллической структуры в напыленном слое напрямую подтверждается вычисленными значени- 
ями размеров областей когерентного рассеивания (ОКР). Наличие значимых величин угла разориентирования 
структуры кристаллов подтверждается величиной ОКР для исследованного рефлекса 110, а также заметно для 
рефлекса 220, что подкрепляется высоким процентным содержанием аморфной фазы. Похожие явления могут 
быть вызваны также высокими значениями углов разориентирования кристаллитов относительно друг друга. 

Подводя итоги, можно констатировать, что все вышеперечисленные явления (в том числе размеры об- 
ластей когерентного рассеивания) прямо подтверждают наличие наноразмерного эффекта [1]. 

Обсуждение и заключения. Итак, метод электроакустического напыления перспективен для получе- 
ния нанокристаллических структур в поверхностном и приповерхностном слоях напыленных образцов. Изме- 
нение технологических параметров ЭЛАН приводит к появлению указанного в работе разброса значений пара- 
метров кристаллической решетки и областей когерентного рассеивания без определенной выявленной тенден- 
ции. В дальнейшем указанная проблема будет решена. Этому может способствовать экспериментальное опре- 
деление оптимальных режимов напыления, обеспечивающих образование нанокристаллических структур. Кро- 


ме того, структуру напыляемого слоя можно наблюдать и оценивать с ПОМОЩЬЮ электронной микроскопии. 
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Введение. Обоснованы рациональные параметры и 
режимы наклонного шнекового смесителя периодического 
действия для достижения наименьшей энергоемкости 
приготовления кормосмеси при условии соблюдения 
зоотехнических требований к качеству приготовляемых 
кормов по неравномерности смешивания. Установление 
функциональных зависимостей между параметрами и 
режимами позволяет проектировать энергоэффективное 
оборудование для внутрихозяйственного комбикормового 
производства. 

Материалы и методы. Экспериментальные исследования 
процесса приготовления кормосмеси проводили на 
наклонном одновальном шнековом смесителе кормов 
периодического действия. План эксперимента включал 
варьирование четырьмя независимыми факторами: 
частотой вращения вала смесителя, коэффициентом 
заполнения камеры смесителя, продолжительностью 
смешивания и углом наклона смесительной камеры. В 
качестве критериев оптимизации, характеризующих 
эффективность смешивания, были приняты 
неравномерность смешивания и удельная энергоемкость 
процесса. 

Результаты исследования. Построены графические 
зависимости критериев оптимизации от уровня 
варьирования факторов, представляющие собой 
двумерные сечения поверхностей отклика второго 
порядка. Рациональные значения при неравномерности 
смешивания менее 5% составили: частота вращения вала 
— 27,5-36,5 мин |, коэффициент заполнения камеры 
смесителя — 0,43-0,51, продолжительность смешивания 
— 3,0-4,2 мин.; угол наклона смесительной камеры — 
227—25°. При таких значениях параметров неоднородность 
смешивания будет минимальной и составит 4,10—4,18%, а 
энергоёмкость процесса составила от 2,08 до 2,16 кВт:ч/т 
Обсуждение и заключения. Зависимости, полученные в 
результате экспериментальных исследований, позволили 
установить области рациональных конструктивных 
параметров и режимов наклонного  одновального 
шнекового смесителя периодического действия. 
Полученные результаты могут быть использованы в 
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дальнейших исследованиях при разработке исходных 
требований на создание новых технических средств с 
гравитационным эффектом интенсивного смешивания. 
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Введение. Актуальной задачей при проектировании нового кормоприготовительного оборудования 
является обеспечение более низкой себестоимости приготовления кормов [1-3]. Для этого необходимо снизить 
энергобёмкость процессов при сохранении необходимого качества комбикорма. Качественное приготовление 
кормосмеси наиболее целесообразно выполнять смесителем периодического (порционного) действия [4, 5]. 
Среди них широкое распространение получили смесители с лопастными и шнековыми рабочими органами [6]. 
Это обусловлено тем, что состав комбикорма характерен наличием сухого измельчённого зерна (более 90%), не 
изменяющего своих свойств в процессе смешивания с добавками. То есть реологические свойства смеси можно 
считать неизменными в течение всего опыта [7]. 

Одной из перспективных моделей циркуляционных смесителей является модель с наклонным 
бункером, внутри которого установлены один или два шнековых рабочих органа. При этом снижение 
энергоемкости процесса достигается вследствие того, что принудительная подача смешиваемого материала на 
верхнем горизонтальном уровне заменяется гравитационным его осыпанием из верхней части бункера [8]. В 
отличии от горизонтальных смесителей, фактическое отсутствие «придавливающих» слоёв корма в верхней 
части шнека значительно улучшает диффузионное перемешивания продукта [9]. Интенсивное движение 
материала и выровненный гранулометрический состав сырья, при применении современных дробилок [10, 11], 
позволяет минимизировать возникающие при осыпании эффекты сегрегации [12]. 

Материалы и методы. Экспериментальные исследования проводили на наклонном одновальном 
шнековом смесителе кормов периодического действия АКМ-3, разработанном в АНЦ «Донской». Смеситель 
кормов АКМ-3 (рис. 1) вместимостью 2,5 м’ предназначен для получения однородной смеси из сухих 
измельченных компонентов (зерновое, мучнистое, белково-минеральное сырье) при порционном 
(периодическом) режиме работы. Он состоит из рамы, корпуса со смесительной камерой, в которой размещен 
вал со спиральной ленточной противоточной навивкой (шнек), и лопастями в верхней части бункера. 


Предусмотрена возможность установки корпуса на тензометрические датчики. 





Рис. 1. Общий вид наклонного одношнекового порционного смесителя 
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В таком смесителе при вращении вала перемещение групп частиц материала из одного места в другое 
(конвективное смешивание) с наклонным бункером осуществляется в процессе транспортирования его шнеком 
в верхнюю часть наклонного бункера и гравитационного осыпания материала в нижнюю часть бункера. 
Поэтому основным преимуществом наклонного смесителя, по сравнению с вертикальными и горизонтальными, 
является отсутствие мёртвых зон между рабочими органами и стенками смесительной камеры [13]. 

В качестве критериев, характеризующих эффективность смешивания, были приняты неравномерность 
распределения контрольного ингредиента в смеси и удельная энергоемкость процесса. В качестве основного 
ингредиента (наполнителя) смеси в экспериментах использовали измельченную пшеницу влажностью не более 
15% с насыпной плотностью 750 кг/м”. В результате приготовили двухкомпонентную смесь в составе: 
измельченная пшеница — 99%, измельченная поваренная соль — 1% (по массе). 

В качестве контрольных ингредиентов для определения неравномерности смешивания использовали 
измельченную поваренную соль. Её распределение определяли, отбирая по 20 проб массой 50 г каждая из 
различных точек смесительной камеры, после завершения смешивания [14]. Отбор конкретной пробы 
производился по ГОСТ 13 496.0-80. Содержание сухого контрольного компонента (поваренной соли) в пробах 
определяли ионометрическим методом по ГОСТ 13 496.1-98. В качестве показателя неравномерности 
смешивания использовали коэффициент вариации фактического распределения контрольного компонента в 
пробах у, % (у1). 

Удельную энергоемкость процесса смешивания №,„(у2) определяли как общую энергоемкость, 
отнесенную к массе кормосмеси, измеряемую посредством трехфазного электрического счетчика [15]. 

Частоту вращения вала смесителя изменяли частотным преобразователем РЕГТАУЕЯа-075Е путем 
изменения частоты тока асинхронного электродвигателя. 

При проведении эксперимента был реализован трёхуровневый план второго порядка. Опыты при 
реализации плана проводили в трехкратной повторности [16,17]. План проведения экспериментов 
предусматривал варьирование четырьмя независимыми факторами, влияющими на процесс смешивания 
(таблица 1). 

Установлено, что значимость показателя неравномерности смешивания значительно выше, чем у 
последнего фактора, который должен рассматриваться как вспомогательный. 


Таблица | 
Факторы и уровни их варьирования 


Ф Кодированное Диапазон Интервал 
акто 
ы обозначение варьирования варьирования 
Е 


Коэффициент заполнения камеры смесителя 


Результаты исследования. После обработки результатов экспериментов были получены уравнения 





регрессии в закодированном виде и определены соответствующие коэффициенты корреляции: 

у|= 4,06 — 0,96х, + 0,29х› + 0,16хз — 0,88х4 — 0,14х1х› — 0,19х1хз +0,54жхл + 0,28хохз — 0,39хохд — 0,20 хзхд + 
0,56х°, + 0,27х*› + 0,3473 + 0,64х°, 

коэффициент корреляции К! = 0,9607; 

у2= 2,26 + 0,99, + 0,42х› — 0,07хз — 0,28х4 + 0,22х\х›— 0,15х хз +0,35хх4 + 0,1Зхохз + 0,47 хохд — 0,08 хзх4 + 
0,73%“, + +0,37х^› + 0,14х^з + 0,84х°, 

коэффициент корреляции К. = 0,9781. 
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Адекватность модели оценивали по критерию Фишера. Табличное значение критерия Фишера с 


410 уровнем значимости Д=0,05 равно Р=2,8. Значение критерия Фишера в моделях: неоднородность смешивания 
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Е=2,51; энергоемкость процесса смешивания Ё=2,65. Сравнение расчетных значений с табличными показало, 
что расчетные их значения меньше табличных. Следовательно, регрессионные модели адекватно описывают 
исследуемый процесс. Экспериментальные значения критерия Кохрена не превышают табличных. Дисперсии 
являются однородными. 

Переходя от кодированных значений факторов (Х\, Х5, Хз, Х4,) к натуральным (и, К, 1,Ф), получили 
зависимости показателей неравномерности смешивания (9, %) и энергоемкости процесса смешивания (М, 
кВтч) от основных факторов в следующем виде: 

— неравномербность смешивания: 

д = 3,76 + 11и- 0,15К — 1,021- 2,35ф + 16ик — 21: -— 59пф — 0,0641 + 0,98Кф + 2,35 + 2,15и° + 4,33К? + 
3,651” +0,75Г; 

— энергоемкость процесса смешивания: 

М= 2,19 +23п +4,13К + 6,25ф— 13и^—З.Аиф— 5,23 + 0,997 + 0,48 + 0,39%. 

В результате обработки экспериментальных данных построены графические зависимости критериев 
оптимизации от уровня варьирования факторов, представляющие собой совмещенные двумерные сечения 
поверхностей отклика второго порядка. 

На рис. 3 и 4 показаны некоторые двумерные поверхности отклика влияния факторов на процесс 
смешивания. В частности, показаны зависимости влияния частоты вращения вала и коэффициента заполнения 


смесительной камеры на неравномерность смешивания и энергоемкость процесса. 
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Рис. 3. Влияние частоты вращения вала и коэффициента заполнения смесительной камеры на неравномерность смешивания 
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Рис. 4. Влияние частоты вращения вала и коэффициента заполнения смесительной камеры 
на энергоёмкость процесса смешивания 
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В результате проведенных исследований определены рациональные значения параметров и режимов 
работы наклонного одновального порционного смесителя вместимостью 2,5 м” со шнековым рабочим органом 
при приготовлении кормосмесей, отвечающие зоотехническим требованиям по неравномерности смешивания 
при условии минимальной энергоемкости процесса. Рациональными значениями для достижения 
неравномерности смешивания не более 5% являются: частота вращения вала и= 27,5-36,5 мин"; коэффициент 
заполнения смесительной камеры К = 0,43-0,51; продолжительность смешивания 1= 3—4,2 мин; угол наклона 
смесительной камеры ф = 22—25°. При этом энергобёмкость процесса составила от 2,08 до 2,16 кВт-ч/т. 

Помимо этого, в результате экспериментальных исследований установлена наименьшая 
неоднородность смешивания 3,2 % при следующих режимах и параметрах смесителя: частота вращения вала 
п=35 мин"; коэффициент заполнения смесительной камеры ^=0,5; продолжительность смешивания #=3 мин; 
угол наклона смесительной камеры ф=30°. При этом энергоёмкость процесса составляла 2,2 кВт:ч/Т. 
Наименьшая энергоёмкость 2,1 кВт-ч/т при неоднородности смешивания 5% достигается при следующих 
значениях параметров и режимов работы смесителя: частота вращения вала и=30 мин"; коэффициент 
заполнения смесительной камеры К=0,55; продолжительность смешивания [= 3 мин.; угол наклона 
смесительной камеры Фф=25°. 

Обсуждение и заключения. Зависимости, полученные в результате экспериментальных исследований, 
позволили установить области рациональных конструктивных параметров и режимов наклонного одновального 
шнекового смесителя периодического действия. Полученные результаты могут быть использованы при 
разработке технических средств для производства комбикормов, а также для модернизации существующих 
комбикормовых производств с целью увеличения их энергоэффективности. 
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Совершенствование методов контроля процессов коррозии на объектах 
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Введение. Цель работы — повышение надёжности методов 
прогнозирования поведения гидроксида лития в пароводя- 
ном тракте на ТЭС и АЭС и оперативного контроля водо- 
родного показателя паровых растворов. Разработан метод 
оперативного контроля на основе кондуктометрических 
измерений водородного показателя паровых растворов ин- 
гибитора коррозии конструкционных материалов гидрокси- 
да лития, применяемого на ТЭС и АЭС. 

Материалы и методы. Для практической реализации высо- 
котемпературного оперативного контроля водородного по- 
казателя паровых растворов использован метод математиче- 
ского моделирования. 

Результаты исследования. Разработан метод контроля во- 
дородного показателя паровых растворов гидроксида лития, 
основанный на определении концентрации в паре путём 
конденсации пара в охлаждаемом кондуктометрическом 
датчике, размещённом в паровом пространстве парогенера- 
тора, что позволило значительно повысить точность опреде- 
ления концентрации гидроксида лития. Предложены урав- 
нения, описывающие изменение предельной молярной экви- 
валентной электропроводности и констант диссоциации 
гидроксида лития в широком диапазоне изменения парамет- 
ров состояния. 

Обсуждение и заключения. Предложенные метод оператив- 
ного контроля водородного показателя паровых растворов и 
математические модели для расчёта предельной эквивалент- 
ной молярной электропроводности и констант диссоциации 
обеспечивают приемлемую для практических расчётов по- 
грешность и возможность автоматизации измерений. При 
увеличении температуры пара до 573,15 К возникает необ- 
ходимость в увеличении концентрации гидроксида лития в 
паре до 10`* моль/Кг. 


Ключевые слова: машиностроение, коррозионные процес- 
сы, автоматический контроль, паровая фаза, кондуктометри- 
ческие методы, мониторинг. 
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Введение. При организации и ведении водно-химических режимов (ВХР) на ТЭС и АЭС одной из 
главных задач является сведение к минимуму коррозионных повреждений теплосилового оборудования в зоне 
контакта с водой и паром [1]. Дозирование в питательную воду ОН в необходимых количествах обеспечивает 
создание на поверхности металла прочной защитной литий-ферритной плёнки, замедляющей коррозию и кор- 
розионное растрескивание теплообменных труб парогенераторов АЭС с ВВЭР-1000 [2, 3], а также значитель- 
ное увеличение времени наработки оборудования на отказ [3]. При дозировании [ДОН в питательную воду ба- 
рабанных котлов ТГ-104 на ТЭС также наблюдается формирование защитной литий-ферритной плёнки [4]. На 
АЭС ЦОН применяют в России и за рубежом [4—6]. В перспективных атомных реакторах ТУ поколения 5СУЮ 
(бирегсииса[ исиет-соо[еА теасюг), обеспечивающих увеличение термического КПД по сравнению с существу- 
ющими (РУ) с 33 до 44 %, при давлении Р=25 МПа температуру Т водяного пара планируют увеличить до 
953 К [6]. Для создания на ТЭС и АЭС оперативного контроля значений рН водных и паровых растворов МОН, 
определяющих эффективность защиты от коррозии, необходимы данные о предельной эквивалентной электро- 
проводности Л. и константах диссоциации Ка МОН при высоких Р и Т, представленные в виде, удобном для 
практических расчётов. Методика расчёта рН растворов [4ОН на основании данных об удельной электропро- 
водности ), предложенная авторами для жидкой фазы [7], неприемлема для паровой фазы при малых значениях 
плотности р и концентрации С в связи с высокой погрешностью определения хи Ка [8]. 

Целью настоящей работы является повышение надёжности методов прогнозирования поведения [4ОН 
в пароводяном тракте на ТЭС и АЭС и оперативного контроля рН паровых растворов. 

Для достижения поставленной цели созданы уравнения для расчёта Л. и Ка на линии насыщения и в 
перегретом паре при Т до 773,15 К, а также разработана методика оперативного контроля рН паровых раство- 
ров [4ОН. 

Материалы и методы. По сравнению с традиционным контролем, связанным с отбором и охлаждени- 
ем проб, оперативный контроль рН при помощи кондуктометрических датчиков, установленных в теплоноси- 
теле парогенераторов и работающих при Ти Р, близких к рабочим, является более объективным и практически 
безынерционным [7, 8]. В расчётах рН для жидкой фазы в [7] были использованы экспериментальные данные о 
Ка МОН [8], ионном произведении К, [9], и диэлектрической проницаемости & для воды [10]. Определение кон- 
центрации С водных растворов [4ОН осуществлялось на основании зависимостей С МОНЕ} (хх МОН)т-сопзь ПО- 
лученных на основании экспериментальных данных о % растворов. 

При измерении Х растворов в паровой фазе (особенно при малых ри С), когда полярные свойства рас- 
творителя ослаблены, степень диссоциации [4ОН мала и сопротивление межэлектродного пространства изме- 
рительной ячейки кондуктометрического преобразователя может стать соизмеримым с сопротивлением элек- 
трического изолятора, что приведёт к большим погрешностям определения Хх [8]. 

Для практической реализации высокотемпературного оперативного контроля рН необходимые данные 
о Л. и К. могут быть получены в исследуемом диапазоне изменения параметров состояния методом математи- 
ческого моделирования и представлены в виде уравнений. 

Результаты исследования. На основании анализа существующих методов в данной работе предложен 
метод контроля рН паровых растворов /[4ОН, основанный на определении С ШОН в паре путём конденсации 
пара в кондуктометрическом датчике с охлаждаемым капилляром [11], размещённом в паровом пространстве 
парогенератора в потоке теплоносителя. В этом случае С МОН в паре равна С [ЛОН в конденсате пара и может 
быть определена из зависимостей С МОН=Г (Хх МОН)т-сопзь полученных для жидкой фазы [7]. Конденсация пара 
позволяет значительно повысить точность определения 9 и С [8]. Вычисления рН для парорастворов МОН про- 
изводились по методике, описанной в [12], с использованием данных [8] 0 К... 

С целью оценки эффективности применения [4ОН для защиты от коррозии в атомных реакторах ТУ по- 
коления 5СУ/К авторами [6] на основании экспериментальных данных о К/ [13—15] выполнен расчёт рН водных 


Машино строение и машиноведение 


и паровых растворов [4ОН при Р=25 МПа в широком диапазоне изменения р, Т, С и предложено эмпирическое 
уравнение для вычисления моляльной константы ассоциации К. [16] в диапазоне изменения р от 50 до 1000 
кг/м” и Тот 373 до 873 К: 

[0810К. = а! + а>/ + + (аз + а4/т)-10 8106 - (а + аб/т)-(10°6)}, 15 
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где д=р/рс; т=Т/То; ре, Тс — плотность и температура раствора в критической точке; ат, а, аз, ад, а5, ав — кон- 
станты. 

В перегретом паре с уменьшением ри Г К. [4ОН уменьшается, уменьшается и разность между рН НО 
и РН парорастворов [4ОН, что снижает эффективность защиты от коррозии элементов конструкций при высо- 
ких Ри Т. Авторами вычислены значения рН парорастворов [4ОН при Р=25 Мпа: Т=673, 723 и 773 К для С=10°, 
10, 10", 10? и 10 * молыкг. В расчётах были использованы экспериментальные данные о К. [4ОН [8], ионном 
произведении К, [9], и диэлектрической проницаемости воды & [10]. На рис. | представлена зависимость рН 
водных и паровых растворов /4ОН от р при Р=25 МПа для С=10 °-10?, моль/кг. 

В области р<330 кг/м” авторы [13-15] не выполняли экспериментальные исследования, поэтому на ри- 
сунке представлены только результаты расчёта рН ОН с использованием данных о К. МОН, полученных по 
уравнению авторов [16 ]. 


773 ТК 
723 , 
873|| 673 623 573523473 373 





100 200 300 400 500 600 700 800 900 1000 
р, кг/м3 


Рис.1. Зависимость рН водных и паровых растворов 4ОН от плотности р при Р=25 МПа: 
1 — #0; 2, 3, 4, 5, 6 — С МОН=10®, 10°, 10%, 10°, 107 молыкг; 
А,о — расчёт рН по данным о Ка [13-15] и [8]; 
линии 2, 3, 4, 5, 6 — расчёт рН по данным о К, полученным из уравнения [16] 





Полученные авторами ранее на основании измерений электропроводности высокотемпературных рас- 
творов Г4ОН, значения Лои Ка [8], необходимые для создания уравнений, описывающих изменение Ло и Ка при 
изменении р и Т, приведены в таблицах 1 и 2 в сопоставлении с данными других авторов. Для создания уравне- 
ний, описывающих поведение Л. на линии насыщения, с целью уменьшения погрешности аппроксимации дан- 
ных диапазон изменения р был разбит на 2 части (от 100 до 500 и от 500 до 958 кг/м”), для каждой из которых 
предложено уравнение вида: 

Ло=Бо:р” + Вр” + Бор + Ь, 
где Бь, Бу, Б>, Бз — определяемые константы; р — плотность раствора, кг/м”. Значения Бо, бл, Б>, Бз для р=(100— 
500) и (500-958) кг/м” равны —0,000004; 0,005878; —2,810820; 1942,943226 и —1,25711:10?; 0,024177411; — 
16,58559801; 5334,222549 соответственно. Для значений рК.=0е1о Ка вычислены значения коэффициентов 
1—4 в уравнении, аналогичном представленному в [16] и приведённому выше: 
рКа=а! + 4 + (4 + а/®)-10910б + (45 + а5/®)- (1085), 

где д=р/р‹; т=Т/То; ре, Т. — плотность и температура раствора в критической точке; 41, 4, 43, 44, 45, ав — кон- 
станты. 

При вычислении констант 4—4 были использованы полученные авторами ранее экспериментальные 
данные о К. [8]. На линии насыщения авторы учитывали только свои экспериментальные данные. Расчёты для 
перегретого пара при Т=673,15 и 773,15 К выполнены в двух вариантах: в первом рассматривались только ав- 


торские данные о К\ [8], во втором — все данные, приведённые в табл. 2. 
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Таблица 1 
Зависимость предельной молярной электропроводности Л. [4ОН 
от плотности р на линии насыщения и в перегретом паре 


р, кг/м” ИА 10: См-м”/кмоль 


Ами 
Е ОИ ПОС СР ООО ООО 
[пя [вв | вЕВТ 
ООС И ПРЕДО 0 ООО ООО О 
ОИ ЛЬ ОЗЕСЕ 2 ООО ОО О 
ПОС ПОС ЕОЬООС ООО ОЯ ООИ 
ООО ПО ОО О ПО 
И м 
С О ОИ ПО ОИ 
ОС ОЕ ООО ПО ОО 

Для линии насыщения: 4=119,50984; 4=-115,33347; 4:=—174,20233; а.=167,90533; &=—143,78368; 
@4=183,09356. 

Для первого варианта расчёта для перегретого пара: 4=7,6323975; 4=-3,214114; 4.=10,9289; 
44=2,5463453; @=-—10,34132; ф=15,190972; для второго варианта расчёта 4=10,113404; 4=—-6,133638; @&=- 
17,85984; 41=8,888042; 45=8,405401; 4=14,648395. 

Отклонения в % экспериментальных значений Ло и рК( от вычисленных по предложенным уравнениям 
представлены в табл. 3. 

Таблица 2 


Зависимость отрицательного логарифма констант диссоциации ОН рКа 
от плотности р на линии насыщения и в перегретом паре 
рРКа, молыкг 


Км 
о 
о 
о 
о 
О О ЕТО ООО 
| З —_ чвьзошя —_ 
0 

0 , 
2 


р 

3 

3 
5 

РЖ 


3,86[8] 3,98[8]; 2,82[16]; 3,14[13] 
3,1115] 
2,1313] 2,25[13] 


А ИЕ 
3,3818]; 2,81[13] 3,4318]; 3,1513] 


И ОТО ООО ПО 
р р м 
м 80 |  ЗЕОЯЯ 


1,4115] (= 1.02 МНа) 


БО КРЕМА | 


1,0315] (Ф=4,88 МНа) 


0 
3 
6 
0 ; 5, 
0 `_ 

И 

0 4,3118]; 4,0116] 
0 ) . ) 

0 

1 

0 ем, 
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Таблица 3 
Отклонение (в %) экспериментальных значений предельной эквивалентной 
электропроводности (0Л.) и отрицательного логарифма моляльных констант диссоциации (брКа) 
от значений, вычисленных по уравнениям 
ОЛ. ОрКа 
Т=673,15 К 1=773,15 К 





1=673,15 К 1=773,15 К 
Линия Линия 


1 вариант [вариант 2 вариант 2 вариант 


насыщения насыщения 





расчета расчета расчета расчета 


0,36[8] 0,04[8] 1,20[8] 1,36[8] 
5,86 [16] 15,5[16] 
Боя | сям розни пин ея | яя _— 


0,52[8]; 5,95[8]; 
200 0,34[8] 0,11[8] 0,92[8] 0,39[8] 
0,98[16] 10,5[16] 
3,86[8]; 1,9918]; 
0,09[8] 2,50[8] 0,29[8] 
9, вия — 5,96[16] 


14 и — 
6,44[8]; 
0,03[8] 0,12[8] 4,99[8] 4,07[8] 20,2[16]; 
0,56[16]; 
7,9613] 


Е О О И По И 


500 0,09[8] 0,16[8] 8,95[8] 2,40[8] 
5,90[13] 0,67[13] 
13,72[8]}; 5,24[8]; 
0,10[8] 0,177[8] 4,56[8] 6,11[8] 8,89[13]; 0,6013]; 
67,8[16] 24,1[16] 


ЕЕ зан 
0,44 [8] 0,05[8] 7,82[8] 8,0313] 
5,95[13] 


О О О ПО С С ПОЕТ 
ПС С О ОИ ПО МЕС ОС 
О Е О О О О И 
ПС СО О О О ООН О 


Обсуждение и заключения. Как очевидно из рис. |, полученные авторами значения рН МОН при 
Т=673 К ниже по сравнению с полученными путём расчёта [6] в среднем на 0,2; при 1=723 К выше на 0,1-0,2; 

















при 1=773 К выше на 0,1-0,7 единиц рН. 

Таким образом, величина рассогласования данных о рН МОН в перегретом паре, полученных авторами 
и изложенных в [6], возрастает с увеличением Ги С. Такое рассогласование можно считать удовлетворитель- 
ным, если учесть, что с уменьшением р до 200-100 кг/м” погрешность экспериментальных данных [8] о Ка 4ОН 
возрастает до 180 % и, кроме того, при составлении уравнения авторами [6] не использовались эксперимен- 
тальные данные о К„ при р<330 кг/м”. 

Отклонение экспериментальных значений Л. и рКа от вычисленных по предложенным в настоящей ра- 
боте уравнениям на линии насыщения не превышает 1 % и только при р=917 и р=958 кг/м” достигает 2,27 и 
1,51 % для Л. (табл. 3). В перегретом паре при Т=673,15 и Т=773,15 К для первого варианта расчёта с использо- 
ванием только авторских экспериментальных данных о К„ [8] наибольшее отклонение экспериментальных дан- 
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ных рК. от расчётных не превышает 2,5 % при р=100-300 кг/м” и возрастает до 8,95 % с увеличением р до 400- 
700 кг/м”. Для второго варианта расчёта с учётом данных всех авторов, представленных в табл. 2, наибольшая 
погрешность расчёта по уравнению при р=100-300 кг/м” и 7=673,15 К; Т=773,15 К достигает 15,5 % для данных 
авторов [16], а при р=100-300 кг/м” и Т=673,15 К; Т=773,15 К — 67,8 %. 

На основании анализа полученных результатов можно сделать вывод о том, что использование /[4ОН в 
качестве корректирующей добавки для регулирования рН с целью предотвращения коррозионного разрушения 
элементов конструкций перспективных атомных реакторов ГТУ поколения 5СУК не вызывает затруднений при 
1<673,15 К и Р=25 МПа. С ростом температуры уже при 1=773,15 К для эффективной защиты от коррозии 
необходимо увеличивать С [4ОН до 10” молыкг, что выдвигает задачу исследования растворимости [4ОН в 
паре при этих параметрах состояния. 

Предложенные в настоящей работе метод оперативного контроля рН паровых растворов [4ОН и мате- 
матические модели для расчёта Л. и Кл обеспечивают возможность определения рН при помощи измерительно- 
вычислительного комплекса, работающего непрерывно в автоматическом режиме, имеющего малое время за- 
паздывания сигнала и обеспечивающего приемлемую для практических расчётов погрешность результатов из- 
мерения. Использование таких устройств в современных системах химико-технологического мониторинга на 
ТЭС и АЭС позволит значительно повысить надёжность методов прогнозирования поведения [4ОН в парово- 
дяном тракте и оперативного контроля рН. 
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Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация 


Раа теазигетепе зузет Гог поп-ЧезгисНуе диаеу 1е5Нпе оЁ Вага аПоу5 ^^ 


Г. К. Твубчу', У. Г.. УаНКоу”, У. 1. епаепКо? ” 


1.3 Роп $(ае Тесытиса1 Отгуегзиу, Козюу-оп-Роп, Кизз1ап Недегайоп 


Введение. Известно, что качество изделий из спеченных 
твердых сплавов, при отсутствии нарушений технологии 
изготовления, определяется фазовым составом и средним 
размером зерен карбидной фазы. Однако твердые сплавы 
обладают недостатком, присущим в той или иной степени 


всем изделиям порошковой 


структуры и 
вариацией свойств. Традиционные методы контроля 
структуры и фазового 
выборочных разрушающих испытаний, не позволяют 


металлургии — 
неоднородностью соответствующей 
состава по результатам 


оперативно в условиях производства получать 
информацию о качестве и не гарантируют идентичности 
свойств всех изделий исследуемой партии. Основным 
методом неразрушающего контроля качества твердых 
сплавов является измерение коэрцитивной силы. Однако 
отечественных коэрцитиметров в настоящее время не 
производится. В связи с этим целью настоящей работы 
является создание отечественного ’информационно- 
измерительного комплекса для неразрушающего контроля 
качества изделий из спеченных твердых сплавов и оценка 
достоверности полученных с его помощью результатов. 
Материалы и методы. В работе использованы 
цилиндрические и сфероконические стержни диаметром 
от 4 до 10 мм, изготовленные из сплавов марок ВК8 и 
ВК!10ХОМ. Для определения коэрцитивной силы в 
информационно-измерительном комплексе использован 
метод сравнения с аттестованным изделием. 

Результаты исследования. Проведены — измерения 
коэрцитивной силы двух партий изделий из твердых 
сплавов марки ВК1ОХОМ и марки ВК8, выполненных с 
помощью разработанного информационно- 
измерительного комплекса и коэрцитиметра КОЕВУЛМАТ 
1.097 Нсу. 
Обсуждение и заключения. Проведенные 
экспериментальные исследования показали, что значения 
коэрцитивной силы образцов, полученные на разном 


оборудовании, имеют примерно одинаковый уровень. 


* “ 
Работа выполнена в рамках инициативной НИР. 
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Подтверждена, как следствие технологии производства, 


характерная для спеченных — твердых сплавов 
неоднородность структуры. Учитывая, что относительная 
погрешность измерений коэрцитивной силы для каждого 
образца должна находиться в пределах + 6 %, можно 
считать полученные результаты удовлетворительными. 
Показано, что предложенный и реализованный в 
информационно-измерительном комплексе алгоритм 
измерения путем сравнения характеристик испытуемых и 
аттестованного образца с известной коэрцитивной силой, 
позволяет производителю изделий из твердых сплавов 
расширить выборку контролируемых изделий вплоть до 
100% и существенно снизить затраты на неразрушающий 


контроль продукции. 


Ключевые слова: твердые сплавы, неразрушающий 
контроль качества, коэрцитивная сила, информационно- 
измерительный комплекс. 
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Введение. Среди материалов, получаемых методами порошковой металлургии, широкое 
распространение получили твердые сплавы, базовый состав которых представляет собой мелкодисперсные 
карбиды вольфрама и титана на кобальтовой связке. Применение изделий из твердых сплавов обусловлено 
рядом ценных свойств, таких как твердость, сопротивление износу при трении, износостойкость. В то же время 
нарушение технологии производства порошков и особенностей спекания вызывают появление структурной 
неоднородности внутри одной партии готовых изделий и, как следствие, неоднородности эксплуатационных 
характеристик в широких пределах. Поэтому выборочный контроль нескольких твердосплавных изделий из 
партии с применением, например, металлографического анализа или механических испытаний не дает гарантии 
идентичности структуры и свойств всей партии. 

Известно, что качество изделий из спеченных твердых сплавов определяется, главным образом, 
фазовым составом, т.е. соотношением количества ферромагнитной кобальтовой фазы и парамагнитных 
карбидов и средним размером зерен карбидной фазы [1-3]. При фиксированном содержании кобальтовой фазы 
для каждой марки твердосплавного материала прочность готовых изделий практически полностью 
определяется зернистостью карбидов [4]. 

Альтернативой разрушающим испытаниям является неразрушающий контроль качества твердых 
сплавов по коэрцитивной силе, которая в первом приближении обратно-пропорционально зависит от 
зернистости карбидов. Значение коэрцитивной силы, наряду с прочностными характеристиками твердых 
сплавов, введено в действующие до настоящего времени стандарты [5, 6]. Однако отсутствие в настоящее 
время производства отечественных коэрцитиметров вынуждает производителей изделий из твердых сплавов 
применять для неразрушающего контроля зарубежные аналоги. 

В связи с этим, целью настоящей работы являлось создание информационно-измерительного 
комплекса для неразрушающего контроля качества изделий из спеченных твердых сплавов и оценка 
достоверности полученных с его помощью результатов контроля путем сравнения с результатами, 
полученными при использовании коэрцитиметра КОЕВЛМАТ 1.097 Нс) фирмы «зави Ог.Еоегуег ОшЬН & 
Со. КО» [7], включенного в Государственный реестр РФ средств измерений [8]. Результаты измерения 
коэрцитивной силы на коэрцитиметре КОЕКУЛМАТ 1.097 Не] предоставлены ОАО «Серпуховский 
инструментальный завод «ТВИНТОС». 

Материалы и методы. Образцы для исследований представляли собой цилиндрические и 
сфероконические стержни диаметром от 4 до 10 мм, изготовленные из сплавов марок ВК& и ВК10ХОМ. 

Разработанный информационно-измерительный комплекс [9,10], структурная схема которого 


представлена на рис. 1, состоит из шести функциональных блоков. 


Цыбрий И. К. и др. Информационно-измерительный комплекс для неразрушающего контроля качества твердых сплавов 


ТзуБту [. К. апа йе ошету. аа теа5игетепЕ ;убет ог поп-аейгисйНуе диаШу (езИпэ о] ага аПоу5 


Измеритель- Устройство Устройство 


ный блок масштабирования математической обработки 


Стабилизатор 


напряжения Блок питания и Органы управления и 


ВЕ Я индикации лицевой панели 


2208 50Гц автоматики 





Рис.1. Структурная схема информационно-измерительного комплекса 


Стабилизатор напряжения собран на повышающем трансформаторе и предназначен для увеличения 
напряжения питания до 300 В и подачи его на катушку намагничивания измерительного блока. Включение 
трансформатора производится магнитным пускателем с искрогасящим контуром, что обеспечивает мягкое 
включение и выключение измерительного блока. Управление магнитным пускателем осуществляется 
микроконтроллером через симисторную схему включения. 

Измерительный блок представляет собой соленоид с намагничивающей катушкой и закрепленным 
внутри него измерительным контуром (рис. 2). Параметры катушки намагничивания позволяют получать 
напряженность магнитного поля внутри соленоида до 350 кА/м. Величина магнитного поля, создаваемого в 
соленоиде блока намагничивания, позволяет измерять коэрцитивную силу близко к полю насыщению 
материала. 

Устройство масштабирования предназначено для настройки уровня измеряемого сигнала под 
возможности микроконтроллера. Блок питания и релейной автоматики предназначен для подачи различного 
вида и уровня напряжений на составные части устройств масштабирования, математической обработки и 
органы управления и индикации. 
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Рис. 2. Измерительный блок 


Устройство математической обработки предназначено для вычисления и запоминания уровня сигнала 
соответствующего коэрцитивной силе аттестованного образца и вычисления методом сравнения уровня 
сигнала, соответствующего коэрцитивной силе проверяемого изделия. 

Алгоритм работы информационно-измерительного комплекса следующий: 

— аттестованный образец помещается в измерительный блок, значение коэрцитивной силы выводится 
на индикатор лицевой панели; 

— измеряется и запоминается напряжение 01, соответствующее известному значению коэрцитивной 
силы; 

— аттестованный образец удаляется из измерительного блока; 

— испытуемый образец помещается в измерительный блок; 

— измеряется и запоминается напряжение (5, соответствующее значению коэрцитивной силы 
испытуемого образца; 

— полученное путем сравнения напряжений О\ и > значение коэрцитивной силы испытуемого образца 
выводится на индикатор лицевой панели; 

— испытуемый образец удаляется из измерительного блока. 

При необходимости измерения партии изделий одной марки и одного типоразмера с аттестованным 
образцом процесс измерения повторяется, начиная с пункта 4). 

Результаты исследования. Результаты сравнительного анализа, полученные путем измерения 
значений коэрцитивной силы 24 образцов изделий из твердых сплавов марки ВК10ОХОМ с помощью 
разработанного информационно-измерительного комплекса (ряд 1) и коэрцитиметра КОЕК МАТ 1.097 Нс] 
(ряд 2), представлены на рис. 3. Аналогичные результаты для 46 образцов изделий из твердых сплавов марки 
ВК& представлены на рис. 4. 
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Рис. 3. Коэрцитивная сила образцов марки ВК10ХОМ 
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Рис. 3. Коэрцитивная сила образцов марки ВК8 


Обсуждение. Проведенные экспериментальные исследования показали, что значения коэрцитивной 
силы образцов, полученные на приборе «КОЕК МАТ 1.097 НсТ» и информационно-измерительном комплексе, 
имеют примерно одинаковый уровень. Учитывая, что относительная погрешность измерения коэрцитивной силы 
для коэрцитиметра «КОЕВ/ЛМАТ 1.097 Нс.» находится в пределах = 6 %, можно считать полученные результаты 
сравнительных испытаний удовлетворительными. 

Вариация значений коэрцитивной силы образцов в рамках одной марки до 10-12% может быть вызвана 
следующими причинами: 

— использованием при спекании порошковых смесей разных производителей; 

— несоблюдением требуемых режимов производства; 

— влиянием на показания приборов конфигурации образцов, т.к. пониженные значения коэрцитивной 
силы наблюдались на образцах сфероконической формы, как следствие большего размагничивающего фактора, 
характерного для такого типа образцов. 

Заключение. Таким образом, сравнительные испытания информационно-измерительного комплекса 
показали принципиальную возможность его использования для контроля качества изделий из спеченных 
твердых сплавов, а предложенный и реализованный в комплексе алгоритм измерения путем сравнения 
характеристик испытуемых образцов и аттестованного образца с известной коэрцитивной силой, позволяет 
производителю изделий из твердых сплавов расширить выборку контролируемых изделий вплоть до 100%, 
сократить время на контрольные операции и существенно снизить затраты на неразрушающий контроль 
продукции. 
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Введение. Статья посвящена теоретическому описанию 
процесса переноса ионов в мембранных системах в галь- 
ваностатическом режиме. В качестве мембранной системы 
рассматривается канал обессоливания электродиализного 
аппарата. Цели работы: создание и верификация двумер- 
ной математической модели стационарного переноса 
ионов соли в канале обессоливания электродиализного 
аппарата для гальваностатического режима. 

Материалы и методы. Предложена новая модель перено- 
са ионов. Она основана на системе уравнений Нернста — 
Планка — Пуассона для электрического потенциала и на 
уравнении для функции электрического тока. Получено 
численное решение краевой задачи модели методом ко- 
нечных элементов с помощью программного пакета 
Сот50о1 МиШрйу51с$. 

Результаты исследования. Разработанная математическая 
модель позволяет описать стационарный перенос ионов 
бинарной соли в канале обессоливания электродиализного 
аппарата. При этом учитываются нарушение электроней- 
тральности раствора и формирование расширенной обла- 
сти пространственного заряда при сверхпредельных токах 
в гальваностатическом режиме. Об адекватности постро- 
енной модели свидетельствует хорошее совпадение физи- 
ко-химических характеристик переноса, рассчитанных по 
моделям для гальваностатического и потенциостатическо- 
го режимов. 

Обсуждение и заключения. Разработанная модель позво- 
ляет интерпретировать результаты экспериментальных 
исследований переноса ионов в мембранных системах, 
если данный процесс протекает в гальваностатическом 
режиме. Некоторые электрокинетические процессы связа- 
ны с появлением расширенной области пространственного 
заряда при сверхпредельных токах. Описывая формирова- 
ние указанной области, можно выяснить, каким образом 
зависящие от нее процессы влияют на перенос ионов в 
гальваностатическом режиме. 
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Введение. Мембранные системы составляют основу электродиализных аппаратов, нано- и микрофлю- 
идных устройств, которые применяются при очистке воды, переработке сельскохозяйственной продукции (мо- 
лока, вина и т. д.), выполнении химических анализов и в других сферах деятельности [1-4]. В многочисленных 
математических моделях процессов массопереноса в мембранных системах для потенциостатического или по- 
тенциодинамического режимов электрический режим задается как скачок потенциала между двумя эквипотен- 
циальными плоскостями, параллельными мембранам. Подробный обзор работ, посвященных моделированию 
для потенциостатического режима, представлен в [5—7]. 

В то же время в практике электродиализа, электрохимической характеризации мембран (хронопотен- 
циометрия, импедансометрия и др.) часто используется гальваностатический режим, при котором на межфаз- 
ной границе поддерживается постоянная средняя плотность тока. Об этом режиме собран огромный объем экс- 
периментальных данных, которые необходимо интерпретировать [8—10]Ошибка! Закладка не определена.. 
Исследования в области математического моделирования гальваностатического режима ведутся по нескольким 
направлениям. 

Первое направление — метод обратной задачи. Как ясно из названия, речь идет о решении обратной 
задачи: для заданной плотности тока на межфазной границе «раствор — мембрана» находится соответствую- 
щий скачок потенциала, а далее рассматривается задача для потенциостатического режима [11]. Низкая эффек- 
тивность данного метода обусловлена тем, что его реализация требует многократного решения задачи в потен- 
циостатическом режиме для одного заданного значения плотности тока. 

Второе направление — метод декомпозиции. При этом система уравнений Нернста — Планка — Пуас- 
сона заменяется системой декомпозиционных уравнений [12—16]. Предположение о квазиравномерном распре- 
делении заряда позволяет получить модель для гальваностатического режима в приближении закона Ома [17— 
20]. 

Третий подход можно назвать прямым методом. В этом случае для плотности тока в канале обессоли- 
вания выводится уравнение, заменяющее уравнение Пуассона [21]. 

Гальваностатический режим можно описать иначе — с помощью численного решения уравнений 
Нернста — Планка — Пуассона для электрического потенциала со специальным граничным условием, которое 
позволяет установить плотность тока как параметр, задающий электрический режим в системе. В [22, 23] для 
одномерного случая производная по времени градиента электрического потенциала определялась как явная 
функция плотности тока. Это отличает подходы авторов указанных работ от потенциостатических моделей, в 
которых задается разность потенциалов. 

В данной статье представлена стационарная модель процесса переноса ионов в мембранных системах 
для гальваностатического режима. Она основана на системе уравнений Нернста — Планка — Пуассона с гра- 
ничным условием, которое позволяет установить плотность тока как параметр, задающий электрический режим 
в системе. В этом решение аналогично [22, 23]. Отличие же состоит в том, что предлагаемая модель двумерна и 
учитывает непостоянство плотности тока по длине канала. 

Материалы и методы. Под мембранной системой подразумевается канал обессоливания электродиа- 
лизного аппарата (ЭДА), образованный анионообменной (АОМ) и катионообменной (КОМ) мембранами. Через 


него со средней скоростью У прокачивается раствор бинарного электролита. На рис. | х — нормальная к по- 


верхности мембраны координата, изменяющаяся от 0 (граница с АОМ) до й (граница с КОМ); у — тангенци- 
альная к поверхности мембраны координата, изменяющаяся от 0 (вход в канал) до [ (выход из канала). 
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Рис. 1. Схема канала обессоливания ЭДА. Показаны концентрационные профили катионов С\ (сплошная линия) и 


анионов С› (пунктирная линия). Скорость вынужденного течения У 


Система уравнении. Рассмотрим двумерный стационарный случай системы уравнений, описывающей 
перенос бинарного электролита в отсутствие химических реакций [24]: 


Е. Е а 

Л 51. С; УФф-Р; УС; + СУ, У ВАА (1) 
КТ 

—@У 7 =0, 1=Ь2, (2) 
508,Аф =-Р(аС1 + 52С.), (3) 
ТЕР +22),). (4) 
Здесь |, О; зи С, — соответственно поток, коэффициент диффузии, зарядовое число и молярная концен- 
трация {1-го иона; ф — электрический потенциал; У — скорость потока раствора электролита; &0 — электриче- 


ская постоянная, &, —— относительная диэлектрическая проницаемость раствора электролита (предполагается 


постоянной); [ — плотность тока; ЕЁ — постоянная Фарадея; К — газовая постоянная; ТГ — абсолютная темпе- 


ратура; Л, /2, 1, Ф,С., С› — неизвестные функции координат х иу. 


Уравнения Нернста — Планка (1) описывают поток ионов, обусловленный миграцией в электрическом 
поле, диффузией и конвекцией; (2) — уравнение материального баланса в стационарном случае; (3) — уравне- 
ние Пуассона для потенциала электрического поля; (4) — плотность тока в растворе электролита. Будем счи- 
тать, что распределение скорости в канале соответствует течению Пуазейля [24]: 


у. =0, У, =6% > а | (5) 


Для определения перечисленных выше неизвестных функций необходимо задать краевые условия для 
системы (1)—3). Рассмотрим два электрических режима: потенциостатический, когда скачок потенциала на 
границах системы задается постоянным, и гальваностатический, когда постоянна плотность тока, протекающе- 
го через межфазную границу. 

Граничные условия для моделирования потенциостатического режима. Будем считать поверхности 
ионообменных мембран эквипотенциальными. Система (1)—(4) включает потенциал электрического поля толь- 
ко в форме производных по пространственной координате. В этом случае существенным является только ска- 


чок потенциала Аф = ф<(й, у) — Ф<(0, у) ‚ где Аф — известная функция, поэтому положим, например: 
<(0, у)=0. (6) 
Тогда 
До = соп5Е . (7) 
Условия (6) и (7) определяют потенциостатический режим. Остальные граничные условия приводятся 
ниже. 


На межфазной границе «АОМ — раствор» (х = 0) концентрация коионов (катионов) определяется из 
условия непрерывности их потока у границы «мембрана — раствор» с учетом селективных свойств АОМ [16]: 
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[7 Е =\ о (1-7>)1, (0,5) (8) 


х ВТ ах эВ 


где Т; (= 1, 2) — эффективные числа переноса противоионов в мембране (КОМ и АОМ соответственно); Т; — 


числа, близкие к 1, причем для идеально селективной мембраны /; = 1, и условие (8) превращается в условие 


непроницаемости мембраны для коионов. 
Концентрация противоионов (анионов) зависит от обменной емкости КОМ, что можно задать в виде: 
С. (0,у) =С.„ =М,С.. (9) 
Здесь постоянная №, показывает, во сколько раз эта концентрация отличается от концентрации в объеме раство- 
ра [25]. 
На межфазной границе «раствор — КОМ» (х = й) для концентраций ионов приняты условия, аналогич- 
ные условиям на границе «АОМ — раствор» (х = 0): 


С(, у) = Ст = № Со, (10) 
ОС Е д 1-7) Г. (®, 
+ м (у) = (27215 >) :) х( У (10) 
Ох ЮТ Ох <о ЕР 
На входе в канал (у = 0) принято равномерное распределение концентрации ионов: 
С: (х.0) = Су, #=12. (12) 


Условие для потенциала электрического поля получено из уравнений (1) и (4) с учетом отсутствия про- 
текания тока через вход /, (х, 0, В) =0: 


Оф(х,О КТ ОС1 (х,0 ОСо (х,О 
оо Е ЗЕ [а РН, а] (13) 
ду ЕГО + &22.5)Со ду ду 
На выходе из канала (у = Г) ионы свободно выносятся потоком раствора: 
. = 5 2 - | 
(п, ]; ) = = 50; С, Уф-Б; УС; +СМ) = СИ, 1=Ь2. (14) 


Система (14) означает также, что сумма диффузионной и миграционной тангенциальных составляю- 
щих потоков анионов и катионов равна 0: 


бе ЙО 
| _ итас, |) 0,5=1,2. (15) 


Для потенциала принято «мягкое» условие, означающее отсутствие резких изменений потенциала на 





выходе из канала: 


= 0. (16) 
ду 


Краевая задача, включающая уравнения (1)—(4) и краевые условия (6)—(16), моделирует потенциостати- 
ческий режим, причем его определяющее ключевое условие — (7). 

Граничные условия для моделирования гальваностатического режима. При моделировании гальвано- 
статического режима условие (6) остается, а (7) необходимо заменить условием, связанным с заданным значе- 


нием средней плотности тока 1, на межфазной границе «раствор — КОМ» (х = Й). Для вывода такого гра- 


ничного условия подставим соотношения (1) в (4) и выразим градиент потенциала электрического поля: 
АТ = - 
Уф = — р р. ([ + Е(аПУС я 520>У С) == Е(аС = 52С5)У). (17) 
Е (3 ВС -Ё 5210№С.) 
Полагая в (17) Хх = ‚, получаем соотношение, связывающее градиент электрического потенциала с за- 


данным значением плотности тока на границе, то есть граничное условие на межфазной границе «раствор — 
КОМ»: 


ас, вс, 
Г; + РО — + Е 5205 Е 
Ви = ВКТ дх ах 


(й, у) (18) 
Ох Е2 2 О:С! ЕН &>02С. 


При этом плотность тока /,. должна удовлетворять условию: 
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1! | 
- [1.4 = 4. (1) 
0 


Математическая модель гальваностатического режима состоит из системы уравнений (1)—(4). Гранич- 
ные условия (18) и (19) заменяют условие (7). Остальные граничные условия совпадают с условиями для по- 
тенциостатического режима. 

Скачок потенциала в гальваностатическом режиме является вычисляемой величиной. 

Преобразование граничных условий для моделирования гальваностатического режима. Условие (19) 
неудобно для численного решения, так как содержит интеграл. Ниже приводится один из вариантов преобразо- 


вания этого условия. 


В стационарном случае плотность тока 1 является соленоидальным вектором. Действительно, если 


умножить (2) на 5; и сложить, то ЧУ/ = 0. Следовательно существует такая функция 1, что 


Е, 


| 20 
9х ” ду ь — 
С использованием функции 1 условия (18) и (19) переписываются следующим образом: 
| ы = + Ес О, ет + Ре,Ь, — 
х х 
ААУ) ===> 2 2 (п, у), (20) 
Ох Е Са. С, 
1! 1-0 1 | 
= [1 =-- [уу = @)-п@0)) =. (22) 
в. Г. ду | 
Уравнение (22) можно переписать: 
(й,1)-п(ф,0) ==. (23) 


Чтобы замкнуть систему выражений, необходимо получить уравнение для функции 1. С этой целью, 


как и в работах [15, 16], введем в рассмотрение линейный дифференциальный оператор, который является 
функцией вихря (ротором) в двумерном случае, для произвольного двумерного вектора У: 
и _ СУ 


И) = —* | 
(и) х в (24) 





Несложно проверить, что: 
1) г(Уи) = 0 для любой гладкой функции И; 
2.) (и) = (Уи,\)1 + иг(’) для любой гладкой функции и и любого гладкого вектора У. 


р ди ди ы 
Здесь (Уи, И)! = = — — И. — кососимметричное скалярное произведение векторов Уи и И ,‚ причем 
х 


У ду 
(а,а)| =0 для любого вектора а. 
Применяя (24) к уравнению для плотности тока (4), получаем: 
"(Г = Еаг(П) + Ее"). (25) 


Используя формулу потоков (1), получим соотношение: 


_ Е ы 
Г(];) —-— ариг(С;Уф) — Б.У С;) + (СУМ), а (26) 
Отсюда с учетом свойств оператора Г: 
- Е - 2 
1: )= и — 5:0; (УС; Уф), + (УС, У) + С;т(У), =1,2. (27) 
- ОГ, 9, 
Учитывая (27) и г(Г) = = = Ди ‚ уравнение (25) можно записать: 
х 
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2 
ыы /Зщ |[ 2р А +12 ор, 0с, \0ф _ ар, А +2, 0, 00 |. 
ВТ ОХ @х / 0 ду ду ты 


дс д 0 Ос, ОУ, _ 0, 
+Е не ИЕ а— “+: У +Р (ас+5с,) | —^ 
Ох Ох ду . бх 0у 


Из (28) и (20) следует, что функция П определяется с точностью до константы, поэтому можно допу- 


(28) 





стить: 
1(й,0) =0. (29) 

Тогда из (23) получаем 

(А, =. (30) 

Условия (29) и (30) — краевые для функции 1. 

Математическая модель для гальваностатического режима в безразмерном виде. Для численного ис- 
следования краевых задач удобно перейти к безразмерному виду. Так можно упростить уравнения и выяснить 
фактическое количество и набор параметров, определяющих поведение системы. Безразмерные переменные 
описывают класс подобных процессов, характеризующихся одинаковым значением безразмерных чисел. 

Характерные величины, описывающие задачу. При моделировании процессов массопереноса в камере 
обессоливания ЭДА принимается ряд характерных значений: 

— для пространственных координат — межмембранное расстояние й; 


— для концентраций ионов — объемная концентрация электролита Со; 

— для скорости — средняя скорость вынужденного течения У, ; 

— для коэффициентов диффузии — коэффициент диффузии электролита р =р.О, ($, -5, )/ Ву, 
— для электрического потенциала — тепловой потенциал Фф, = КТ/ Иа 

— для плотности тока — величина &, = ЕРСу/й (аналог предельной плотности диффузионного тока); 
— для потока ионов — диффузионный поток ®=В9( ‚ й. 


Формулы перехода. Переведем уравнения в безразмерную форму с помощью следующих соотношений 
(индексом (и) обозначены безразмерные варианты величин): 


— 


ох У ЮГ У со 


Хх = у , 1 
й й й \ Со 
1 1 р к 
060 =. =, о р =. 
Фо й ' Л 
Система уравнений в безразмерной форме имеет вид (индекс (И) для упрощения записи опущен): 
1; =-=), С, Уф-Б, УС; + РеСУ, 1=12, (32) 
—ау =0, 1=Ь2, (33) 
дс 
Ат=-|| В, — р, 9% ыы гр ба :2р, 9% 2 |8 |+ 
де сх ду ты - Ох (35 
д д д д 0У, ду 
+Ре ча, — У, -Ре ры НМ У. +Ре(ас + 5,с,)| —^-—^ |, 
Ох Ох а. ду 0х ду 
ТЕ, + 5212 . (36) 
Система уравнений (29)—(35) содержит два безразмерных числа: Ре =Ууй/Рр — число Пекле и 


в й) ыы <“ 2 

&=8 „воВТ/ (Сий” Е”). Физический смысл параметра = состоит в том, что это удвоенный квадрат безразмерной 
ь | у, 

дебаевской длины р: &=2(1р /й)" [5]. Оценка величин параметров показывает, что при естественных для 


электродиализа условиях число Пекле имеет порядок 10? —10°, число & имеет порядок 10! —107” ‚то есть 


может считаться малым Параметром. 
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Для удобства численного решения преобразуем систему уравнений, подставив плотность потока (32) в 


уравнения (33) и (36): 
4’ (-=.0, С, Уф- Б,УС, +РеСУ)=0, #=1,2, (37) 
- 2 и 
ГУ: [ <), С. Уф-Б, УС, + Рес (38) 
= 


Таким образом, система уравнений содержит следующие неизвестные функции х, у: С, С2,Ф, [х,1.. 

Поля концентраций С\, С. и потенциала ф определяются решением уравнений (37), (34) соответственно. Ком- 

поненты плотности тока /,,[, вычисляются с помощью (38). Распределение скорости (5) в безразмерной фор- 
ме: 

У =0, У, =6х (1-х). (39) 


Присутствие малого параметра # в уравнении Пуассона (34) означает, что краевая задача является син- 
гулярно возмущенной. Это значительно усложняет ее численное решение, поскольку такие задачи относятся к 
жестким [26]. Потенциал электрического поля ф и концентрации ионов С, С> изменяются очень быстро в уз- 


ком пограничном слое, толщина которого равна длине Дебая [> [5]. Для решения этой проблемы целесообраз- 
› р 


но уплотнить вычислительную сетку в пограничном слое и использовать специальные методы решения жест- 
ких задач [26]. 
Граничные условия в безразмерной форме. На межфазной границе «АОМ — раствор» (х = 0): 





ас! дФ (-7>)1,(0, у) 
И + С ЕЕ 0, р ЗЕ ь 40 
Ох т Ох } У) 52 ) 
С> (0, у)= М. (41) 
Ф(0,у)= 0, (42) 
21 (0,5) =0 (43) 
ОХ 
На границе «раствор — КОМ» (х = |): 
С(Ь У) = №, (44) 
9С> 0Ф (7: )7,,(, у) 
А 45 
Ох 622 ОХ } У) 5205 ) 
— О" +; р ОСТ +2 р ОС 
дФ ду ах “^^ а 
(= о (, >) (46) 
2 2 
Ох &1 ВАС1 + 5262С> 
21, у)=0. (47) 
ОХ 
На входе в канал (у = 0): 
С: (х0)=1 #=12, (48) 
МТ [рр ры 49) 
ду = +5505 ду ду 
1(х.0)= 0. (50) 
На выходе из канала (у = [): 
- р 5;С; 5 ый = 0, = 1,2, (51) 
ду ду 
м = 0, (52) 
ду 
1(х,1) =-, у = с0п5 . (53) 


После численного расчета системы (34), (35), (37)—(53) скачок потенциала Аф в канале обессоливания 


определяется по формуле: 
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й 17 
Ар = | ®(ьу)4у. (54) 
0 


Численное решение найдено методом конечных элементов с помощью пакета Соп1$01 Ми!арВу$1с$ на 
неравномерной вычислительной сетке (плотность элементов сетки увеличена у границ «раствор — мембра- 
на») [27]. 


Результаты исследования. Вычисления проведены для &=1,9.10”, Ре= 2355 , что соответствует 


следующим значениям параметров системы: 
— входная концентрация раствора электролита МаС Си = 0,1 моль/М’; 
— температура 7 `= 298 К; 
— коэффициенты диффузии катионов и анионов соответственно О; = 1,33:10`° м”/с, 2» = 2,05:10° м”/с; 
— числа переноса противоионов в мембранах Т! = 0,972, Т. = 1; 
— зарядовые числа ионов 1 = 1, > =-1[; 
— отношение концентрации противоионов на границе с мембранами к ее значению на входе в канал № = М, = 1; 
— й= 10 м — ширина канала: 
— 1=2.10 м — длина канала; 
— И = 3, 8-10 м/с — скорость прокачки раствора. 
На рис. 2 приведены поля концентраций С! и С›, потенциала ф и функции 1, рассчитанные при плотно- 
сти тока („= 1,5йт» где пм — это предельная плотность тока, определенная по формуле Левека (в безразмерной 
форме) [28]: 


1/3 


2 
. 24 =0,2. |. (55) 


т — › 7 
1 -& 





Здесь й = 0,395 — число переноса катионов в растворе [9]. 
й 








в) г) 


Рис. 2. Поля концентрации катионов С! (а) и анионов С. (6), потенциала ф (в) и функции т (г). 


Расчет по модели для гальваностатического режима при плотности тока й, =Ъ5й„ 
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В канале действует вынужденное течение, поэтому в областях у мембран истощение концентрации 
ионов увеличивается при удалении от входа в канал (вдоль направления, тангенциального к поверхностям мем- 
бран). Соответственно, толщина расширенной области пространственного заряда увеличивается вдоль кана- 








ла (рис. 3). 
01 
$ \2 1 
® 
5 0,05 
С. 
ыы 
х 0,98 0,99 1,00 
х 
а) 6) 
0,01 
= 0,005 
О 
0,96 0,98 100 


ви 
в) 
Рис. 3. Результаты численного расчета по модели для гальваностатического (сплошные линии) и потенциостатического 


(пунктирные линии) режимов при плотности тока Г, =1,51„, в сечениях у = 0, (1), у=0,41 (2), у=0,91(3): 


концентрационные профили С! и С. (а); увеличение рис. 3, а (6); плотность пространственного заряда р = <1С] + <2С. (в) 


Вольтамперная характеристика (ВАХ, кривая | на рис. 4) рассчитана на основе гальваностатической 


модели. При этом задаваемая плотность тока изменялась: („= 0; 0,0051; ...; 1,ЭЙрт. 


40 
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д к. 
— 20 
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== 


О 10 20 30 40 
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Рис. 4. ВАХ, рассчитанные по моделям для гальваностатического (1) и потенциостатического (2) режимов; 
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Предельная плотность тока определена по точке пересечения касательных к ВАХ в начальной части и 
на наклонном плато [5]. Здесь она совпадает со значением й„, оцененным по приближенной формуле Леве- 
ка (55). Также ВАХ была рассчитана на основе потенциостатической модели при изменении скачка потенциала 
Дф = 0; 0,4; ...; 40 (кривая 2 на рис. 4). Как видно из рис. 4, ВАХ Ти 2 совпадают. Таким образом, существует 
однозначное соответствие каждому скачку потенциала некоторой плотности тока, и наоборот. 

Рис. 3 и4 демонстрируют достаточно хорошее совпадение различных физико-химических характери- 
стик переноса, рассчитанных по моделям для гальваностатического и потенциостатических режимов. Это дока- 
зывает адекватность построенной авторами модели переноса в гальваностатическом режиме. 

Вычислительные затраты предлагаемой модели оценивались следующим образом. Фиксировалось вре- 
мя, затраченное на решение краевых задач для гальваностатического и потенциостатического режимов. При 


этом значения задаваемой точности вычислений, параметры системы были одинаковыми, а 1, = 1,5 Г„. Затем 


показатели времени сравнивались. Таким образом выяснилось, что расчет по гальваностатической модели тре- 
бует в 1,6 раз больше времени. Это связано с тем, что: 
— гальваностатическая модель содержит дополнительное уравнение для определения распределения плотности 
тока; 
— для потенциала на границе х = 1 установлено условие второго рода (46). 

Комплексный потенциал электромембранной системы для гальваностатического режима. В рабо- 
те [25] предложено обобщение импеданса электрохимической системы с помощью функции \, введенной в 


условиях электронейтральности. Аналогичное обобщение приемлемо и в данном случае. Отметим, что при ис- 
пользовании электрохимического импеданса объект рассматривается только как «черный ящик», и его внут- 
ренние свойства определяются косвенно. Полученные выше результаты позволяют ввести понятие комплексно- 
го потенциала электромембранной системы: Р=Фф-+1.1. Комплексный потенциал — это функция координаты 


любой точки внутри объекта, поэтому, в отличие от электрохимического импеданса, он позволяет исследовать 
внутренние свойства объекта. 

Обсуждение и заключения. Описан новый метод математического моделирования стационарного 
процесса массопереноса в гальваностатическом режиме для мембранных систем. При этом рассмотрен двумер- 
ный случай с использованием специального граничного условия, позволяющего задать плотность тока в систе- 
ме. Представлены уравнения для электрической функции тока. Результаты численного решения по потенцио- 
статической и гальваностатической моделям хорошо согласуются. Это показывает адекватность предлагаемой 
модели переноса в гальваностатическом режиме. 

Разработанная модель позволяет интерпретировать результаты экспериментальных исследований пере- 
носа ионов в мембранных системах, если данный процесс протекает в гальваностатическом режиме. Некоторые 
электрокинетические процессы связаны с появлением расширенной области пространственного заряда при 
сверхпредельных токах. Описывая формирование указанной области, можно выяснить, каким образом завися- 
щие от нее процессы влияют на перенос ионов в гальваностатическом режиме. 
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1.3 Роп $(ае Тесытиса1 Отуегзиу, Козюу-оп-Роп, Кизз1ап Недегайоп 


Введение. Рассмотрен частный случай синтеза одной и той 
же электромеханической системы управления методом 
максимума Понтрягина и методом синергетического 
синтеза. В качестве задачи была определена задача 
синтеза оптимальной по быстродействию 
электромеханической системы позиционирования, при 
этом момент сопротивления движению линейно зависел 
от выходной координаты системы. Этот подход к выбору 
задачи синтеза был обусловлен тем, что синтез 
оптимальных по быстродействию систем является одной 
из самых широко распространенных задач, которая 
решается при повышении эффективности действующих 
систем управления. 
Материалы и методы. Синтез оптимальной по 
быстродействию линейной системы управления на основе 
принципа максимума — широко распространенная задача 
в современной теории управления. Однако процедура 
синергетического синтеза такой формализации не имеет. 
Исходя из этого, в статье предложен подход, солижающий 
эти два метода, который, по мнению авторов, позволит 
повысить эффективность метода  синергетического 
синтеза, добавив в него некоторые особенности 
методологии синтеза оптимальных систем. 

Результаты исследования. В работе сформулированы два 
основных научных положения. Первое - применение 
принципа максимума для объекта класса двигатель 
постоянного тока при синтезе алгоритма 
позиционирования в условиях линейной нагрузки, 
функционально зависящей от угла поворота двигателя, 
позволяет оптимизировать систему по быстродействию. 
Второе - синтез системы управления на основе 
синергетического подхода позволяет получить систему, 
близкую к оптимальной (квазиоптимальную), но уже 
после модификации самого метода синергетического 
синтеза. Сформулирована гипотеза о возможной связи 
процедуры 
методе 


между 
расширения 


ВВОДИМЫМИ, при реализации 


пространства — состояния в 
синергетического синтеза, постоянными времени с 


определяемым В методе максимума оптимальным 
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Введение. В 20-м столетии потребности техники, в частности космической, выдвинули серию задач, 
для решения которых была разработана новая теория, получившая название теории оптимального управления 
[1]. Одним из основных методов при синтезе оптимальных систем управления является разработанный в 50— 
60 гг. двадцатого века метод максимума, над которым работали советские математики Л. С. Понтрягин и его 
ученики [2]. Применение этого метода основано на формализации задачи синтеза с переходом к форме задачи 
Майера и последующем решении систем линейных или, в некоторых особых случаях, нелинейных 
дифференциальных уравнений [1]. Стоит отметить, что необходимость решения систем дифференциальных 
уравнений, а в нелинейном случае это не всегда возможно, во многом является ограничением, как самого 
метода максимума, так и всей концепции синтеза оптимальных систем управления. 

В последние десятилетия широкое распространение получил новый подход к синтезу систем 
управления, который в своей основе опирается на синергетическую концепцию анализа и синтеза систем. 
Используемый метод, автором которого является А. А. Колесников, получил название метода аналитического 
конструирования агрегированных регуляторов (АКАР) [4—5]. В основе предложенного подхода лежит 
концепция синтеза нелинейных обратных связей. Они обеспечивают асимптотическую устойчивость системы 
управления по отношению к требуемому движению аттрактора в пространстве состояний [4]. Отличием этого 
метода от методов синтеза оптимальных систем управления является то, что здесь отсутствует как критерий 
оптимизации системы управления, так и констатация оптимальности синтезированного управления. С точки 
зрения реализации процедуры синтеза, метод АКАР обладает несомненным преимуществом перед методами 
синтеза оптимальных систем, которое выражается в отсутствии ограничений на нелинейность системы 
дифференциальных уравнений [4—5]. 

С практической точки зрения синтез управления процессами или объектами должен обладать 
возможностью ответа на вопрос, а есть ли другое управление, которое обладает свойством превосходства над 
всеми другими. Так современное массовое производство, в целях обеспечения конкуренции на рынке, 
постоянно требует повышения эффективности. Одним из наиболее востребованных способов такого 
повышения эффективности является оптимизация процессов управления. Этот подход применим как к 
системам автоматизированной сборки оборудования [6-7], так и к металлорежущим системам, в частности к 
системам управления сверлением в металлорежущих станках [8—15]. Идея совмещения этих подходов к синтезу 
систем управления обладает, с точки зрения авторов, несомненной практической ценностью. С научной точки 
зрения актуальным является совмещение преимуществ метода АКАР, которые выражаются в возможности 
учета при синтезе сложной нелинейной динамики процессов в машиностроении [16—19], с четкой и ясной 
формализацией постановки задачи синтеза и оценки ее достижения в методе максимума [1—2]. 

1. Синтез базовой математической модели и постановка задачи исследования 

В современной экономике широкое распространение получило направление синтеза оптимальных по 
быстродействию систем, которое позволяет добиться существенного роста эффективности систем в целом [7— 
8]. Так, при решении задач автоматизации процессов сборки различного оборудования, часто возникает задача 
крепления различного вида частей друг к другу [6-7]. Здесь от скорости проведения этой операции зависит 
экономическая эффективность всего процесса производства. Эта же ситуация наблюдается в системах 
обработки металлов резанием на металлорежущих станках [9-19], в которых чем быстрее осуществляется 
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процесс обработки деталей, тем меньше издержки всего производственного процесса. Исходя из рассуждений, 
примем в качестве допущения, что требуется синтезировать оптимальную по быстродействию систему 
управления, которая с учетом схожести процессов автоматизированного крепления деталей и сверления 
глубоких отверстий может быть проиллюстрирована следующей схемой (рис. 1). 


С 





дпт =) 


Рис. 1. Иллюстрационная схема 











На рис. 1 представлен двигатель постоянного тока, осуществляющий либо операцию крепления 
оборудования, либо подачу в сверлильном станке. В случае операции крепления — С суммарная 
характеристика жесткости сжимающегося при скручивании материала. Для случая сверления — С 
характеризует линейное нарастание сопротивления подаче при накоплении стружки в стружкоотводящих 
каналах сверла. 

В основе исполнительной подсистемы системы управления в обоих рассматриваемых случаях лежит 
электропривод, обеспечивающий преобразование электрической в механическую энергию движения 
исполнительных механизмов системы. Примем, что в рассматриваемом случае такое преобразование 
осуществляется электродвигателем постоянного тока с коллекторным управлением, который описывается 


следующей системой уравнений [20]: 


ео а и 


И 
сС=У ее +М, 
ти 
где И — напряжение, подаваемое на коллектор двигателя; г — ток, потребляемый двигателем; К, [ — 
параметры электрической части двигателя; / — параметр, характеризующий инерционные свойства ротора 
двигателя, приведенный инерционный момент всех вращающихся масс; х — частота вращения ротора 
двигателя; М, — внешний приложенный момент сопротивления; с„, с. — механическая и электрическая 


постоянные двигателя. В рассматриваемом случае момент представлен как функция угла поворота ротора 


двигателя, то есть М, = Со. С учетом этого (1) перепишем как: 


ооо 8 














р ИИ 
© 
сСй=У—+Са . (2) 
И 
(о 
— = 
ИИ 
Выразив во втором уравнении значение тока и подставив полученное в первое уравнение, получим: 
1 [7 45 ВЛ ах 1© КС 
И—-®= - р 04 
Г. сс, @ с.с, @ СС, Со. (3) 
44 
— = 
(А 
Разрешив первое уравнение относительно старшей производной и переписав его вторым, получим: 
44 
щ “ 
| 4 
Ио Ка®х сс, С с, (9) 
Е —)®-———@9+— И 
И га Ш Л 
46 
Обозначим переменные и=хХ,®=Х,, а - ре В качестве постоянных введем 
1 
К ге. ©. &С Е 
— = а, [м +—)=а,,—=а., —“ =Б. Тогда система (4) примет вид: 
Л Л ЕЛ ЕЛ 
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и 

и“ 
> 

Я - (5) 
1х 

> =-ах -а»х, -азх,+БИ 


Или в матрично-векторном виде: 


х=Ах+ Ви, (6) 
0 1 0 0 рр 0 
где А=| 0 0 1 |, В= О], х=х, | и=|0 
Ч: Ч» 4 р хе й] 


Примем в качестве привода, обеспечивающего движение рассматриваемой системы, электродвигатель 
постоянного тока, параметры которого обеспечивают следующие значения постоянных: аз1=4,65, аз›=4,6, азз=2. 


0 1 0 
Матрица А, с учетом этого, примет вид: А=| 0 0 1 |. Зная матрицу А, найдем собственные числа /^, 
—,65 —,6 —2 


матрицы А: Л = —1,2666 + 0,00001, ^› =-—0,3667 + 1,8809, /з = —0,3667 — 1,88061. 

Как видно из полученных собственных значений матрицы, система управления асимптотически 
устойчива по Ляпунову [21]. Более того, можно утверждать о колебательном характере процессов протекающих 
в системе, так как собственные значения не только содержат отрицательные вещественные части, но имеют 
ненулевые мнимые части. 

2. Синтез управления методом максимума Понтрягина 

Для начала сформулируем задачу синтеза оптимального управления в следующем виде: 

—щ используя принцип максимума для объекта, описываемого Системой (5), определить алгоритм 
оптимального уравнения, который обеспечивает перевод объекта из начального состояния х/(0)=50, х2(0)=0, 
хз(0)=0 в конечное состояние х(Т)=0, х>(Т)=0, х(Т)=О за минимальное время Т. При этом необходимо 
определить число переключений, моменты переключений и построить кривые управления и(Т) и координат х/(1), 
х2(1), Хз(1). 

На управляющее воздействие наложено ограничение |и<Иук. Параметры системы: 
а =Ъа, =12, а, = 20,6 =1; Окк=440 В. 


Р ешение: Запишем Ггамильтониан: 
Н=У чу, (7) 


где у; и } — координаты векторов уи /. В дополнение к системе (5), составим систему уравнений для 


вспомогательных переменных \у, ‚, где 1=1,2 , используя следующую зависимость: 





-, 1=1,2. 


Или в раскрытом виде: 








(5) 
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С учетом (5), система уравнения (8) примет вид: 











а$ 
-=4,65$. 
а 
г =—$, +4,6Ф, (9) 
а$ 
Е -=-ЧФ, +24, 
1 
Матрица собственных значений для случая (9): 
1,2666 0 0 
А"? = 0 0,3667 + 1,8806 0 
0 0 0,3667 + 1,8806 
Решение системы (9) для диагонализированного случая матрицы А%-? примет вид: 
уе = Фа 
р — С або (10) 
2 2. . 


р =С ао 
3—3 


Г. учетом матрицы уз ‚В ИСХОДНОМ базисе решение примет ВИЛ: 
4, =-0,9474С 61255 + (0,5812 —0, 3347 Себе + (0,5812 0,3347) С,6 3657188050! 
м7 — = | 1893С в 700 +0 6883С ао _0 6883С вов 
2 2 1 2 и) › 3 , 
Ч, =-0,2581С. 62589 + (0,1812+0,20860 Сер + (0,1812 0,2086) Се“? вр, 


где Су, (>, С, — постоянные интегрирования. 
Общее выражение, описывающее функцию Гамильтона: 
ПАР Зе (11) 

В выражении (11) важную, с точки зрения метода синтеза, роль играет член, включающий в себя 
управление: 

Н' =4.0 =(—-С.0,2581е'°*' + С. (0,1812-+0,20867 + С, (0,1812—0,2086е°°°°7` 9). 

Для того, чтобы гамильтониан Н, вычисляемый по формуле (11), принимал максимальное 
положительное значение, слагаемое Н” должно быть всегда положительным и наибольшим. Для этого алгоритм 
оптимального управления должен имеет вид: и(1)=сИмах, где 

с = 5191(—0,2581Се'*°%"' + (0,1812+0,2086)е°°°°"*'°6Р! + (0,1812 —0,2086)е°°°7`'°65г). 

Оптимальное управление и(Г) является кусочно-постоянной функцией, принимающей значения = Онах, И 
имеет не более двух интервалов постоянства, так как нелинейная функция 
—0,2581С ее" + (0,1812+0,20860е°°°'*°°! + (0,1812—0,2086е°°°7”`88 6" , 

меняет знак не более одного раза. При этом возможная смена знака происходит с плюса на минус, то есть для 
выполнения принципа максимума требуется сначала подавать на двигатель Ч= +0, а затем 9= —Энах. 
Проверим эти рассуждения, построив полученную функциональную зависимость 4’, в пакете МаЙаь (рис. 2). 





Рис. 2. График изменения \Р. 
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Как видно из рис. 2, Ч, действительно меняет знак с (+) на (—) только один раз. С учетом этих 


рассуждений определим алгоритм оптимального уравнения, который обеспечивает перевод объекта из 
начального состояния х1(0)=0, х›(0)=0, хз.(0)=0 в конечное состояние х!(Т)=50, х›(Т)=0, х.«(Т)=0 за минимальное 
время Т. Здесь отметим, что в современных условиях нет необходимости получать аналитическое решение 
исходной системы уравнений. Можно, используя имеющиеся возможности современных пакетов прикладных 
программ типа МаЧаб, достаточно просто и наглядно получить численное решение для рассматриваемого 
случая. 

Результаты численного моделирования для расчетного оптимального по быстродействию алгоритма 
управления представлены на рис. 3. Параметры модели подобраны таким образом, чтобы она выполняла 
требуемое по правой границе граничное условие. 


Х!, рад 
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Рис. 3. Графики выходной координаты (а) и переключения управления при # = 0,5 (Ъ) 


Как видно из рис. 3, время переключения управления с (+ вх) на (-Овах) выбрано как 1 =0,5 секунды, 


а общее время регулирования составило и Т=|1,3 секунды. При этом система управления выходит на 


требуемый нам уровень выходного сигнала, то есть достигается правая граница краевых условий. 

3. Синтез системы управления методом АКАР 

Для синтеза системы управления методом АКАР можно не осуществлять переход к абстрактному 
случаю пространства состояния (см. (5)), а удобно использовать исходную систему (2). Однако требуется 
преобразовать ее к следующему виду: 


2 
—_ =-а хм таз: > (12) 


ПМ) 


3 —_—_ — — 
——=-фа»х, -а,х, БИ 


ПМ) 


[@] С К 1 


С 
где переменные &«=х,®=х,,1=х,; в качестве постоянных введены —“ =а.,—=а,—=а, —=а.,—=Ь. С 
У `Ё, Г. 


учетом определенных ранее значений, получим: а, =1,515, а», = 2,325, а», = 2,6 =1 и система (12) 


__ 
ь 
з 
| 


примет вид: 
Ах, 
р вы” 
4% 2 325х ВО (13) 


1 
Ч 
Е О 
ИИ 
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Матрица А, с учетом этого, примет вид: 


0 1 0 
А=| -2,325 0 1,515 
0 —15 —2 
Зная матрицу А, найдем собственные числа /, матрицы А: Л =-1,2674 + 0,00001; Л› = —0,3663 + 1,88011; 
Лз = -0,3663 — 1,88011. Как видно из этих результатов, речь идет именно о рассматриваемом при синтезе 


методом максимума случае системы управления. 
Координата, характеризующая угол поворота ротора (х!), является выходной координатой системы. 
Поэтому для формирования требований к желаемому поведению системы в пространстве состояний введем 


макропеременную порядка т. =х —х =>0, где х„ — заданное требуемое значение координаты х,. Значение 
скорости изменения угла должно обеспечивать стремление угла поворота ротора к заданному выше значению. 
Для этого введем еще одну макропеременную Ч, =х, —$,(х)>0, где ф,(х) — некоторая функция, 
описывающая стремление х, > ф,(х/) в стационарном состоянии. Дальше вводим макропеременную порядка 
по координате Ч. =х, —-$.(х,х,) > 0, где ф.(х,х,) — некоторая функция, описывающая связь между 
координатами в стационарном состоянии системы и, следовательно, х, = ф.(х,^х,). 

Для вновь введенных макропеременных Ч,, Ф,, ЧФ, потребуем асимптотически устойчивый закон 


изменения, то есть выполнение следующей системы дифференциальных уравнений: 








аф 
т +4 =0 
д 
р т 4, =0, (14) 
и 
ТА +У, =0 
М 
где Г, Г, ТГ, — постоянные времени, определяющие скорость стремления введенных макропеременных к 
| 2 3 р р р р ропер 


нулю, либо, иными словами, скорость стягивания произвольной траектории движения системы в пространстве 


состояния к требуемому аттрактору. 

Введенная в выражении (14) система уравнений расширила пространство состояний системы с 3-го до 
6-го, что является одним из основных положений метода АКАР. Дальнейший синтез системы управления 
сводится к поэтапному процессу декомпозиции (сжатия) этого пространства до начального уровня [1-2]. 
Начнем процесс декомпозиции с подстановки в первое уравнение системы (14) вместо макропеременной Ч! ее 
значения в координатах управляемого процесса. В результате получим: 


Ах 
Паб ж) =0, (15) 
С учетом (13), уравнение (15) примет вид: 
1х +2 -х)=0. (16) 


Учитывая следующий шаг алгоритма синтеза, в результате которого приняли, что х => ф,(х), найдем 


значение ф,(х), как: 


И) 
фо = > о (17) 
1 
Учитывая это, макропеременная Ч. примет вид: 
Ф, ое —\. (18) 
1 


Справедливость выражения (17), с точки зрения целей управления, подтверждает тот факт, что в стационарном 
состоянии значение х, > ф,(х) > 0. Тогда, с учетом (18), второе уравнение системы (14), примет следующий 


вид: 


Хх аф. (х 
Пт, к, $.) =0. (19) 


Или, с учетом значения ф,(х): 
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0. 20 
и Та Т —_ 


да: 
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На следующем шаге алгоритма синтеза, пользуясь теми же рассуждениями, что и раньше, и тем 
фактом, что х, => ф.(х,х,) , определим из (20), с учетом (13), значение ф.(х,х,) в координатах управляемого 
процесса: 


Г -+Т 
5.) 535% —0,66 212, 0,66 


1 21 


Получив значение ф.(х,х,) в координатах управляемого процесса, можно определить значение 


№0. (21) 


макропеременной Ч. в координатах состояния системы и раскрыть третье уравнение системы (14). 


4х. _ ЧФ. (х,х,) 
И 


т 1+х, -$.(х,х,)=0. (22) 


После подстановки в (22) полученных ранее значений ф.(х,х,) из (21) и использования системы (13) 


на этом шаге декомпозиции пространства состояния системы управления, определим значение управления в 
координатах управляемого процесса: 


ЕЁ. АРЬЕ Е Е 
И =1,5х, +2х, +1,535х, ЕЕ Я ВИА 
1151 1151 


со, ббх, ПЕТЬ 0,66 
17 2-3 17 2-3 
Выражение (23) и определяет асимптотически устойчивое управляющее воздействие на систему 
управления, описываемую уравнениями (13). 


Результаты моделирования системы дифференциальных уравнений (13) с требуемым значением угла 


(23) 


поворота ротора двигателя х!0-50 рад и значениями введенных постоянных времени Г =0,1, Г, =0,2, 


Т, =0,3 представлены на рис. 4. 
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Рис. 4. Результаты моделирования с управлением при Г =0,1, Т, =0,2, Т; =0,3: 


переходный процесс по выходной координате (а); график изменения управления (Ъ) 


Как видно из рис. 4, в сравнении с рис. 3, структура управления в обоих случаях одинакова, то есть 
речь идет о том, что в начале подается положительное по знаку управление, а затем его знак меняется на 
противоположный. Время регулирования резко увеличивается до значения близкого к 2,5 секунде, но при этом 
максимальное управление достигает значений больших 5 000 В. Это недопустимо с точки зрения метода 
максимума. Для ограничения максимального значения управления введем в контур управления нелинейное 
звено типа «насыщения» с таким же порогом, как и в случае максимума. Результаты моделирования системы с 


учетом этого порога и с подобранными значениями постоянных 1, Т,, Г, представлены на рис. 5. 
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Рис. 5. Результаты моделирования с управлением при 1, = 0,15, Г, = 0,16, Г, =0,17: 


переходный процесс по выходной координате (а); график изменения управления (Ъ) 


Как видно из рис. 5, введение нелинейного ограничения на управление типа «насыщения» позволяет 
получить требуемый с точки зрения допустимости управления результат. При этом фактически соблюдается 
определенная методом максимума структура управления, но управление не является оптимальным по 
быстродействию. Здесь необходимо отметить, что такая модификация метода АКАР может приводить и к 


потере устойчивости системы (рис. 6). 
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Рис. 6. Результаты моделирования с управлением при 1, =0,015, Т, =0,016, Т, =0,017: 


переходный процесс по выходной координате (а); график изменения управления (Ъ) 
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Однако, несмотря на результат, представленный на рис. 6, дополнение системы управления полученной 
методом АКАР, ограничениями на амплитуду управляющих сигналов существенно сближает ее с решениями, 
получаемыми при реализации процедуры синтеза методом максимума. 

Обсуждение и заключения. Таким образом, для базовой модели системы управления (2) было 
синтезировано оптимальное по быстродействию управление методом максимума. Для этого же случая, в 
третьей части работы с этими же параметрами системы и граничными условиями было синтезировано 
управление методом АКАР. Оба эти управления внешне несхожи. В случае метода максимума — это 
программное кусочно-постоянное управление, где важен оптимальный момент времени переключения знака 
управляющего сигнала. В случае АКАР — это непрерывная гладкая функциональная зависимость управления 
от значений координат состояния системы. Однако, при их реализации видна структурная схожесть, 
выражаемая в смене знака управления. После модификации метода АКАР структурная схожесть стала еще 
более видна. Отметим, что добиться параметрического подобия реализуемых управлений можно за счет 
подбора вводимых при синтезе методом АКАР постоянных времени. Это позволяет авторам сформулировать 
следующую гипотезу. 

— Подбор вводимых в процедуре синтеза методом АКАР значений постоянных времени Т,Т.,Т, ‚в 


случае модификации полученного управления способом предложенным в статье, позволяет добиться от 
синергетического управления свойства оптимальности по быстродействию. Авторам не удалось получить 
общее доказательство этой интуитивно понятной гипотезы. Однако такая задача не стояла в рамках этой 
работы. В результате численных экспериментов с вариацией Т,,Т,,Т, удалось добиться квазиоптимальности 


управления, полученного методом АКАР. 
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Параллельное построение двоичного дерева на основе сортировки" 


Я. Е. Ромм', Д. А. Чабанюк”” 


`? Таганрогский институт имени А. П. Чехова (филиал) Ростовского государственного экономического университета (РИНХ)", г. Таганрог, 


Российская Федерация 


РагаПе! сопзегисНоп оЁ Ытагу гее Базе оп зогии?”” 


Уа. Е. Вопит', О. А. Свабапуик”” 


`? Тазапгое Свеквоу шзийке, Возюу Зе ОшуегзИу оЁ Есопопис$ (ВПМН) Бгапсв, Возюу-оп-Роп, Визаи Еедеганоп 


Введение. Разработаны алгоритмы параллельного постро- 
ения двоичного дерева. Алгоритмы выполнены на основе 
сортировки и описаны в конструктивной форме. Для мно- 


№ 


жества из элементов временная сложность имеет оцен- 


Т(Ю=О(Ф и Т(Ю)= Оо, №) 


КИ где число процессо- 


2 
в Дерево строится со свойством един- 


ро 
ственности. Алгоритмы инвариантны относительно вида 
входной последовательности. Целью работы являлась 
разработка и исследование способов ускорения процесса 
организации и преобразований древовидных структур 
данных на основе алгоритмов устойчивой максимально 
параллельной сортировки для их применения к базовым 
операциям информационного поиска в базах данных. 
Материалы и методы. Взаимно однозначное соответ- 
ствие множества входных элементов и построенного для 
него двоичного дерева устанавливается при помощи 
устойчивой адресной сортировки. Сортировка обладает 
максимальным параллелизмом, в операторной форме 
устанавливает взаимно однозначное соответствие вход- 
ных и выходных индексов. На этой основе разрабатыва- 
ются методы взаимного преобразования двоичных струк- 
тур данных. 

Результаты исследования. Получен эффективный парал- 
лельный алгоритм построения двоичного дерева на основе 
адресной сортировки с 


временной сложностью 


Т(№) = Оо, №). От известных аналогов алгоритм отли- 


чается структурой и логарифмической оценкой временной 
сложности, позволяющей достигать ускорения аналогов 


порядка О(№), а>1. В качестве усовершенствованного 


варианта предложена модификация алгоритма, обеспечи- 
вающая максимально параллельное построение двоичного 
дерева на основе устойчивой адресной сортировки и 
априорного вычисления хранимых индексов корней под- 
деревьев. Алгоритм отличается структурой и оценкой 
временной сложности Т(1)=0О(). Аналогичная оценка 


достигается в последовательном варианте модифициро- 


“Работа выполнена в рамках инициативной НИР. 
" Е-тай: гот @ Н6.ги, депсБафапуйк @ отай.сот 
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Обсуждения и заключение. Полученные результаты 
направлены на организацию эффективных способов дина- 
мической обработки баз данных. Предложенные способы 
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Введение. В области современных высокопроизводительных вычислений имеется тенденция к конверген- 
ции технологий параллельной обработки информации и различных архитектур процессоров. Несмотря на многооб- 
разие архитектур процессоров и способов представления информации, для повышения скорости обработки 
данных одной из важнейших задач информатики является идея параллельной обработки. Для ускорения обра- 
ботки данных авторы предлагают использовать алгоритм устойчивой адресной сортировки, обладающей мак- 
симальным параллелизмом. 


Метод параллельного построения двоичного дерева. Для массива А = ( ба а, матрица сравнений 


строится в соответствии с [1, 2]. Элемент а,, этой матрицы определяется как 


+ @, а: 
а,, = 121 (а, а, = 0, а, =а, 
=: бо Я, 
) 
Де = 2: И 
п 
Элемент а, в отсортированном массиве С - в. Ре о получает номер К = > а ‚ где а,, 20 
1=1 


при 1 <], а,, >0 при 1> 7. Все сравнения взаимно независимы, сортировка устойчива и максимально парал- 


2 
лельна с оценкой временной сложности Т о = О(1) . На этой основе можно выполнить параллельное 


построение двоичного дерева [3, 4]. Пусть дано множество из М№ элементов Х, ‚, все элементы которого пред- 
1 


ставлены в виде одномерного массива. На множестве предполагается заданным отношение порядка <. Требу- 


ется преобразовать массив в двоичное дерево. Для этого выполняется описанная сортировка массива. Середин- 


ь | № 
ный элемент массива С имеет индекс }., = Е и принимается за корень дерева [3]. Все элементы массива С 


слева от С, образуют левое поддерево (левый подмассив). Элементы справа от С, образуют правое поддере- 
ср ср 


во (правый подмассив). Левый подмассив интерпретируется как новый массив. В нем аналогично находится 


индекс корня ] нь 2 = - — — — 5 ‚ При этом С, „„— ближайший слева потомок корня дерева 

С т Все элементы подмассива слева от С ‚ .› Н@ превосходят С, ‚> Все элементы подмассива справа не 
ср ср. лев. ср. лев. 

меньше Одновременно определяется индекс корня правого подмассива 
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| 1 
МТМ |=, + Е 


Лоо прави: = 5 т я |= . При этом С. — ближайший справа потомок корня дерева 


Уер. прав. 1/2. 


С, . Процесс рекуррентно возобновляется в каждой паре прилегающих подмассивов: 
ср 


: __ Че лев. 1/21 1 
ыы лев. 1/2', 1 = 2 2 
ы —_ Е? ны лев. 1/27! — 1 
Тор. лев. 1/2',2 — ВВ, о 2 › 
; г. 1 ср. прав. 1/21 | Л ср. прав. 1/2”? —1 
Лео. прав. 1/2',1 — 9, прав. 1/27 2 ) 
7 Еф ^^ ] Ри) —1 
7 гы: ср. прав. 1/2 ср. прав. 1/2 г. 
Лор. прав. в — НА, о 2 , № р веаь 105, №. 


В результате за время О(1) формируются все элементы нижестоящего уровня двоичного дерева. Про- 


цесс можно продолжать до исчерпания 105, М уровней двоичного дерева. 


Число шагов алгоритма построения двоичного дерева в параллельной форме складывается из шага сор- 


тировки И последовательности шагов при расчете индексов корней поддеревьев. Отсюда 


Т(К)=105, Мя+т=0(108, №), где К — число процессорных элементов, т — время бинарного сравнения, 1 


— время вычисления одного индекса корня. Число процессоров А определяется максимально параллельной 


сортировкой Л входных элементов, а затем вычислением индексов с удвоением по числу уровней дерева. При 


> 1 М№М-—1 > 
вычислении индексов это число не превзойдет 2° =М№/2, поэтому числа процессоров, задействованных 


2 


о ы М 
сортировкой, достаточно. В итоге А не превзойдет — [3]. Окончательно, временная сложность парал- 


лельного алгоритма построения двоичного дерева составит 


2 
Т —— =0(10з, М). 


Пример [3]. Двоичное дерево ДлЯ массива из 15 элементов 
Х = (14, 9, 24, 7, 11, 20, 28, 3, 68, 10, 13, 17, 21 25, 30) строится следующим образом. 
Результатом сортировки является массив 


С = (3, 7, 8, 9, 10, 11, 13, 14, 17, 20, 21, 24, 25, 28, 30). 


15 


Корнем двоичного дерева является серединный элемент массива С: },= 5 


=8, С, =М. Левый 


| 8—1 ь 
подмассив имеет корень Л, лев. 12 = ет =4, элемент С, =9 — корень левого поддерева, который является 


ближайшим слева потомком  серединного элемента С, . Правый подмассив имеет корень 
ср 


8—1 


ие = В+ — =12 , элемент С, = 24 является корнем правого поддерева и ближайшим справа потом- 


С ыы. 2 С. =7 б И 
ком корня С, . Далее, ем: о: элемент С, = 7 — корень поддерева слева и является олижаи- 


шим слева потомком корня поддерева С . В правом поддереве корень имеет номер 


Уср. лев. 1/2 


4—1 о 
И =4-+ и =6, элемент С; =1| — корень правого поддерева и ближайший справа потомок 


Информатика, вычислительная техника и управление 


ты 
л 
= 


Вестник Донского государственного технического университета. 2016. Т. 18, №4. С. 449-454. 155М№ 1992-5980 е155М 1992-6006 
Уезник о} оп Уие Тесйтса! Иттуегуйу. 2018. Уч. 18, по. 4, рр. 449-454. 155М№ 1992-5980 е135М 1992-6006 





| 12-—8—1 
. Аналогично, слева от С, определяется корень Л.‚ прав. 1/41 =12— о =10, элемент С, =20 


Уср. лев. 1/2 Яр. прав. 1/2 


— корень левого от него поддерева и ближайший слева потомок корня поддерева С. . Для смежного с рас- 


7 ср. прав. 1/2 


12 — 14 28 
смотренным правым подмассивом корень имеет номер /о, прав. 114.2 217+ - =14 , элемент С = НЕ 


ближайший справа потомок С. и корень правого поддерева. Нижний уровень дерева сформируют потом- 


|, ср. прав. 1/2 


ки, оставшиеся слева и справа от каждого из 4 идентифицированных корней (рис. 1): 


= 8 


0 уровень 






Ур. лев. 1/2 — ыы Ур. прав. 1/2 =12 


| уровень 








ий) = 14 


Уср. прав. 1/4,2 





Уср. лев. 1/4,1 Уср. лев. 1/4,2 = Уср. прав. 1/4,1 — 












=7 И = 15 


И прав. 1/8,1 — ыы прав. 1/8,2 — Уср. прав. 1/8, 3 прав. 1/8,4 


3 


Ур. лев. 1/8,1 =1 Уср. лев. 1/8,2 Ур. лев. 1/8,3 =5 Ур. лев. 1/8,4 


Рис. 1. Пример построения двоичного дерева на основе сортировки 


Имеет место 
Теорема 1 [3]. Для одномерного массива из № элементов двоичное дерево может быть построено па- 


2 
раллельно при помощи сортировки с временной сложностью Т > О (105, №). 


Использованная сортировка устойчива, как следствие двоичное дерево строится с единственностью. 
Индексы всех серединных элементов (всех корней поддеревьев) можно идентифицировать [3]. С учетом этой 
модификации все индексы из предыдущего примера для № значений поддеревьев можно вычислить синхронно 
и взаимно независимо. Это приводит к единичной оценке времени построения двоичного дерева. Для каждого 
конкретного № все значения индексов узлов дерева можно вычислить априори и хранить в памяти компьюте- 
ра. С их помощью отсортированные элементы можно синхронно и взаимно независимо адресовать по всем уз- 
лам дерева. Формулы вычисления индексов узлов зависят только от общего количества № входных элементов 
и никак не зависят от их взаимного расположения после устойчивой сортировки. Для упрощения адресации 
памяти вычисленные индексы можно упорядочить на каждом уровне и расположить по возрастанию уровней. 
Тогда по ключу № считывается вся совокупность упорядоченных индексов узлов. Остается только по счита- 
ным адресам расположить отсортированные элементы дерева. На основании изложенного имеет место 


Теорема 2. Для одномерного массива из № элементов двоичное дерево может быть построено парал- 


2 
лельно при помощи сортировки и априорного вычисления индексов с временной сложностью Т - = О(1) 


Ниже представлена единая таблица, содержащая формальные оценки временной сложности последова- 
тельных и параллельных алгоритмов построения двоичного дерева в сопоставлении с предложенными алгорит- 


мами. 
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Таблица 1 
Сравнительные оценки временной сложности последовательных 
и параллельных алгоритмов построения двоичного дерева 
в сопоставлении с предложенными алгоритмами 


Алгоритм Временная Ускорение при использовании | Ускорение при использовании 
построения сложность алгоритма с единичной алгоритма с логарифмической 
двоичного дерева алгоритма временной сложностью временной сложностью 


А]еогифт о? Гагапа в 
А., Китаг У. (2004) Т=ОМ* | [5] 
[5] 
А] отт оЁ 
СраегиазооК Р. 
(2015) [6] 


Полиномиальный 
алгоритм (2016) [7] 


Алгоритм «ей св! 
— пор $1» (2014) 


[8] 


Райегп-Базед 
а]>огифит (1991) [9] [9] 


Представленный 
алгоритм с 
логарифмической Т=0(ю5, М) [3] 
оценкой временной 
сложности (2015) [3] 
Представленный 
алгоритм с 
единичной оценкой ТГ” = О(1) 
временной 
сложности (2015) [3] 


В таблице 1 Р — мощность словаря шаблонов, № — число входных элементов двоичного дерева, К 
— размерность пространства, в котором выполняется сортировка. 


Из таблицы видно, что предложенный алгоритм с логарифмической оценкой временной сложности аб- 
страктно улучшает оценки известных алгоритмов. Минимальное ускорение достигается по отношению к алго- 


— 


Г ТГ 
ритму из [5]: —- о(м . 12) ‚или, == О (№), а максимальное ускорение достигается относительно полиноми- 


Т № Т Е . 
ального алгоритма из [7]: = =О ИЛИ т = О(м | В случае с предложенным алгоритмом с единичной 





105, № 
оценкой временной сложности также улучшаются оценки известных алгоритмов. В данном случае минималь- 
Т 
ное ускорение достигается относительно алгоритма из [5]: И (№), максимальное ускорение достигается 


— 


Г 3 
относительно полиномиального алгоритма из [7]: т = О(м ) 





- = 

Заключение. Разработанные алгоритмы отличаются от известных способов [5—7, 10, 11] построения 
двоичного дерева тем, что используется максимально параллельная сортировка для вычисления индексов узлов. 
При этом для построения дерева либо затрачивается логарифмическое число шагов, либо вообще не затрачива- 
ется дополнительное время, если значения индексов априори рассчитаны для всех значений М№ в некоторых 
реальных границах и хранятся в памяти компьютера. Предложенный параллельный алгоритм построения дво- 
ичного дерева может использоваться с целью организации эффективных способов динамической обработки баз 
данных. 
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